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Final Report for NASA NAG3-2179 


Ian A. Waitz, S. P. Lukachko, and A. Chobot, Massachusetts Institute of Technology 
R. C. Miake-Lye and R. Brown, Aerodyne Research Incorporated 

Prompted by the needs of downstream plume-wake models, the Massuchusetts Institute of Technology 
(MIT) and Aerodyne Research Incorporated (ARI) initiated a collaborative effort, with funding from the 
NASA AEAP, to develop tools that would assist in understanding the fundamental drivers of chemical 
change within the intra-engine exhaust flow path. Efforts have been focused on the development of a 
modeling methodology that can adequately investigate the complex intra-engine environment. Over the 
history of this project, our research has increasingly pointed to the intra-engine environment as a possible 
site for important trace chemical activity. Modeling studies we have initiated for the turbine and exhaust 
nozzle have contributed several important capabilities to the atmospheric effects of aviation assessment. 
These include a more complete understanding of aerosol precursor production, improved initial 
conditions for plume-wake modeling studies, and a more comprehensive analysis of ground-based test 
cell and in-flight exhaust measurement data. In addition, establishing a physical understanding of 
important flow and chemical processes through computational investigations may eventually assist in the 
design of engines to reduce undesirable species. 

One key area of investigation has been into the intra-engine oxidation of S0 2 to SO 3 and H 2 S0 4 . The 
amount and nature of aircraft sulfur emissions that are deposited in the atmosphere (Fahey et al., 1 995; 
Curtius et al., 1 998) depend on how much fuel sulfur leaves the engine in a condensable form. The 
combustion process oxidizes fuel bound sulfur to mainly S0 2 , which will not condense, but post- 
combustor chemistry can further oxidize the S0 2 to SO 3 and H 2 S0 4 , S(VI), which can participate in 
condensation processes as the wet exhaust gases cool. These newly formed volatile aerosols then 
contribute to the atmospheric aerosol loading and can lead to various potential atmospheric impacts, 
including direct radiative forcing, contrail formation, and cirrus cloud perturbations. In order to 
understand the aerosol emissions from aircraft, the conversion of fuel sulfur to condensable sulfur species 
is necessary. Because of the difficulty and cost of engine measurements, we have developed modeling 
tools to simulate the post-combustor chemistry and gain insight into the conversion processes. A list is 
included in Appendix A of publications and major presentations; Appendix B contains a student thesis 
(Chobot). Jointly these publications and the thesis that serve as a more complete account of the MIT/ARI 
work described below. 
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One-dimensional models have been used to follow both homogeneous (Brown et al., 1996; Tremmel and 
Schumann, 1999) and heterogeneous (Brown et al., 2000) chemistry along a temperature/pressure 

trajectory through the turbine. Our effort has been to 
develop additional capabilities to extend such one- 
dimensional studies. In order to determine those regions in 
the turbine where trace chemistry is most active, and where 
more detailed multi-dimensional modeling is appropriate, a 
tool was developed to evaluate the potential for oxidation 
within trace species families. The analysis compares 
appropriately defined characteristic flow and chemical time 
scales, that, with knowledge of the thermodynamic 
reference state, can be used to understand conversion 
potential in pressure-temperature coordinates. Figure 1 
shows an example space for the conversion of S0 2 to SO 3 . 
The contour lines represent levels of formation through a 
single turbine blade row. The red dot places the inlet 
condition to that blade row, here for the first in the high- 
pressure turbine, on the temperature and pressure map. The 
temperature and pressure flow path through the entire 
engine is shown as the blue line. 

It is worth noting that the flow will encounter cooling as it passes through the turbine blades and the 
shaded region extending to lower temperatures from the high pressure turbine point on the trajectory 
indicates a range of temperatures associated with such cooling. Clearly, the effects of cooling move a 
given parcel of fluid into regions of higher conversion relative to the uncooled trajectory at the same 
pressure. Parallel analyses can be performed for different engine cycles, to examine the effects of engine 
cycle parameters, or for subsets of the complete chemical mechanism, to evaluate the relative 
contributions of identified chemical pathways. 

Phenomena or locations in the engine to which trace chemistry is found to be sensitive can be investigated 
with higher-dimensional tools. For detailed investigations of the impacts of turbomachinery fluid 
mechanics, we employ a CFD tool developed for this program called CNEWT. CNEWT combines a well- 
established 3-D turbomachinery code and a chemical kinetics solver within a structure currently capable 
of calculating passively reacting internal flows where chemistry has no influence on the course of the 



Figure 1. Characteristic Time 
Scale Analysis for SO, Chemistry 
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flow solution. An additional 1-D kinetics capability, used for solving along averaged profiles of relevant 
flow parameters, extends the CNEWT code to allow approximations at several levels of model detail. 
Currently, CNEWT is best used for understanding trends and influential fluid and chemical processes 
rather than as a predictive tool since empirical data does not presently exist to validate CNEWT for the 
applications of interest. 



Figure 2. 1-D Parametric Analyses. Cases shown for an advanced subsonic engine cycle developed 
by NASA. Baseline case (black lines) initiated at a combustor exit temperature of -1600 K and 
-15 atm. Comparative case (red lines) shown for combustor exit temperature of 1390 K. 

One-dimensional calculations performed with this tool constituted the first broad parametric 
investigations of chemical processes in the turbine and exhaust nozzle. Using a chemical mechanism 
developed through the work of Robert C. Brown of ARI and Fred L. Dryer of Princeton University, 1 -D 
parametric analyses were conducted for the entire turbine and exhaust nozzle flowpath of a typical 
advanced subsonic engine to understand the effects of important flow and chemistry variations. An 
example of the results from these efforts is shown in Figure 2. These calculations were initiated at the 
combustor exit plane and continued to the exit plane of the exhaust nozzle. The same engine operating 
data and chemistry is used as in the time scale analysis of Figure 1. 

It is evident that many trace species undergo significant change within the engine. Total 
(S0 3 +H 2 S0 4 )/S0 x at the engine exit for a baseline case was found to be 8.3% for a representative, modem 
subsonic engine. Comparatively, a total of 1 .2% NO y oxidation to HN0 2 and HN0 3 was computed. A 
variety of combustor exit speciation prescriptions were used to illustrate the competition among NO y , 

SO„, and CO x species for free radicals. From these studies, it was found that engine exit composition is 
very sensitive to the use of equilibrium versus nonequilibrium assumptions in establishing species initial 
conditions. Equilibrium specifications resulted in a reduced level of chemical activity through the turbine 


4 





and exhaust nozzle, a smaller amount of oxidized S0 2 as SO 3 and H 2 S0 4 (5.3% versus 8.3% using non- 
equilibrium NO x and CO), and a smaller amount of NO x as FTN0 2 and HNO 3 (approximately 0.2% versus 
1 . 2 %). 


In response to this sensitivity, a methodology for specifying initial conditions at the combustor exit was 
developed. This procedure utilizes empirical information in combination with combustor fuel type, 
operating point, thermodynamic specifications, and 1-D kinetics to estimate trace species concentration at 
the combustor exit. Recent combustor measurements of various trace species supported by the 
NASA/DERA collaboration will be used to assess the 


validity of this methodology. In our effort to improve 
these models, we have also updated reactions for SO x 
chemistry included in our analysis based on the work of 
the Dryer group at Princeton University. In comparison 
to the chemistry used for the 1-D results shown in Figure 
2 as well as the 2-D results in Figure 3, which will be 
discussed below, SO x production and destruction 
mechanisms are actually more active. In the balance, 
there is less oxidation than indicated by the values 
presented above. For the same conditions, total 
(S 03 +H 2 S 04 )/S 0 X at the engine exit is instead 4.5% 
rather than 8.3% for the baseline case, and 3.6% rather 
than 5.3% for the equilibrium initial condition case. This 
underscores the need for continued fundamental kinetic 
investigations and input as a complement to the 
proposed tasks. 

In the 1-D calculations, flow variables were also varied 
to show the importance of effects such as combustor exit 
temperature and cooling flow mass addition on the exit 
plane result. Such effects were relatively less important 
for these 1-D, flow-averaged studies than were changes 
in the method by which species initial conditions were 




1307 K 1620 K 0 37e-06 0.44e-0 


Cooled Blade Surfaces 


established. This is not the case if we consider the details Figure 3. 2-D Steady Analysis of Trace 

r . , + r-, , . .. r, Chemistry through a Turbine Blade Row. 

of the combustor exit profile and turbine fluid J 
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mechanics. Time scale analyses make it clear that the active conversion of sulfur species to S(VI) (i.e. 

SO3 and H2SO4) occurs primarily in the first stages of the high pressure turbine and can be perturbed by 
large temperature gradients in the flow in a nonlinear fashion. In general, the identified regions of active 
sulfur chemistry can be explored in more detail using the multi-dimensional code CNEWT. 

Representative 2-D, single turbine blade row simulations were performed to determine the potential 
impact of flow non-uniformities that cannot be. captured directly or modeled simply through 1-D analyses. 
Temperature non-uniformities that result from the use of an internal blade cooling strategy were 
investigated and revealed a significant impact in SO s chemistry. A comparison of results for cooled and 
uncooled blades is shown in Figure 3, where it can be noted that oxidation of SO 2 to SO 3 is enhanced in 
the cool regions of the blade boundary layer and wake. 

Comparisons of 1-D approximations to the 2-D turbine solutions were then carried out to help determine 
the extent to which current 1-D modeling capabilities can resolve changes in chemical composition. An 
example of such a 1-D reduction for the 2-D results of Figure 3 is shown in Figure 4. The results call into 
question the validity of 1-D averaged flow analysis for the highly non-uniform, unsteady flow fields of 
the turbine and exhaust nozzle. Because of the need for a high density computational grid, we have 
partially parallelized the CNEWT code for faster run times and have refined the calculation methodology 
in order to accommodate larger domains. 



Figure 4. Comparison of 1-D and Mass-Averaged 2-D Results for Similar Calculations. 

Further analyses of the impact of non-uniformities were undertaken to understand the persistence of the 
blade cooling impact through downstream blade rows and to determine the magnitude of that impact in 
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comparison to the much more significant circumferential and radial variations at the combustor exit. An 
unsteady, full-stage row calculation was performed for the nozzle guide vanes and rotor (together forming 
the first stator-rotor pair of blade rows) for a representative engine, but older and cooler than the engine 
example used in the discussion to this point. We used a different, older engine in support of the 
NASA/DERA measurement effort and expect, from our analyses as described above, that trace chemistry 
would be less active. 

From a technical standpoint, the goal of this full-stage analysis was to capture the unsteady interaction of 
the cooled wakes leaving the NG V as they flow over the rotor blades. These results were compared to the 
corresponding steady calculation where the NGV wakes are averaged before passing on to the rotor. In 
the former case, we allow the mixing processes in the flow to take place without intervention. In the latter 
case, mixing is enforced. The impact of mixing on non-uniform cooled flow sulfur conversion can thus be 
investigated for this specific engine. Figure 5 shows a snapshot of the unsteady, full-stage calculation. 

We find that mixing does not have an impact on the 
amount of sulfur conversion for this engine operating 
point, a total of 0.6% at the rotor exit. However, the 
fact that the steady and unsteady calculations agree so 
closely for this engine may be due to the low overall 
S(VI) conversion for this engine cycle. The results 
may also stem from the competing effects of larger 
regions of enhanced S0 3 for the unsteady case being 
countered by mixing and dilution. For other engine 
cycles, these unsteady effects may impact the SOj 
formation and this is an active area of investigation 
for our group. 

Another active area of investigation is the impact of 
combustor exit non-uniformities associated with the 
radial and circumferential temperature variability of 
the burner exit. Performing a similar comparison of 
cases as for the wake propagation investigation, this 
time for a higher temperature to increase conversion 
potential, our analysis indicates that mixing times are 

Figure 5. 2-D Unsteady Analysis of Trace 
Chemistry through a Turbine Stage. 
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long in comparison to the chemical reaction. As a result, persistent stratification of the flow originating at 
the combustor exit appears to inhibit S(VI) conversion. Mixing-out the flow between blade rows results in 
a conversion of 0.25% compared with a 0.16% conversion if the non-uniformity is allowed to propagate 
through the rotor. Also note that a higher conversion is realized for this higher temperature case for the 
same geometry and less stage flow-through residence time. 

In addition to gaseous trace chemistry, research efforts have also addressed the role of heterogeneous 
processes. One-dimensional flow models and unity probability heterogeneous rate parameters were used 
to estimate the maximum effect of heterogeneous reactions on trace species evolution in aircraft gas 
turbines. The analysis included reactions on soot particulates and turbine/nozzle material surfaces. 

Results for a representative advanced subsonic engine indicated that the net change in reactant mixing 
ratios due to heterogeneous reactions is < 10' 6 for 0 2 , C0 2 , and H 2 0, and < 10 10 for minor combustion 
products such as S0 2 and N0 2 . The change in the mixing ratios relative to the initial values were < 

0.01%. Since these estimates were based on heterogeneous reaction probabilities of unity, the actual 
changes will be even lower. Thus, we concluded that heterogeneous chemistry within the engine cannot 
explain the high conversion of S0 2 to SO ? which some wake models require to explain the observed 
levels of volatile aerosols. Furthermore, turbine heterogeneous processes will not affect exhaust NO x or 
NO y levels. 

From these results, we have proposed a modeling strategy for post-combustion trace species chemistry 
through the turbine and exhaust nozzle as defined in Figure 6. Like other gaseous chemistry, ion 
recombination is influenced by processes occurring through the turbine and exhaust nozzle. A one- 
dimensional model for aircraft turbine flows was used to estimate chemi-ion emission number densities at 
the engine exit plane. Results for representative subsonic and supersonic aircraft indicate that significant 
reductions in ion number density can occur through the turbine and exhaust nozzle. These reductions, 
typically a factor of 10-1000, result in chemi-ion emission indices at the engine exit that range between 
10+16 ions kg-fuel-1 and 1 0+18 ions kg-fuel-1, corresponding to a typical ion number density of 10+8 
ion cm-3. Available ground based measurements are in line with these estimates. Aircraft wake modeling 
suggests that chemi-ion emissions can alter aerosol microphysical processes in such a manner as to 
resolve discrepancies between model predictions and in-flight observations of volatile aerosol emissions. 
However, we predict peak ion number densities at the engine exit that are an order of magnitude below 
the levels required at the onset of ion-induced nucleation in the exhaust wake to completely resolve 
particle size/growth discrepancies. 
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A. Combustor Exit 


1. Steady, but radially/circumferentially variable properties 

2. Temperature defined by CFD, data, or parameterization 
with pressure based on operating point 

3. Species defined by non-kinetic initial condition methodology 
based on ... 

• combustor equivalence ratio and species El's where available 

• equilibrium model for other species (O, HOx, S02/S03, HNOx/NOx 

B. Characteristic time scale analysis to evaluate necessity for detail 

C. High-Pressure Turbine 

1. Model details reacting flow through turbine 

• 2-D CFD using finite-rate kinetics 

• 1-D calculations for freestream, 2-D+ for near-wall, wake 

2. Cycle-based parameter specification anchored by data where available 

D. Inter-turbine duct, low-pressure turbine, and exhaust nozzle 

1. 1-D kinetic calculations using averaged properties 


Figure 6. Proposed Modeling Strategy* 
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Abstract 

A consensus among scientists has formed that human activities have an influence on the atmosphere and are 
contributing to global warming. Recently, aviation has garnered a significant amount of attention since it is the 
fastest growing sector of transportation and it is thought to be an important component of such impacts. Trace 
species emitted by aircraft directly into the upper atmosphere are believed to have an enhanced influence 
relative to other common emissions generated by human activity near the surface of the earth. Prior work has 
concluded that trace species can change significantly in the engine post-combustor gas path. Research efforts to 
understand the role of aircraft gas turbine engine emissions in global climate change and the mechanisms by 
which their emissions are produced have been intensified. 

A modeling methodology was developed which can be used to predict the evolution of trace species in the post- 
combustor flow path of gas turbine engines. The modeling methodology incorporates various levels of detail to 
accurately and efficiently predict levels of intra-engine trace species by considering the key parameters 
affecting their evolution, specifically temperature, pressure, residence time, and species concentration. The core 
of the model is a computer code which performs a numerical solution to the fluid mechanics and chemical 
kinetics problems. The model is intended to improve the overall understanding of the fundamental physical 
processes that affect trace species evolution and to serve as a predictive design tool which can direct the 
development of new engine technologies which reduce undesirable aviation emissions. 

The development effort consisted of several refinements aimed at increasing the accuracy and usability of the 
model. The improvements included an improved species initial condition specification procedure, improved 
grid generation, capability for multiple inlet/exit planes, addition of parallel chemistry subroutines, revision of 
chemistry convergence criteria, addition of circumferentially-varying unsteady inlet conditions, and 
incorporation of a new chemical mechanism. 

Several validation exercises were performed which benchmark the capabilities of the model and test the added 
features. Past simulations of a duct and rotor were used as a baseline to verify several of the improvements. 
These simulations were used to improve the chemistry convergence criteria, perform a comparative study of 
several chemical mechanisms, and verify the periodic boundary conditions. The Princeton University Variable 
Pressure Flow Reactor was also simulated as a validation exercise in an attempt to evaluate the overall 
effectiveness of the model in predicting the flow field, mixing, and species concentrations in a well known and 
controlled environment. 

The model was used to perform the first complete post-combustor engine simulation in support of an engine 
measurement campaign. The results provided guidance for test parameters and measurement strategies. This 
simulation was further used to assess the effects of flow non-uniformity on the evolution of trace species in a 
typical aircraft engine and to refine the current modeling practice. Engine operating conditions, multi- 
dimensional non-uniformities, and the unsteady interaction of non-uniformities with subsequent blade rows 
were all found to influence trace species evolution. 

Thesis Supervisor: Professor Ian A. Waitz 

Title: Associate Professor of Aeronautics and Astronautics 
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Nomenclature 
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ROMAN: 


c 
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c p 

Specific heat at constant pressure 

D 
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ds 

Differential length along curve 

L 

Length 

M 
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m 

Mass flow rate 
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1 Introduction 


Aviation is currently the fastest growing transportation sector. Trace species emitted by aircraft directly into the 
upper atmosphere are believed to have an enhanced influence relative to other common emissions generated by 
human activity near the surface of the earth. Therefore, research efforts to understand the role of aircraft gas 
turbine engine emissions in global climate change and the mechanisms by which these emissions are produced 
have been amplified. Of particular interest are trace chemical constituents, which are understood to be 
particularly influential in affecting local air quality and the global atmosphere. The Atmospheric Effects of 
Aviation Project (AEAP), created by NASA to study the influence of aviation on the global atmosphere, 
initiated the effort to develop a model to predict the evolution of trace species in the post-combustor flow path 
of gas turbine engines. The work contained in this thesis is part of an ongoing research effort aimed at gaining 
accurate and informative data on these trace chemical species. Ultimately, a modeling methodology was 
developed and applied to a research engine in support of an engine test program being conducted for the project. 
A computational model which can predict the evolution of trace species within gas turbine engines is an 
essential tool which can be used to help understand intra-engine trace species evolution. The model can be used 
to investigate how the operating parameters and design characteristics of gas turbine engines affect their 
emissions. It can also help provide a physical understanding of trace species evolution and could be used to 
direct the development of new engine technologies that reduce aviation emissions. 

1 . 1 Chapter Overview 

The purpose of this chapter is to give a short introduction to aircraft emissions and their potential impacts on the 
local and global atmosphere. A brief history of the research on atmospheric effects of aviation is described. 
The composition of aircraft exhaust emissions is explained, the origins of a few trace species are presented, and 
a description of the exhaust deposition process is given. A few climatic and ozone effects are highlighted, with 
further details for subsonic and supersonic aircraft given. An overview of prior work to characterize trace 
species emissions is presented, including work done in near field plume and wake modeling and turbine and 
exhaust nozzle modeling. The importance of trace species modeling in the turbine and exhaust nozzle is 
established by several previous studies, which gives motivation for the current research. This chapter also lists 
the overall objectives and contributions of the current research effort which serves as a thesis roadmap. 
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Section 1.2 gives an overview of aircraft exhaust emissions, where these emissions are deposited, and the 
effects of both subsonic and supersonic aircraft exhaust emissions on the environment. Section 1.3 reviews 
some prior work on emissions characterization, specifically near field plume and wake modeling and turbine 
and exhaust nozzle modeling. Section 1.4 discusses the motivation for the work and lists the contributions of 
this thesis. Section 1.5 is the chapter summary. 


1.2 Atmospheric Effects Of Aviation 

Aircraft emissions can impact the environment both locally and globally. The aviation industry is becoming a 
more significant part the world economy, in terms of both commercial and military activity, and it contributes to 
these atmospheric effects primarily via the emissions released from gas turbine engines. Future growth in 
aircraft fleets is expected to increase these effects [3]. Therefore, a complete understanding of aircraft 
emissions and effects of these emissions on the atmosphere is necessary in assessing the current and possible 
future effects of aircraft. The need to characterize trace species evolution has been heightened because trace 
species concentrations are known to change through the post-combustor gas path of the engine (turbine and 
exhaust nozzle), exhaust plume, and aircraft wake prior to deposition in the atmosphere. These processes 
collectively govern the impact of engine exhaust which can affect the global climate. Although there is still 
uncertainty surrounding the effects of emissions impacts on the global climate, a number of scientific and 
technological investigations have helped clarify the issue. The following sections discuss the reasons for 
focusing on aviation emissions, provide a brief overview of potentially influential aircraft emissions, describe 
atmospheric deposition, and highlight impacts of both subsonic and supersonic aircraft on the atmosphere. 

The history of research on this issue dates back to the 1972, when the US Department of Transportation 
commissioned the Climatic Impact Assessment Program (CIAP) to record all the scientific research on the 
effects of aviation on the atmosphere up to that date. After CIAP was completed in 1975, NASA became the 
primary agent for atmospheric research and its Upper Atmosphere Research Program has been active ever since. 
NASA's research has also been incorporated into reports of the Intergovernmental Panel on Climate Change 
(IPCC) [2]. In 1988, the AEAP was established to comprehensively predict the atmospheric impacts of future 
aircraft. The program is a collaborative effort between several government agencies, academic institutions, and 
industry which was created to study the influence of aviation on the global atmosphere. 

The IPCC report [2], entitled “Aviation And The Global Atmosphere," is the most comprehensive source of 
detailed background information related to this thesis topic. This chapter highlights a few of the relevant topics, 
refer to the original report for more information. 
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Operating Condition 

Species 

Idle 

Take-Off 

Cruise 

co 2 

3160 

3160 

3160 

h 2 o 

1230 

1230 

1230 

CO 

25(10-65) 

< 1 

1-3.5 

HC (as CH 4 ) 

4 (0-12) 

<0.5 

0.2-1. 3 

NO x (as N0 2 ) 

3-6 

10-65 

7.9-15.4 

SOx (as S0 2 ) 

1 

1 

1 


Table 1.1 Typical emission index levels (g/kg fuel) [3] 


The exact concentrations of the species emitted from the engine are difficult to determine since each combustor 
behaves differently, and an efficient scheme for collecting accurate data is yet to be established. The levels of 
the trace species NO y and SO x emitted at the engine exit plane are kinetically determined by the amount of time 
spent in each section of the flow path and the time variability of the temperature and pressure histories realized 
through the engine (from combustor to the exhaust nozzle). Therefore, the engine specifications and 
thermodynamic cycle have a direct effect on the emissions. For example, higher temperatures result in greater 
radical concentrations at the engine exit [1], 

1.2.2 Altitude Deposition Process 

The area of the atmosphere in which emissions are deposited determines the effects they have upon the 
atmosphere. Aircraft emissions are released at altitudes ranging from the lower troposphere into the middle 
stratosphere. The troposphere is the region of the atmosphere located between the ground and about 10 to 15 
km (thinner in polar regions and thicker of the equator), the stratosphere is located between the altitudes of 10 to 
15 km to about 50 km, and the tropopause is the area dividing these regions. Typical subsonic aircraft operate 
primarily in the upper troposphere and the lower stratosphere, while supersonic aircraft, such as the Concorde or 
military aircraft, fly in the middle stratosphere, the exact height being a function of the design Mach number. 

When exhaust first exits an aircraft engine, it enters a region termed the near-field plume region, where the 
engine plume and the aircraft wake do not interact. The initial, turbulent shear induced dispersion of the 
exhaust jet is more important than any aerodynamic perturbation from the aircraft. Eventually, as the 
momentum of the exhaust plume is diffused, the vorticity shed from the aircraft wings begins to influence the 
jet plume. In this wake vortex region, the plumes are entrained and confined in the two counter-rotating tip 
vortices formed by the self-induced vortex sheet roll-up of the aircraft wake. This structure, at cruise altitude, 
sinks about 100 m for subsonic aircraft and several hundred meters for supersonic aircraft through a self- 
induced downwash. After a certain amount of time, the structure begins to deteriorate and break up due to 
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Radiative forcing is a measure of the importance of a potential climate change mechanism. It expresses the 
perturbation or change to the energy balance of the earth-atmosphere system in watts per square meter (W/m 2 ). 
Positive values of radiative forcing imply a net warming, while negative values imply cooling. Figure 1.2 
shows the radiative forcing from different aircraft emissions in 1992 and a projection for 2050. 


o.io 

0.08 

r* 

E 0.06 

£ 

M 0.04 

| 0.02 

u. 

► 0 00 
*<3 

a 

fiS 

- 0.04 

- 0.06 


1992 


& 

C02 


3.5 %of TOTAL 
FORCING DUE 


. 


TO MAN 


\ 


ch 4 


T 




Dtrtct 

Sulfate T 

. _ 

KSO Contrails Cirrus X Direct Total 

Clouds Soot (without 

cimjs 
clouds) 


good fair poor poor fair very fair fair 
poor 



Figure 1.2 Radiative forcing of various trace species [2] 

The impacts of C0 2 are well known due to its long atmospheric residence time (-100 years). It becomes well 
mixed throughout the atmosphere and thus the effects of its emission from aircraft are impossible to 
differentiate from the same quantity of carbon dioxide emitted by any other source [3]. Emissions such as NO x , 
SO x , and water vapor have shorter residence times and remain concentrated near the aircraft. Thus they lead to 
changes in radiative forcing localized near the flight route as opposed to emissions that are globally mixed. 

Increases in particles emitted from aircraft have mixed effects: soot tends to warm the earth’s surface, while 
sulfate particles tend to cool it. While direct effects of these particles are believed to be small, increases in their 
emissions by jet aircraft may potentially influence the formation of clouds and contrails, thin white-line clouds 
seen behind jet aircraft in the upper atmosphere, which may contribute to future climate change [4]. 

It is difficult to discern the impact of aircraft emissions upon changes in ozone concentration since a decrease in 
the ozone in the upper and middle stratosphere lead to global warming whereas decreases in ozone in the upper 
troposphere/lower stratosphere lead to global cooling. In addition, supersonic aircraft, flying in the middle 
stratosphere, emit species that lead to the destruction of ozone whereas subsonic aircraft, flying in the upper 
troposphere/lower stratosphere, emit species that lead to the formation of ozone. This formation and destruction 
of ozone along with emission of species such as C0 2 and H 2 0 impact the radiative balance of the earth- 
atmosphere system in the direction of global warming [1]. 
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Global warming due to ozone and water vapor in the upper troposphere/lower stratosphere altitudes is greater 
than at the surface due to low background concentrations, longer residence times, and large radiative efficiency 
near the tropopause as compared to the surface. 

1.23.2 Supersonic Aircraft Effects 

Supersonic aircraft fly in the middle stratosphere, near the peak concentration of atmospheric ozone, the exact 
height being a function of the design Mach number. Figure 1.3 shows the variation of ozone concentration with 
altitude. The peak ozone range of Mach 2.4 to 4 is where the most recently proposed US high speed-speed civil 
transport (HSCT) designs operate. 



Figure 1.3 Variation of ozone concentration with altitude [7] 

NOx is the trace species that is the most significant contributor to the reduction of ozone, by the following 
catalytic cycle [7]: 

NO + 0 3 — » N0 2 + 0 2 
N0 2 + hv — » NO +0 2 
and 

N0 2 + 0 — > NO + 0 2 

In this cycle, no NO is destroyed, but 0 3 is combined with O to form 0 2 . This problem was first publicized in 
1971, when Johnston (1971) calculated that emissions of oxides of nitrogen from a fleet of five hundred 
supersonic transport aircraft flying at an altitude of 20 kilometers could reduce stratospheric ozone by a global 
average of 10 to 20 percent [8]. 
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The chemistry that occurs during the deposition process in the plume and wake regions occurs at higher 
temperatures (500 to 600 K) than atmospheric (200 to 300 K), but lower than the temperatures encountered 
during the combustion process. At these temperatures, the inclusion of heterogeneous chemistry and dynamics 
for aerosol formation within the overall processing mechanism is necessary. Primary pollutant and trace species 
are highly active in the atmosphere, involving both gas-phase and heterogeneous chemistry. The hydroxyl 
(OH) radical, emitted either with the exhaust stream or formed by photochemical breakdown of atmospheric or 
exhaust constituents (0 3 , H 2 0) initiates and propagates several reactions including gas phase conversions of 
NO, to NO y which occur in general though the reactions: 

NO + OH + M <— » HN0 2 + M 
N0 2 + OH + M <— > HNO, + M 
N0 2 + O <r -> N0 3 :N0 2 + N0 3 + M 4- -> N 2 0 5 
Water vapor also has an important role. The primary S0 2 gas phase oxidation reactions are. 

S0 2 + O + M <r -> S0 3 + M 

S0 2 + OH +M 4~» HS0 3 + M:HS0 3 + O, S0 2 + S0 3 

so 3 +h 2 o <—> h 2 so 4 

Competition for the O and OH radical also occurs with the conversion of CO to C0 2 via the reactions: 

CO + OH 4~» C0 2 +H 
CO + O + M<~ » C0 2 + M 

The gas phase chemistry consists of other conversion processes, most notably the inevitable conversion of NO 
to N0 2 . Oxidation with the HO, family can occur in many different ways including reactions with O and H. 

Oxidation of S0 2 to SO, and H 2 S0 4 is of greater interest to modeling and exhaust sampling now since studies 
have shown that increases in S0 3 concentration lead to higher total aerosol surface area, and greater activation 
of soot particles. The enhancement of aerosols, as mentioned before, can have important consequences for 
ozone and climatic impacts. Gas phase species such as N 2 O s and N0 3 and S0 3 are readily converted in the 
presence of water containing aerosols to form aqueous HN0 3 and H 2 S0 4 . 

1.3.2 Near Field Plume And Wake Modeling 

Most in-situ and stationary measurements have been done for subsonic aircraft. These efforts concentrated on 
chemical activity occurring outside the engine, while only a limited understanding could be gained about the 
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The research by Lukachko [i ), entitled "Research on the Science And Politics of the Atmospheric Effects of 
Aviation Debate," is a primary source of work on turbine and exhaust nozzle modeling. The research 
performed in this thesis is an extension of the work by Lukachko. Lukachko also presents a more detailed 
summary of other prior work on intra-engine modeling. Refer to the original document for more information 
related to this section. 


1.4 Motivation 


As discussed in Section 1.2, important chemical reactions occur in the post combustor flow path which 
potentially alter the trace species concentration at the nozzle exit plane. Trace species play an important role in 
the problem of emissions characterization because SO* and NO y families have acids and radicals that influence 
aerosol formation and oxidation processes, which effect the processes that occur outside the aircraft and on a 
global scale. 

The two principal components in the efforts to characterize aircraft emission are physical tests and analytical 
modeling. Combustor emissions data is gathered using uninstalled combustor -only and full-scale engine tests. 
Full-scale engine tests at altitude conditions are the most relevant, so some tests have been performed in-flight 
(in situ). However, in-flight experiments are limited in accuracy due to the distances that must remain between 
the two aircraft. From combustor-only tests and full-scale engine tests, useful data can be gathered to determine 
the products that travel though the turbine and exhaust nozzle. There has been few stationary tests preformed to 
measure NO*, CO and HC [10J, [11]. Corrections such as those for humidity and other ambient conditions must 
be applied to make the test results comparable to each other. 

All the above approaches are valid, but there are even more complexities within a combustor. A diversity of 
both reactants and products of combustion are involved in hundreds of simultaneous reactions. So while the 
aforementioned species are the majority of the resulting compounds, there are many minor constituents that 
exist in very small, trace amounts. These species form because of the non-equilibrium chemical state of the 
exhaust gas and include the hydroxy family (HO*), nitrogen compounds (NO y ), and the sulfur oxide family 
(SO*). These species are emitted at the engine exit plane in levels that are kinetically determined, governed by 
the time variability of temperature and pressure histories realized within the combustor and through the turbine 
and exhaust nozzle of the engine as well as the amount of time spent in any particular section of flow path. 
Therefore, the thermodynamic cycle on which a particular engine is based will have a direct effect on these 
emissions. Higher temperatures generally result in greater radical speciations and less oxidation progress at the 
engine exit. 


Due to lack of suitable test opportunities and adequate instrumentation only limited data has been acquired. 
Since experimental testing is expensive and time consuming, computer modeling is an efficient method for trace 
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1 .4.2.1 Modeling Methodology 

A modeling methodology was developed to simulate trace species evolution in the post-combustor flow path of 
a gas turbine engine. The methodology involves a fundamental time scale analysis which can help guide the 
overall modeling strategy and provide insight into the mechanisms influencing trace species evolution. The 
methodology includes specifying chemical initial conditions, fluid boundary conditions, and performing 
numerical solutions to the fluid mechanics and chemical kinetics problems at various levels of detail using 
either high fidelity, low fidelity, or a combination of modeling techniques. Key parameters affecting levels of 
intra-engine trace species were identified and incorporated into the model thereby improving the accuracy and 
increasing the efficiency of the model. 

1.4.2.2 Modeling Tools 

Usability and modeling accuracy improvements to the existing modeling tools were made, which included: 

• incorporating a 3D fully unstructured grid generation tool which allows simulations of complex 
arbitrary geometries 

• adding multiple inlet/exit capability which allows simulation of mixing flows 

• adding parallel chemistry routines to improve execution time 

• refining circumferentially varying unsteady inlet conditions to improve modeling accuracy and further 
investigate the effects of flow non-uniformity on chemistry 

• improving chemistry initial condition specification procedure 

• adopting improved chemistry convergence criteria 

• selecting a new chemical mechanism which incorporates new kinetic data for sulfur chemistry 
provided by Princeton University 

1.4.2.3 Validation 

Some model limitations were established through a validation exercise that attempted to model the Princeton 
University Variable Pressure Flow Reactor. The new grid generator and multiple inlet/exit modifications were 
demonstrated, but the fluid, mixing, and chemistry modeling could not be validated due to compressibility 
limits of the code, as well as, the possible unsteady flow regime of the reactor. 

Further validation efforts were performed on a representative ID engine cycle and 2D turbine geometry. They 
included evaluation of various chemical mechanisms, species convergence criteria, periodic boundary 
conditions, and the parallel code modifications. 

1.4.2.4 Modeling The NASA/DERA Engine Test 

The modeling methodology was applied to an engine to perform the first complete post-combustor engine 
simulation. Trace species concentrations were predicted at the nozzle exit for several test cases based on an 
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• The region of the atmosphere where emissions are deposited determines their effects on the atmosphere. A 
decrease in the ozone concentration in the upper and middle stratosphere leads to global warming whereas 
a decrease in ozone in the upper troposphere/lower stratosphere leads to global cooling. Overall, this 
formation and destruction of ozone impact the earth-atmosphere system in the direction of global warming. 

• Subsonic aircraft affect chemical and radiative processes mainly through emission of NO x , but also 
H : 0, SO x , and soot. These emissions lead to the formation of ozone in the upper troposphere/lower 
stratosphere. 

• Supersonic aircraft fly in the middle stratosphere near the peak concentration of atmospheric ozone. 
Through NO x emissions they lead to the destruction of ozone 

• Prior work has involved modeling and measuring trace species in the near field plume and wake of an 
aircraft and computational fluid-chemical modeling in the intra-engine environment. Results suggest that 
most trace species chemistry occurs early, thus there is a need for better characterization of trace species in 
the intra-engine environment. 

• See the thesis by Lukachko [1], “Research on the Science And Politics of the Atmospheric Effects of 
Aviation Debate,” for more background and preliminary work on intra-engine modeling. 

• Computational modeling is an efficient means to acquire detailed information about intra-engine trace 
species chemistry since experimental techniques are complicated, time consuming, and expensive. 

• The contributions of this thesis include: 

• Continued development the post-combustor trace species modeling methodology 

• Improvements to modeling tools 

• Continued validation of the model and improvements to the modeling tools 

• Performed the first complete post-combustor engine simulation to support an engine test campaign 

• Investigation of a few physical phenomena which influence trace species evolution 

They are aimed at further characterization of aircraft emissions, in particular modeling the evolution of 
trace species in the post-combustor flow path of gas turbine engines. 
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2 Modeling Methodology 


The outcomes of this work are several tools which can be used to evaluate and improve the design of gas 
turbine engines with respect to trace species emissions as well as investigate the effects of various fluid or 
chemical parameters on the evolution of trace species. Section 1.4 detailed the motivation and advantages of 
modeling the evolution of species in the post-combustor flow path. In many cases, the numerical tools can be 
used as a fast and inexpensive means to evaluate potential designs or to probe specific research questions. The 
procedure developed is adaptable to several levels of modeling detail (or accuracy) and its capability can easily 
be extended beyond its current state. Many aspects of the models have been previously validated. The flow 
solver itself is a well-established turbomachinery code [51]. Validation exercises done by Lukachko [1] 
verified proper calculation of kinetics as specified in the chemical mechanism, species convection, species 
diffusion, mixing, influence of heat transfer or boundary layers, and some numerical issues on the chemistry 
solution. The validation efforts are ongoing, with several code improvements verified through test cases in this 
thesis (see Section 4) and the ensuing effort to benchmark the code against actual engine test data in the 
NASA/DERA engine test simulation. The modeling methodology is a collaboration of work from many 
researchers [3], in particular it is a continuation of work by Lukachko [1]. 

2.1 Chapter Overview 

The purpose of this chapter is to provide an overview of the intra-engine species evolution modeling 
methodology presented in this thesis. Various physical processes that can be incorporated into the models, as 
well as their potential effects, are listed. The basic procedure for a time scale analysis is defined. The relevance 
and advantages of performing a preliminary time scale analysis are presented, using several examples to 
demonstrate the utility of the time scale analysis. This chapter also discusses some of the details associated 
with the specification of the initial chemical conditions at the combustor exit plane. Finally, a broad overview 
of both the low and high fidelity modeling techniques is presented. 
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Modeling Element Or Physical Process 

Potential Effects On Species Evolution 

Chemical mechanism 

Important to capture appropriate reactions for accuracy 
[Section 4.3] 

Boundary layers and wakes 

Chanaes residence time and thermodynamic state 

Engine operating condition 

Higher combustor exit temperatures enhance SO 3 
conversion [3], [Section 5.6.1] 

Combustor temperature non-uniformities 

Induces iarqe scale non-uniformity [Section 5.6.2] 

Downstream persistence of non-uniformities 

Various implications [Section 5.6.3] 

Blade cooling (thermal) 

Enhances SO 3 conversion [1], [Section 5.6.2] 

Heat release or secondary reactions 

Can be important locally, details unknown [65] 

Chemiions 

Neqliqible effect [58] 

Heterogeneous chemistry 

Negligible effect [58] 

Blade cooling (mass addition) 

Unknown 

Turbulence 

Unknown 


Table 2.1 Some possible modeling elements influencing intra-engine chemistry 


Table 2.1 is not meant to be a comprehensive list of all physical processes and modeling elements which can 
impact intra-engine chemistry and their effects, but rather a list of a few items which were brought out in the 
course of this research. It is important to note that the effects of many of the modeling elements or physical 
processes are often inter related and may not be important in all situations. Hence, the term “potential effects” 
meaning that the effect was observed for a given set of conditions modeled, however, is not necessarily always 
the case nor the only possible effect. The modeling methodology and tools developed in this thesis incorporate 
many of the elements listed in Table 2. 1 or can be readily adapted to include them. The modeling exercises in 
Section 4, Validation, and Section 5, Modeling the NASA/DERA Engine Test, explore the impacts of several 
these modeling elements and physical processes. The exercises were used to investigate the magnitude of the 
impacts and develop appropriate models that can be used to simulate the relevant processes. 

The main effects typically derive from the interaction of fluid mechanics and chemistry. The fluid mechanical 
phenomena can alter the chemistry through parameters involved in the Law of Mass Action or Arrhenius 
Equation, specifically through the temperature, pressure, residence time, or chemical composition which will be 
discussed in more detail in Section 3. The scale (i.e. size) of the fluid mechanical phenomena can have an 
impact on the magnitude of the effects. 


2.3 Procedure 

The steps outlined in Sections 2.3.1 through 2.3.3 constitute the modeling methodology for studying the trace 
species evolution in the post-combustor flow path of gas turbine engines. Figure 2.1 is a graphical depiction of 
the procedure. 
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analysis is benificial. Furthermore, the time scale analysis can be used to assess the impact of detailed flow 
processes and influence of various chemistries. The time scale analysis includes three fundamental controls on 
chemical evolution; thermodynamic potential, chemical kinetics, and residence time. 

The concept and method of the time scale analysis were developed by Lukachko [31] and are further detailed in 
unpublished internal research documentation. An endeavor to formally publish the work is planned in the near 
future. 

2.3.1.1 Basic Procedure 

The thermodynamic potential is the difference between a given composition and the equilibrium composition. 
This potential change in a given species and the chemical time scale required for the change can be mapped out 
over a large temperature and pressure parameter space using only the key reactions of interest. Figure 2.2 
shows a typical plot of the thermodynamic potential, expressed as % SO^SO* at equilibrium, over a parameter 
space relevant to gas turbine engines. Several additional simplifying assumptions about the chemistry can be 
made in this preliminary analysis. In this case, SO 3 chemistry is being investigated using only the O and OH 
pathways from [28], using only forward reactions, and assuming equilibrium radical levels. 


— 0.001 ms 

— 0.01 ms 

— 0.1 ms 

— 1 ms 

— 10 ms 

— 1 00 ms 

10 20 30 40 50 

Pressure (atm) 

GREY = contours of Tso3 
BLUE = lines of constant % S 03 /S 0 *| eq 
RED = combustor exit conditions for various engines 

Figure 2.2 Typical thermodynamic potential and chemical time scale plot underlying a time scale analysis 

for S Oj (courtesy of Stephen Lukachko) 
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Figure 2.3 Severity parameter for S0 3 from typical time scale analysis at two engine stations (courtesy of 

Stephen Lukachko) 

A curve representing the post-combustor flow path can be plotted on the map. The example in Figure 2.3 
shows a curve (blue) for a representative engine at cruise. The severity parameter plot can be used to point out 
several aspects of the selected chemistry. First, notice that the contours of constant AS(VI) shift as one 
traverses through the engine. The shift is due to the amount of conversion which previously occurred, the new 
thermodynamic state, and to a lesser extent the change in residence time at the specific engine station. Second, 
the contours indicate a zone of maximum SOj formation. In this case, there is more than 10% conversion for 
temperatures of 1,000-1,300K and pressures over 20 atm. Finally, the severity parameter plot at the low 
pressure turbine entrance indicates that there is a negligible amount of further conversion possible (i.e. the 
active chemistry has completed in the high pressure turbine). This would suggest that the need for higher 
fidelity modeling should be concentrated in the high pressure turbine region. 

Finally, by integrating along the path shown in the figure preliminary estimates can be made for the evolution of 
the species. The simplifying assumptions inherent in the time scale analysis must be considered when making 
quantitative estimates. Thus, the time scale analyses are primarily used as a guide for further modeling. 
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Figure 2.5 Investing the effects of operating condition on chemistry using time scale analysis (courtesy of 

Stephen Lukachko) 

As a final example of the utility of the time scale analysis, the influences of different chemistries can be 
assessed. Figure 2.6 shows a comparison of two time scale analyses, one using only the S0 2 + O pathway for 
SO3 formation, and the second using both S0 2 + O and S0 2 + OH pathways. In this situation, which assumes 
an equilibrium level of OH/O of about 10 , it is evident that the OH route is more active than the O route. Thus, 
the time scale analysis can give insight as to which chemical pathways are most important for the species subset 
being studied. 
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Some minor improvements to the code were made to address the rapid initial shift in some species found in past 
simulations. A overview of the initial conditions code procedure is given in Figure 2.7. 



Figure 2.7 Overview of chemistry initial condition specification 

The code starts with the equilibrium composition at the combustor exit conditions. The composition is then 
perturbed using specified emissions indices or species ratios. The addition of the trace species is balanced by 
adjusting the levels of the primary combustion products. Finally, a constant temperature and pressure kinetics 
calculation is run for a very short time to reach a stable non-equilibrium composition. Table 2.2 lists the input 
parameters for the code. 

Initial Conditions Code Input Parameters 

Temperature (K) 

Pressure (atm) 

b (equivalence ratio) 

q (combustion effeciency) 

CO/Hz/HC ratio 

EI(NOx) (g/kg fuel) 

NQ/NQx ratio 

EI(S) (g/kg fuel) 

SO3/SO, ratio 

EI(CO) (g/kg fuel) 

Table 2.2 Input parameters to initial conditions code 

The improvements to the technique in [1] were made by changing some of the input parameters used to perturb 
the equilibrium composition and also in performing the kinetics calculation to stabilize the composition, see 
Section 2.3.2 for more details. The non-equilibrium solution initially output from the code is inherently 
unstable. Although some attempt is made to redistribute the mixture to a “stable non-equilibrium state" when 
applying the EEs, the entire set of reactions in the mechanism can not be considered easily. Running the 
kinetics calculation for a short time allows the species to make the initial rapid re-adjustment to a “stable non- 
equilibrium state" with a negligible change to the overall composition. The duration of the kinetics calculation 
was chosen by plotting each species versus time for the constant temperature and pressure calculation and then 
moving back to the point where each species was deemed to be outside of the “initial shift region." Figure 2.8 
and Figure 2.9 show the mole fraction of a few species and several of the ETs versus time during the stabilizing 
kinetics calculation. The left column of graphs show the result of the entire calculation, w hich was chosen as 5 
ms to exceed a typical blade passage flow-through time. The green lines indicate 0.30 ms, and the graphs on 
the right column are an enlarged detail of the 0 to 0.3 ms region. For this case, the red line at 0.02 ms was 


47 











Mole Fraction, HOSO Mole Fraction, S03 


INI 

o 

co 


o 

e 

E 

Li. 




o 

2 



Time (ms) 


Time (ms) 


Figure 2.8 Representative constant temperature and pressure kinetics calculations for selected species for 

initial conditions formulation 















2.3.3 High And Low Fidelity Modeling 


The modeling techniques are divided up into two categories, high fidelity and low fidelity. The word fidelity 
can mean accuracy in details or exactness, and high or low fidelity represents the degree to which a 
reproduction is faithful to the original. In the sense of the modeling methodology being described here, a low 
fidelity model is a simple ID simulation which can in general only track average quantities. A high fidelity 
model is any model which can directly simulate more detailed multi-dimensional flow features and thus account 
for chemical or flow non-uniformities. The high fidelity model, containing more details of the actual system 
being modeled is expected to be more accurate, although, the accuracy of the high fidelity model is also based 
on its ability to accurately represent the details being simulated. 

In contrast to a simplified time scale analysis which is based on chemical kinetics alone, the low and high 
fidelity models are based on integrating the equations of mass, momentum, and energy over a discrete domain 
(see Section 3.2) coupled with the chemical kinetics solution. For a low fidelity model the integration occurs 
over a ID profile and for a high fidelity model it occurs over a 2D or 3D grid. 

A time scale analysis can be used as a guide to map out the modeling strategy. It is desirable to concentrate 
modeling resources on areas of the engine which have a high degree of chemical activity for the species of 
interest. A 2D or 3D reacting fluid-chemical simulation can be done on the regions of the engine flow path 
targeted for higher fidelity models. A suite of numerical tools to be discussed in more detail in Section 3 was 
further developed in the context of this thesis. It consists of a grid generator, fluid and chemical kinetic solvers, 
and flow visualization codes which can be used to do such higher order modeling. The higher order modeling 
can include the effects of many physical phenomena typically encountered in a turbomachinery environment 
such as complex geometries, viscous walls, cooled surfaces, and unsteadiness (see Section 2.2), These 
phenomena can cause local temperature, pressure, or species gradients, which can influence the subsequent 
evolution of chemical species. Higher fidelity modeling requires more user effort to set up and post-process and 
is more computationally intensive. These simulations take roughly 10 to 100 times an investment in resources 
compared to a low fidelity model. A typical high fidelity simulation may take several days to weeks to setup 
and run whereas a low fidelity simulation can be completed in a matter of hours or days (see Section 3.3.3 for 
more details). However, the high fidelity models can increase the predictive capability of the modeling 
methodology and help provide insight into the mechanisms responsible for altering chemical evolution if 
applied in regions where chemistry is suspected to be more reactive. 

The complete details of the procedure needed to perform a higher fidelity simulation are numerous and highly 
dependent on the particular modeling situation. Many of these details are discussed further elsewhere in this 
thesis (especially Section 3) or in several external references such as [1], [48], and [49]. The basic process 
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• Time scale analyses are useful in formulating a overall modeling strategy by indicating the critical areas of 
the engine flow path to investigate with higher fidelity modeling. 

• A time scale analysis is a simple parametric study which incorporates thermodynamic potential, chemical 
kinetics, and residence time. The output is a plot of severity parameter for a given species of interest, 
although, a preliminary estimate of species evolution can be obtained. 

• Time scale analyses can also be used to provide insight into the physical phenomena influencing trace 
species evolution. A few example time scale analyses were used to show: 

• Blade cooling can enhance sulfur conversion. 

• Engine operating condition can influence sulfur conversion, take-off conditions can have higher sulfur 
conversion than cruise conditions. 

• The OH pathway to SO3 can be more active than the O pathway. 

• A simple method for specifying chemistry initial conditions involves perturbing an equilibrium 
composition , re-balancing the species based on mass, and doing a stabilizing kinetics calculation. 

• High and low fidelity models involve integrating the equations of mass, momentum, and energy over a 
discrete domain coupled with chemical kinetics to simulate intra-engine chemistry. A high fidelity model 
can incorporate many multi-dimensional effects (i.e. non-uniformities) while a low fidelity model uses 
average quantities. 
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3 Modeling Tools 


Underlying the modeling methodology presented in Section 2 are several computer codes which are used to 
generate grids, specify initial conditions, numerically solve the fluid mechanics and chemical kinetics problems, 
and visualize the results for geometries of interest to intra-engine chemistry. 

3.1 Chapter Overview 

The purpose of this chapter is to present a brief overview of the modeling tools and discuss the improvements 
made to these tools during the course of this research. The existing modeling tools consist of the PRE grid 
generation code, PROCESS grid pre-conditioner code, POST grid post-processor code, CNEWT and 
CALCHEM flow-chemistry solvers, and CNEWTVS6 flow visualization code. The improvements made to the 
existing tools allow for more complicated geometries, multiple inlet and exits, improved execution time for 
solutions with chemistry, and a specification of a pressure, temperature, and species concentration deficits 
associated with a wake of an upstream blade or blade row. 

In summary, the following are the main points of the chapter: 

• CNEWT and CALCHEM are useful modeling tools for studying intra-engine trace species flow-chemistry 
problems. Other supporting codes necessary for turbine chemistry modeling are PRE, PROCESS, POST, 
and CNEWTVS6. 

• For more details on the structure and numerical mechanics of the codes, refer to “Research on the Science 
And Politics of the Atmospheric Effects of Aviation Debate" by Lukachko [1] as the primary source of 
information. 

• In order to provide further insight into trace species characterization, improvements were made to existing 
modeling tools which include: 

• Obtaining a new grid generator 

• Adding multiple inlet/exit capability 

• Parallelizing chemistry subroutines 

• Implementing a wake model boundary condition 
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3.3 Improvements 


Several improvements to the modeling tools were made during the course of this research aimed at extending 
the capability and improving the accuracy of the modeling tools. A few extensions are discussed elsewhere in 
this thesis; first, the procedure for specifying chemistry initial conditions was improved and as detailed in 
Section 2.3.2, next, a new chemical mechanism for turbine chemistry modeling was adopted as described in 
Section 4.3, finally, the chemistry convergence criteria were refined as discussed in Section 4.4. This section 
describes several additional improvements to the modeling tools. In particular, a new grid generator was 
obtained, a option for multiple iniets/exits was added, the chemistry subroutines were parallelized, and a wake 
model boundary condition was implemented into the codes. 

3.3.1 CAD Based Fully Unstructured Grid Generation 

Grid construction is important since it determines the resolution (i.e. accuracy) of the flow solution and impacts 
the computational time, therefore an accurate and efficient grid construction is desired. The grid generator used 
for previous investigations, called PRE, is best suited to turbomachinery blade passages. The grids typically 
had four sides, with a blade form either entirely within the domain (block 5 input format) or on the side 
boundaries (block 3 input format). The grids are initiated as a structured mesh composed of hexahedral cells 
and are transformed into tetrahedral grids simply by dividing the hexahedrons. Hence they were “pseudo- 
unstructured,” since although being tetrahedral, the cells were still in an orderly arrangement based on the initial 
hexahedral grid. By piecing together several blocks, a more complicated grid can be obtained, such as a blade 
with tip clearance. 

The structured or “pseudo-unstructured” grid is less flexible in the types of geometry it can mesh and it may 
cause some numerical issues. In previous work, of validation runs using a PRE meshed geometries, the 
temperature and velocity exit plane contour plots showed an undesirable result of hotter, slower flow in the 
corners perpendicular to the mesh diagonal and grid scale oscillations in solution variables which could be 
linked to the structure of the grid [1]. 

New grid generation tools developed at Cambridge University were implemented to overcome these potential 
problems and extend the codes to more complex arbitrary geometries [35] and [66]. This set of grid tools 
requires an initial CAD-generated geometry, which can be created with typical commercial solid modeling 
software. Thus, virtually any arbitrary geometry can be created and meshed. There are essentially three main 
aspects that comprise the new grid generation tools; the solid model, the surface mesh, and the volume mesh. 
Figure 3.1 shows the difference between a pseudo-unstructured grid made with PRE and fully unstructured grid 
made with the new grid generator. 
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Volume Mesh 


The volume meshing routine uses the binary surface mesh file to create the final 3D grid. There are no grid 
controls which can be applied during volume meshing, except for specification of viscous layers. Rather, the 
volume grid generator attempts to make an isotropic mesh with cells that are a uniform size. There are three 
stages involved in the volume meshing which perform the following basic steps: 


Stage 1: 

• creates 8 construction nodes encompassing the surface mesh 

• starts with the initial construction cell and connects surface nodes into tetrahedrons using 3D Delauney 
triangularization by marching across domain, grabbing, and constraining each surface node 

Stage 2: 

• generates viscous layers by extruding nodes near walls flagged as viscous 

• moving front 
Stage 3: 

• in-fills domain by inserting new nodes into cells that fail quality measures 

• locates tetrahedral cell that the new node is in, then divides into many new tetrahedral cells 

• edge swapping based on 3D Delauney to improve local mesh quality 


Figure 3.2 shows a mesh after a stage 1 dump and two meshes, one without viscous layers and one with viscous 
layers, after a stage 3 dump. 



Figure 3.2 Volume mesh stage dumps 


Note that the cells in the cross-section may not look uniform, because the cutting plane intersects some cells 
near a vertex (giving the appearance of a very small cell), some cells near an edge (giving the appearance of a 
sliver shaped cell), as well as some near the center of the cell. The volume mesh routine outputs a mesh 

connectivity variable file (.mcv) which contains the 3D mesh in NEWT format. More details of the new grid 

» 

generation tools can be found in [66], 

During the period of work on this thesis, the new grid generation tools were in a developmental stage. They 
worked well for simpler geometries, however, some difficulties were encountered for more complex 
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case demonstrated that the multiple inlet/exit code additions to CNEWT work properly, however, further 


validation should be pursued. 



Figure 3.3 Multiple inlet/exit test case, X direction velocity on streamlines 

Another implementation of the multiple inlet/exit feature is the Princeton VPFR validation exercise presented in 
Section 4.2 where flow reactor injector inlets are shown in Figure A. 18. 

3 . 3.3 Parallel Chemistry Subroutines 

A high fidelity model of a coupled fluid and chemistry problem can require a large amount of computational 
resources. For a typical problem, the chemistry solution requires about 90% of the computational time. The 
desire to improve accuracy by resolving geometries with more detailed grids and through the use of broader 
chemical mechanisms, as well as, the desire to increase the turn-around time for high fidelity models prompted 
the exploration of methods to improve execution time. The structure of the CNEWT code, with the explicit 
operator-splitting algorithm in which the fluid and chemistry solutions are decoupled lends itself to a 
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Figure 3.4 Run time stastics for parallel CNEWT 

The solutions were also checked for accuracy by comparing a few maximum and minimum flow and species 
variables for the original code and the parallel version across several nodes. The fluid solutions for the original 
and parallel version (any number of processors) of the code were identical. The difference between the original 
and parallel version of the code with regard to the chemistry solution was found to be less than 0.2% in all of 
the species quantities checked. Also, the chemistry solutions from the parallel code were identical when the 
number of nodes was varied. Further instructions on how to compile and run the parallel version of the code are 
given in the README.TXT file included with the CNEWT code. 

3.3.4 Species and Temperature Wake Model 

The original CNEWT code employs a wake model which allowed the specification of a rotating, 
circumferentially varying inlet profile for total pressure. The main purpose of this wake model was to allow the 
specification of a pressure deficit associated with the wake(s) of an upstream blade or blade row without the 
need to do a full unsteady calculation with both blade rows using, for example, a sliding plane technique. Thus 
the code could be run in unsteady mode, with a fixed timestep, and the unsteady effects of the wakes could be 
captured using a single blade row. 

This wake model was extended to include total temperature and chemical species mass fractions. With these 
modifications, the capability of the code was increased by allowing a more accurate specification of the inlet 
conditions and enabling the investigation of the unsteady interaction of chemical and temperature non- 
uniformities with downstream blade rows. 
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3.4 Chapter Summary 


The purpose of this chapter was to present a brief overview of the modeling tools and discuss the improvements 
made to these tools during the course of this research. The existing modeling tools consist of the PRE grid 
generation code, PROCESS grid pre-conditioner code, POST grid post-processor code, CNEWT and 
CALCHEM flow-chemistry solvers, and CNEWTVS6 flow visualization code. The improvements made to the 
existing tools allow for more complicated geometries, multiple inlet and exits, improved execution time for 
solutions with chemistry, and a specification of a pressure, temperature, and species concentration deficits 
associated with a wake of an upstream blade or blade row. 

In summary, the following is a list of all relevant points presented in this chapter. 

• CNEWT and CALCHEM have been useful modeling tools for studying intra-engine trace species flow- 
chemistry problems. CNEWT is a code created by Lukachko [1] which was built upon the NEWT 
turbomachinery CFD code from Dawes [51]. CALCHEM is a simple, ID version of CNEWT. 

• Other supporting codes necessary for turbine chemistry modeling are PRE, PROCESS, POST, and 
CNEWTVS6. 

• For more details on the structure and numerical mechanics of the codes, refer to “Research on the Science 
And Politics of the Atmospheric Effects of Aviation Debate’ by Lukachko [1] as the primary source of 
information. 


• In order to provide further insight into trace species characterization, improvements were made to existing 

modeling tools which include: 

• A new CAD solid model based grid generator which allows for the simulation of more complicated, 
arbitrary geometries. The output is a truly unstructured grid. The tools were in a developmental stage 
and a few difficulties were encountered with the translation step and numerical issues in the volume 
mesher. The VPFR was successfully meshed with the new grid generator. 

• Addition of multiple inlet/exit capability which allows for simulation of mixing flows. A test case and 
the Princeton VPFR were used to validate the multiple inlet/exit code modifications. 

• Addition of parallel chemistry subroutines. A test case showed a factor of 3.4 improvement in time on 
6 compute nodes. The accuracy of the parallel chemistry modifications was verified. 

• Implemented and tested a wake model boundary condition for temperature and chemical species which 
allows for the investigation of unsteady effects including flow-chemistry interactions. 
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4 Validation 


Validation of the modeling tools is essential to provide confidence in the results obtained using the modeling 
methodology and verify new features or modifications. A validation exercise using the Princeton Variable 
Pressure Flow Reactor (VPFR) was pursued in order to prove the accuracy of the analysis tools and modeling 
methodology, in particular the new grid generator and multiple inlet/exit modifications. The chemical 
mechanism is a fundamental element underlying the intra-engine chemistry model. Thus the selection of an 
appropriate chemical mechanism for intra-engine chemistry modeling must also be justified. The convergence 
criteria and periodic boundary conditions are also important aspects of the model that must be proved accurate. 
All of the validation exercise presented help add to the knowledge base available to improve the modeling of 
intra-engine post-combustor chemistry. 

4.1 Chapter Overview 

The purpose of this chapter is to describe the validation exercises performed during the current research period. 
An attempt to benchmark the modeling tools through simulation of a flow reactor experiment is presented. 
Although it was not completely successful, it provided the opportunity to test several code improvements, 
investigate the limitations of the modeling tools, and analyze flow features of the reactor. A chemical 
mechanism study which investigated the effects of the chemical mechanism on intra-engine chemistry modeling 
and selected a mechanism for use in the NASA/DERA engine test modeling is detailed. This chapter also 
discusses a convergence criteria study which resulted in improved chemistry convergence indicators. Finally, a 
periodic boundary condition study is presented. 

In summary, the following is a list of the main points of the chapter: 

• An attempt was made to model a VPFR experiment to benchmark the accuracy of the modeling tools, 
however, a limitation to the flow solver prevented this ultimate objective. The flow features from the 
simulation did not match those expected from experimental data and empirical correlations. 

• The VPFR validation exercise proved useful in validating some improvements to the modeling tools, 
namely the new fully unstructured grid generator and the multiple inlet/exit improvement. 
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This validation effort was conducted in collaboration with Princeton University. Professor Fred Dryer and 
Mark Mueller made important contributions to the work presented in this section. 


4.2.1 Objectives 

The VPFR is ideally suited to evaluate the capability of the modeling tools to represent many flow features of 
interest , mixing, and passively reacting flow. The specific objectives of the VPFR modeling effort were to: 

• Model a VPFR reacting flow experiment in order to benchmark the accuracy of the modeling tools 

• Identify the capabilities and limitations of the analysis tools, specifically with regard to characterizing the 
flow features, mixing, and chemistry 

• Help expand the modeling techniques by evaluating several new code features, in particular, the new grid 
generator, multiple inlet and exit planes, and parallel chemistry subroutines 

• Help establish the modeling methodology applicable to studying intra-engine trace chemistry 

• Provide an analysis that can identify potential improvements to the flow reactor facilities or experimental 
protocol 

The approach was to model a published flow reactor experiment and compare the simulation data to the 
experimental data. The plan was to first simulate the flow-only, then flow with inert species, and finally 
passively reacting flow. In the end, the objectives of the flow reactor were reduced due to the restriction of the 
code to simulation of compressible flows. 


4.2.2 Summary of Results 

Limitations of the flow solver hampered the ultimate objective of validating the modeling tools through a 
successful reacting flow simulation of the Princeton VPFR. Although the validation exercise was not 
completed, many code improvements were tested, a considerable amount was learned about the limitations of 
the tools, and some potential issues with the flow reactor design were identified. 

The new fully unstructured CAD based grid generator was successfully used to mesh the complex VPFR 
geometry. At the time of the work in this thesis, the new grid generation tools were in a developmental state. A 
significant amount of geometry repair was required to obtain a surface mesh and the volume mesh routine had 
trouble with large aspect ratio geometries. The multiple inlet/exit improvement was qualitatively validated. It 
was found that the VPFR grid had to be initialized to a zero initial guess for the velocity to achieve a solution 
due to the complicated geometry. 

Several experimental conditions were considered when applying the modeling tools. The assumptions for 
passively reacting flow, steady or unsteady flow, and compressible or incompressible flow need to be addressed 
prior to embarking on a modeling effort. Two cases were presented to illustrate the major findings of the VPFR 
validation efforts, a low Re case (Run 1) and a high Re case (Run 12). The results were evaluated by 
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4.3 Chemical Mechanism Study 


An important element of the modeling methodology is the chemical mechanism, which describes the reactions 
among the species set. Previous work [1] on intra-engine chemistry modeling used a chemical mechanism that 
was published in [28] and was developed by Aerodyne Research, Inc. [67], The mechanism was created by 
doing a survey of the chemical kinetic literature and piecing together various reactions and sub mechanisms 
which were of interest to intra-engine trace species chemistry. The main species of interest to intra-engine 
modeling are HN0 2 , HNO,, S0 2 , S0 3 , H 2 S0 4 , OH, H 2 0 2 , H0 2 , O, and NO, [67]. However, the task of 
assembling literature data into a mechanism is not easy since much of the chemical kinetics data available is not 
applicable to the ranges of pressure and temperature experienced in gas turbine engines. As the current research 
was being conducted, a new chemical mechanism was developed at Princeton University based on flow reactor 
experimental data. This section contains a study of chemical mechanisms which was conducted to validate the 
prior intra-engine chemistry work [1] and possibly improve the chemical mechanism used in the modeling 
methodology. 

4 . 3.1 Objectives 

The main objective of this study is to investigate the sensitivity of the modeling methodology to the selection of 
chemical mechanism. In essence, the effects of chemical kinetics on trace species evolution were studied. The 
results of this study were used to validate prior work and to select the chemical mechanism for the 
NASA/DERA engine test modeling. The selected mechanism is also recommended for future intra-engine trace 
chemistry modeling. 

4 . 3.2 Approach 

The approach was to repeat prior calculations done by Lukachko [1] using variations of a mechanism which 
was based on recent experimental data. A comparative study was done on the following four turbine chemistry 
mechanisms: 

• Lukachko et al. (1998) [28] (26 species and 74 reactions) 

• Glarborg et al. ( 1996) [24] (56 species and 169 reactions) 

• Mueller et al. (2000) - full [23] (56 species and 169 reactions) 

• Mueller et al. (2000) - truncated [23] (29 species and 73 reactions) 

The Glarborg et al. (1996) and Mueller et al. (2000) mechanisms were based on experimental data and contain 
similar species subsets as Lukachko et al. (1998). The Mueller et al. (2000) mechanism is based on the 
Glarborg et al. (1996) mechanism with modifications in the SO, chemistry at higher pressures. The Mueller et 
al. (2000) mechanism is split into a full and truncated version. The full mechanism contains the complete set of 
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CASE l - NON-EQUILIBRIUM 


CO. OH. and O 



Figure 4.2 Chemical mechanism study for Case 1 



Figure 4.3 Chemical mechanism study for Case 2 
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Table 4.1 Summary of results from chemical mechanism study case 1 



Table 4.2 Summary of results from chemical mechanism study case 2 



Table 4.3 Summary of results from chemical mechanism study case 8 
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Figure 4.6 Example result of 2D calculation on Cambridge No 2 turbine rotor with Lukachko et ai (1998) 

mechanism 


4.3.5 Conclusions 

Overall, the ID ASE engine cycle and 2D Cambridge No 2 turbine rotor blade row calculations both show that 
the species concentration trends are the same and in particular the range of SO x conversion for a given case was 
2.7% to 7.3% for all the mechanisms. This result would imply that the chemical mechanism does have some 
impact on the results of the modeling effort. However, as long as the key reactions are included in the 
mechanism the results are in fairly good agreement. Therefore, the prior calculations [1] were substantiated by 
new experimentally-derived kinetic data. 

The modifications to the SO x chemistry at high pressures in the Mueller et al. (2000) mechanisms were based 
on experimental data and are considered an improvement with respect to the turbine chemistry mechanism 
previously employed. Therefore, the Mueller et al. (2000)-truncated mechanism was selected as the best option 
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Convergence Indicator vs. Number of Iterations 
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Figure 4.7 Original CNEWT fluid and chemistry convergence indicators 

From the Figure, it appears that the maximum nodal species residual converged by 600 iterations and the 
combined average converged by 1,000 iterations for the test case. The fluid solution was already converged at 
the start of this solution. It was suspected that the trace species could still be changing after the average and 
maximum criteria had converged since the trace species would play a minor roll in the average calculation and 
would not necessarily be the maximum residual. In other words, the average and maximum quantities could 
mask the actual changes in the chemistry solution since large portions of the domain could be converged while 
small portions could be changing significantly. Several potentially improved convergence criteria were 
proposed, including monitoring the maximum nodal mass fraction for individual species of interest, the RMS 
average for individual species of interest, or the maximum nodal residual divided by that species mass fraction. 


4 . 4.3 Results 

First, a straightforward visual approach was employed to determine the actual convergence for a few trace 
species of interest. The chemical variables were re-initialized for the test case and it was set-up to re-run with a 
solution file dumped every 100 iterations. Time history plots for the species evolution of O, OH, and S0 3 mass 
fraction were generated for 5,000 iterations. As an example, the mass fraction of O is plotted at several points 
during the solution in Figure 4.8. Also, animations for all three species are included in Appendix F. 


79 



8 






Indicator Conv*r9*nc« Indicator 


Convergence Indicator vs. Number of Iterations 



Convergence Indicator vs. Number of Iterations 



8 8 


S | '§ '§ § i 8 '§ § '§§ 8 § § I 


Number of iteration* 


■Average of 
AJ Species 
Residua)* 


Max Nodal 

Species 

ResiduaJ 


-OH 


I 


0 


— S03 


- MONO 


H2S04 


Figure 4.10 Comparison of original and new chemistry convergence indicators for major species (top) 

and trace species (bottom) 
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4.5.2 Results 


The results of the periodic boundary condition study are summarized in Table 4.5. Excluding V r , the maximum 
error relative to the value at the 0% location was found to be 1.6%. The psuedo-2D grid had only two cells in 
the radial direction, and thus V r had a very small value and was not resolved. 


Summary of Results 
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Table 4.5 Summary of results from periodic boundary condition study 


The detailed results for station 3 are given in Figure 4.12 as an example. The blue circle symbols in the figure 
are the 0%, 33%. 66%, and 100% circumferential locations. 


85 












































































































































































































































presented. Although it was not completely successful, it provided the opportunity to test several code 
improvements, investigate the limitations of the modeling tools, and analyze flow features of the reactor. A 
chemical mechanism study which investigated the effects of the chemical mechanism on intra-engine chemistry 
modeling and selected a mechanism for use in the NASA/DERA engine test modeling was detailed. This 
chapter also discussed a convergence criteria study which resulted in improved chemistry convergence 
indicators. Finally, a periodic boundary condition study was presented. 

In summary, the following is a list of all relevant points presented in this chapter: 

• An attempt was made to model a VPFR experiment to benchmark the accuracy of the modeling tools, 
however, a limitation to the flow solver prevented this ultimate objective. The flow features from the 
simulation did not correlate well with those expected from experimental and empirical data. 

• The VPFR validation exercise proved useful in validating some improvements to the modeling tools. 

• The new fully unstructured grid generator was successfully used to mesh the VPFR geometry 

• The multiple inlet/exit improvement was qualitatively validated. 

• The modeling tools are limited to simulations of compressible flow and should be restricted to situations 
when the bulk of the flow field is above a Mach number of 0.2, an estimate for the compressibility limit of 
the code. 

• The Princeton VPFR may operate with large scale unsteady stall. 

• A comparative study was performed on four chemical mechanisms by repeating a ID and 2D computation 
from prior work. 

• The species concentration trends were the same for all mechanisms studied. 

• The SO* conversion ranged from 2.7% to 7.3% for a given case, which implies that the selection of 
chemical mechanism does impact intra-engine chemistry modeling. 

• Accurately representing key reactions in the mechanism produced reasonable results, thus prior work 
using Lukachko et al. ( 1998) was supported. 

• The Mueller et al. (2000)-truncated mechanism was selected for use intra-engine chemistry modeling. 
It is based on experimental data and is within 1% of the Mueller et al. (2000)-full mechanism. 

• The original chemistry convergence criteria did not represent convergence of trace species well. They were 
improved by using the RMS average of individual species residuals for the species of interest to determine 
solution convergence. 

• Periodic boundary conditions did not introduce any significant error into the computation, selected fluid 
and species variables at the boundary were within 1.6% of their corresponding value with no boundary 
present. 


87 



5 Modeling the NASA/DERA Engine 
Test 


This chapter presents the first complete analysis which applies the intra-engine trace chemistry modeling 
methodology to a real engine. The goal of the simulation is to characterize the evolution of trace species in the 
post combustor flow path in order to provide engine exit plane emissions predictions, help guide the engine test 
plan, and research a few physical mechanisms thought to affect trace species chemistry. 

5.1 Chapter Overview 

The purpose of this chapter is to describe the modeling effort to simulate trace species evolution in the post- 
combustor gas path of an engine in support of the NASA/DERA engine test. The objective of the test is to 
provide measurements to help characterize aviation emissions and benchmark the models. The model is used to 
inform the test effort by providing pre-test predictions which were used to help formulate the test plan and to 
research a few fundamental mechanisms thought to influence trace species chemistry. 

Three conditions are investigated: cruise, max power, and non-uniform max power. Details of the grid 
generation, fluid boundary condition specification, and species initial conditions specification are discussed. A 
time scale analysis is used to guide the modeling efforts. Two types of high fidelity models are used tor the 
HPT1; a mixed-out case and a wake model case. Several properties of the solution are discussed, including: 
convergence, overall flow features, grid resolution, deviation angle, and flow-through times. A sample set of 
detailed flow-chemistry results from the high fidelity model is presented. Low fidelity models are used for 
HPT1 and HPT1 exit to nozzle exit. Low fidelity models for averaged and high and low temperature profiles in 
the HPT 1 are compared to the high fidelity average. A HPT1 exit to nozzle exit profile is used to obtain a 
prediction of engine exit gas composition. A sample set of flow-chemistry results from the low fidelity model 
is presented. 

This chapter also highlights several implications of the modeling results. Specifically, engine operating 
conditions, multi-dimensional non-uniformities, and the unsteady interaction of non-uniformities with 
downstream stations were all found to influence trace species evolution. They can all impact temperature. 
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• The unsteady interaction non-uniformities with downstream blade-rows has a significant impact on intra- 
engine post-combustor chemistry. 

• Local regions of flow or species non-uniformity persist through the subsequent blade rows. 

Section 5.2 outlines the objectives of the engine test and modeling effort, provides some background 
information related to the engine test, and lists the modeling scenarios considered. Section 5.3 gives a time 
scale analysis for the engine test and gives an outline for the presentation of the details of the modeling effort. 
Section 5.4 details the set up and results of the high fidelity analysis. Section 5.5 details the set up and results 
of the low fidelity analysis. Section 5.6 describes the implications of the results from modeling effort. Section 
5.7 is the chapter summary. 

5.2 Objectives, Background, And Modeling Scenarios 

NASA/DERA Engine Test Objectives 

One of the project objectives outlined in Section 1.4 was a validation of engine modeling results through 
measurement studies. A collaboration of many research groups consisting mainly of teams from DERA, 
NASA, MIT/ARI, and UMR have formulated an engine test to address several emissions related phenomena. 
The main objective of the measurement campaign is to carry out ground-based simulation and measurement of 
chemical and physical processes relevant to aviation emissions characterization [8]. The NASA/DERA engine 
test aims to measure the emissions of a typical engine, evaluate the effects of various fuels and operating 
conditions on engine emissions, develop and evaluate various species measurement techniques, investigate 
intra-engine species evolution, and provide inputs to emissions models through characterizing the composition 
of the gas at the combustor and nozzle planes. 

Modeling Objectives 

The NASA/DERA engine test provides a unique opportunity to demonstrate the modeling methodology 
described in Section 2. The models were used to inform the test efforts by providing pre-test predictions for 
intra-engine species evolution. The pre-test predictions assisted in formulating the test plan by helping to 
ensure that operating conditions are selected so as to provide useful information about the intra-engine 
environment. The pre-test predictions are also useful in directing the measurement efforts by indicating the 
important species to monitor and providing estimates for the instrumentation requirements. Modeling of the 
engine test also provides the opportunity to benchmark the accuracy of the modeling methodology against 
measured data. A matched set (i.e. same engine and operating condition) of combustor and nozzle exit plane 
species data is required to evaluate the predictive capability of the intra-engine modeling tools. To date no such 
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Modeling Scenarios 


The engine modeling scenarios that were simulated in the context of this thesis were based on the nominal 
operating points for cruise and max power. Table 5.2 overviews the three simulations conducted which cover a 
range of combustor exit temperatures pertinent to many modern civil aircraft engines. 



Summar 

1/ of NASA/DERA Enc 

line Test Conditions 

Label 

Cruise 

Max power 

Non-uniform max power 

Description 

Nominal cruise 
power setting 

Nominal maximum 
power setting 

Maximum power setting with 
a combustor exit temperature 
profile (nominally 1 ,600K) 

Enaine speed (rpm) 

7,821 

8,864 

8,864 

Ts at combustor exit (K) 

1,123 

1,357 

1 ,430-1 ,840 (average-1 ,600) 

Ps at combustor exit (Pa) 

687,000 

1 ,594,000 

1 ,594,000 

Ts at HPT1 exit (K) 

983 

1,186 

1,186 

Ps at HPT1 exit (Pa) 

428,000 

986,000 

986,000 

HPT1 rotor metal temperature (K) 

826 

999 

999 

HPT1 stator metal temperature (K) 

791 

956 

956 


Table 5.2 Summary of conditions modeled for NASA/DERA engine test [34] 


These three test conditions provide a basis for investigating the effects of operating conditions on trace species 
evolution. As discussed previously in this section, the modeling strategy used for the NASA/DERA engine 
simulations called for high fidelity modeling of the first stage of the high pressure turbine followed by low 
fidelity modeling for the remainder of the post-combustor flow path. These computational simulations were 
extended to investigate the effects of several physical phenomena suspected to influence the evolution of trace 
species, namely circumferential non-uniformities and the unsteady interaction of the non-uniformities on 
downstream stations. 


Two of the simulated conditions were used to investigate two types of circumferential non-uniformity. The 
HPT1 rotor calculation for the max power and non-uniform max power conditions were run both using a steady 
mixed-out initial condition and an unsteady wake model profile initial condition to investigate the effects of 
unsteady interaction of the non-uniformities on downstream stations (see Section 3.3.4 for more details on the 
wake model). Figure 5. 1 gives a overview of all the modeling scenarios explored. 


93 



The profile was clocked such that the peak temperature, corresponding approximately to the injector locations, 
was aligned with a blade. This was considered a worst case scenario since the temperature difference between 
the combustor gas and the cooled blade would be greatest. Given the number of injectors and stator blades it is 
likely that this would occur somewhere within the circumference. The derived profile was then fit to the three 
NGV blade passages as shown in Figure 5.3 (0 to 18 degrees). 


Combustor Exit Temperature Profile 
(max and min not defined, shifted 1.5 deg) 



Figure 5.3 Combustor temperature profile used as initial condition for non-uniform max power case 

Species of importance to evaluating the SO, and NO y evolution are SO3, SO2, H2SO4, O, OH, NO, NO2. NOj 
and HONO. Several instrumentation teams will measure a multitude of primary and trace combustion products 
such as 0 2 , C0 2 , CO, H 2 0 , unburned hydrocarbons, and soot particles. Initially, measurements were focused on 
S0 3 evolution, however, pre-test simulations showed that majority of the S0 2 conversion was completed in the 
combustor for the operating conditions considered (more so for the cruise condition). These results suggested 
concentrating the test efforts on the highest power setting attainable such that the combustor exit temperature is 
highest and the magnitude of change is greatest along the post-combustor gas path (i.e. between the 
measurement stations). Also, the focus was broadened to include HONO evolution which exhibits a greater 
change between the measurement stations for the planned test conditions. Furthermore, since NO y chemistry 
remains active over this range of power settings it was still advantageous to explore multiple operating points. 
These modifications to the test plan would help to magnify the changes in HONO concentration through the 
measurement stations and provide the potential for comparison between two operating conditions. 
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Integration along the path from the combustor exit to the HPT1 exit indicates that there will be less than 4% 
sulfur conversion in that stage. Integration along the path subsequent to this conversion, the severity parameter 
suggests several orders of magnitude less potential for sulfur conversion after the HPT1. Using these findings it 
was decided to use a high fidelity 2D model for the HPT1 combined with a low fidelity ID model from the 
HPT1 exit to the nozzle exit to perform the pre-test modeling for the NASA/DERA engine test. 

5.3.2 Presentation of Modeling Effort 

A complete set of all the specific results accumulated during the course of this research would be overwhelming 
and not add much to the value to this thesis. In order to keep the write-up concise and not to obscure the key 
findings, the specific results were divided into two main parts with two subsections each. The two parts are: 
Section 5.4 High Fidelity Modeling and Section 5.5 Low Fidelity Modeling. The subsections are Set Up and 
Results. Finally, Section 5.6 Implications of Results contains an analysis of the results with regard to their 
impact on the research of intra-engine trace species evolution. The presentation layout gives a representative 
case for each distinct case without listing or plotting every chemical species or fluid variable for each operating 
condition and analysis technique. A diagram showing the important topics in the presentation layout is given in 
Figure 5.5. 



Figure 5.5 Presentation layout of NASA/DERA engine test modeling 
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HPT1 NGV 

Number of Blades = 60 
I. J, K = 25, 80, 3 

Cells = 23,160 

Chord = 31 .44 mm 

Inlet = 10.95 mm 

Axial Chord = 23.09 mm 
Exit = 1 7.82 mm 


Figure 5.6 Single pitch HPT1 NGV grid 

Figure 5.7 shows the rotor grid. The inlet length for the rotor was set as half of the distance from the NGV 
trailing edge to the rotor leading edge. The inlet conditions for the high fidelity rotor calculations were 
extracted from the NGV solution at that plane. Again, the exit length was extended roughly 0.8 axial chords 
downstream. 
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and pressure was derived for each of the three operating conditions investigated. The important input 
parameters used in the initial conditions code are listed in the table below: 



Ctiemistrv Initial Condition Code Inputs 

Condition 

Cruise 

Max power 

Non-uniform max power 

T (K) 

1,122.8 

1,356.6 

1 ,440-1 ,830 

P (atm) 

6.78 

15.732 

15.732 

0 (equiv ratio) 

0.35 

0.35 

0.35 1 

n (comb eft) 

0.9988 

0.99988 

0.99988 

CO/H2/HC ratio 

100/0/0 

100/0/0 

100/0/0 

EI(NOx) (g/kg fuel) 

21.0 

21.0 

21.0 

NO/NOx ratio 

0.99 

0.99 

0.99 

EI(S) (g/kg fuel) 

0.5 

0.5 

0.5 

S03/SOx ratio 

0.06 

0.06 

various equilibrium 
(about 0.03 average) 

EI(CO) (g/kg fuel) 

5.09 

0.51 

0.51 


fable 5.3 Chemistry inlet condition code inputs for NASA/DERA engine test 


The main differences in the chemistry inlet specification technique relative to prior work in [1] and the 
Cambridge No 2 turbine rotor calculations discussed previously are that the CO/H 2 /HC ratio was changed from 
78/22/0 to 100/0/0, the NO/NO x ratio was changed from 0.90 to 0.99, and a kinetics calculation was used to 
resolve rapid numerical changes in the initial conditions due to inconsistencies in the resulting speciations. Due 
to lack of emissions data, the ETs represent a best estimate for an engine typical of this make and era. The 
EI(NO x ) was adjusted to obtain approximately 300 ppm of NO* to coincide with a rough visual average of NO x 
contours of a similar gas turbine engine (see Figure 5.8). The EI(CO) was set by adjusting the combustion 
efficiency to achieve a value close to one interpolated from an emissions data base for a 1970’s era engine with 
a pressure ratio of about 16 at the appropriate power setting [1 1]. The selected EI(NO) was also similar to the 
ones listed in the data base. The SCtySO* was fixed at 6% for the cruise and max power setting based on 
experience from prior work [30] which showed that sulfur conversion in the combustor was limited. In 
contrast, the equilibrium level for these two cases resulted in a S0 3 /SO x ratio of about 38%. The temperature 
for the non-uniform max power condition spanned 1,440-1,830K which gave an equilibrium S0 3 /SO x ranging 
from about 1-8%, thus the equlibrium sulfur ratio was used in this case. 
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5 . 4 . 2.1 Calculation Description And Execution 


Overview Of Analysis Types For Rotor Calculation 

Figure 5.9 depicts the two types of high fidelity analysis were used for the HPT1 rotor in the current research. 
The first analysis uses the mass-averaged fluid and species quantities at the NGV exit plane as inlet conditions 
to the rotor domain and is called the “mixed-out case ” The second analysis uses the wake model described in 
Section 3.3.4 to apply a rotating, circumferentially-varying fluid and species profile from the NGV exit plane as 
the inlet condition to the rotor domain and is called the “wake model case/’ These two techniques were used to 
evaluate the effects of non-uniformities on downstream blade rows. 



Figure 5.9 Two high fidelity analysis types, mixed-out (top) and wake model (bottom) 
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The unsteady HPT1 rotor solutions were deemed to be converged when the convergence criteria became 
periodic about a constant value. Figure 5.12 and Figure 5.13 show example convergence histories for the HPT1 
rotor calculation at the max power conditions. The unsteady solutions converge faster, although to higher 
residuals than the steady solutions. The higher residuals are due the need to resolve the fluid and chemistry 
variables in time (with a finite time-step) as well as the spatial domain. 
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An full unsteady solution consists of the set of one complete cycle after convergence. Figure 5.14 shows the 
locations of a set of solution output files for an example unsteady rotor run at max power. In this case, ten 
evenly distributed locations in time were chosen for output over 275 iterations. All later runs spanned a 
minimum of 400 iterations per cycle for better resolution and typically consisted of 10 to 15 output locations per 
cycle. 


Convergence Indicator vs. Number of Iterations 



Number of Iterations 


Figure 5.14 An example set of solution output files making up one cycle of a full unsteady solution 

S.4.2.2 Flow Features And Analysis 

Overall 

In general, all of the solutions obtained in the high fidelity modeling phase were resonable. The features of the 
flow field and values of the fluid state variables were as expected. Table 5.4 contains a summary of the fluid 
quantities from the velocity triangle analysis and the high fidelity modeling. 
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analysis, the deviation angle is reasonable, and the inlet/exit conditions are as specified. In general, all of the 
solutions from the high fidelity models were probed in a similar manner and the flow results were found to be 
similar. 

Trailing Edge Resolution 

One peculiar feature of the flow field was a very small region near the trailing edge which had a significantly 
higher Mach number. Upon close inspection, it appears that the high speed flow is caused by the tendancy of 
the fluid to accelerate around the TE from the pressure side to the suction side of the blade prior to the point of 
separation. The phenomena is magnified by the faceted shape of the TE and the limited grid resolution. This 
phenomena was much less pronounced in the prior calculations of the Cambridge No 2 turbine rotor because the 
TE was modeled as squared off which fixed the point of separation at the corners. Plots of the Mach number 
and total pressure for both cases are shown in Figure 5.16 for comparison. By modeling the true curvature of 
the TE, the point of separation and flow field in this region are not well resolved for the given grid spacing. The 
effect of this deficiency in the representation of the flow field is likely to be minimal on chemical species 
evolution. Two possiblilities to improve the situation would be to cut ( square off ) the TE or to enhance the 
mesh in this region, however, both methods have drawbacks. Cutting the TE is in essence fixing the point of 
separation and thus dictates the shape of the wake. Reducing the initial grid spacing would start to skew the 
cells in that region and simply enhancing the grid using the boundary layer technique in POST will increase the 
number of cells drastically. 
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Figure 5.17 Streamlines near the trailing edge used to estimate the deviation angle 

The deviation measured from the high fidelity model was between the estimate from Carter’s Rule (7.5°) and 
the assumed value (2°). As mentioned, the flow parameters (including velocity vectors) from the high fidelity 
models at station B were used as input to the rotor high fidelity models, therefore, the assumed deviation of 2° 
was considered adequate to make the initial estimate of fluid boundary conditions. 

Blade-Row Flow-through Time 

The flow-through time for the HPT1 blade rows at cruise condition was estimated by integrating the total 
velocity along three streamlines as follows: 


t = 



(5.23) 


To give a feel for the range of flow-through times for various parcels of fluid, streamlines were chosen near the 
pressure and suction sides of the blade as well as one at mid passage. Figure 5.18 shows the streamlines and 
velocity profiles used for the calculation. 
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rotor leading edge (at X = 0.04515 m in the NGV grid). As an example, the mass-averaged total pressure was 
computed as follows: 


m 



(5.24) 


The column labeled “HPT1 Rotor Exit” is divided into two parts. The “mixed-out steady” column is the mass- 
averaged species concentration at the plane midway between the rotor trailing edge and the leading edge from 
the NGV on the following turbine stage. HPT2, for a mixed-out case (at X = 0.0106 m in the rotor grid). In the 
mixed-out case, the inlet condition for the rotor was simply specified as a constant value equal to the mass- 
averaged species concentration from the “HPT1 NGV Exit” for each species and also the mass-averaged 
temperature and pressure. In effect, all of the quantities were instantaneously mixed-out at the rotor inlet plane 
and thus all non-uniformities were removed. The “time ave. unsteady” column is the time-average of the mass- 
averaged species concentration at the same plane for the unsteady wake model case. As an example, the time- 
averaged mass-averaged total pressure was computed as follows: 



(5.25) 


In the unsteady wake model case, the actual 2D non-uniform species, temperature, and pressure profiles were 
applied to the rotor inlet using the wake model. The wake model which was discussed in detail in Section 3.3.4 
preserves the upstream circumferential non-uniformities. Some example profiles for the max power case are 
given in the graphs on the third row of results in Figure 5.19 and Figure 5.20. The blue symbols on the graphs 
represent the data points input into the wake model. 


Figure 5.19 and Figure 5.20 give example results for total temperature, total pressure, and a few selected species 
mass fractions from the high fidelity modeling for the NGV and mixed-out rotor cases. The first row is a 
spanwise view of the contour plot for the particular quantity where the black vertical line represents the blade- 
row exit plane. The second row is an axial view contour plot on that plane. The third row of the NGV plots is a 
plot of species mass fraction versus circumferential location at the center of that plane. And the last row lists 
the range for the contour plot and the mass averaged value on that plane. 


The total temperature deficit in the NGV wake is about 9%. Also, notice that the SO 3 and HONO mass fraction 
profiles are inverse to the temperature wake profile, suggesting that the temperature gradients have the strongest 
influence on their local concentration and that blade cooling enhances their production. 
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Figure 5.20 SO 3 , HONO, and OH at NGV and rotor exit for max power case 
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The wake model rotor cases are unsteady (i.e. time varying) solutions due to the rotation of the upstream stator 
wakes and combustor profile relative to the rotor blades. Therefore, the specific results consist of a series of 
plots equally spaced out over a time interval equal to one cycle. An example result of the wake model case is 
shown in Figure 5.22 for the non-uniform max-power case. A similar result for the max power case is given in 
Appendix C. The time dependent results can also be viewed in the animations contained in the attached media 
detailed in Appendix E. 

From the figure, it is evident that the localized regions of non-uniform S0 3 mass fraction persist downstream 
past the exit of the blade row. Also, notice that the non-uniformity from the combustor temperature profile is 
much larger than the non-uniformity from the wakes of the cooled blades (i.e. more of the flow is subject to the 
larger scale combustor non-uniformity). 
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A summary of all the key results is given in Table 5.5, Table 5.6, and Table 5.7 for the cruise, max power, and 
non-uniform max power conditions, respectively. The species quantities are given in mole fraction. The 
unsteady wake model case and HPT1 exit to nozzle exit was not run for the cruise condition. Again, for a 
complete set of detailed results for the max power conditions see Appendix C. 


; 


Summw 

ryqfSt-DSbnulation 

Results: CruM* C( 

mdlUon 

V 



Combustor Exit 

HPT1 NGV Exit 

HPT1 Rotor Exit 

Nozzle Exit 

Number 

Species 

Initial Condition 


mixed-out steady 

time ave. unsteady 


1 

CO 

1.21890E-04 

1.21772E-04 

1.21739E-04 

N/A 

N/A 

2 

C02 

4.76468E-02 




• N/A 

N/A 

3 

02 

1. 3321 7E-01 

1.33217E-01 

1.33217E-01 

N/A 

N/A 

4 

H20 

4.78257E-02 

4.78257E-02 

4.78257E-02 

N/A 

N/A 

5 

H 

1.68680E-10 

1.64263E-10 

1.52266E-10 

N/A 

N/A 

6 

0 

4.88767E-09 

4.92741 E-09 

4.53878E-09 

N/A 

N/A 

7 

OH 

3.681 09E-07 

3.16391 E-07 


N/A 

N/A 

8 

H02 

9.90554E-09 

8.70358E-09 

7.39922E-09 

N/A 

N/A 

9 

H202 

4.64661 E-10 

6. 84852 E-10 

7.97481 E-10 

N/A 

N/A 

10 

NO 

3.02737E-04 

3.02550E-04 

3.02462E-04 

N/A 

N/A 

11 

N02 

3.10783E-06 

3.20137E-06 

3.23082E-06 

N/A 

N/A 

12 

N03 

5.26957E-13 

4.23230E-13 

3.37662E-13 

N/A 

N/A 

13 

HNO 

5. 33968 E- 11 

1.38233E-10 

1.63320E-10 

N/A 

N/A 

14 

HONO 

5.3221 4E-08 

1.45864 E-07 

2.04203E-07 

N/A 

N/A 

15 

HN03 

1.50489E-10 

5.45819E-10 

8.45071 E-10 

N/A 

N/A 

16 

S02 

9.83418E-06 

9.83378E-06 

9.83355E-06 

N/A 

N/A 

17 

S03 

6.290 53E-07 

6.29292E-07 

6.29353E-07 

N/A 

N/A 

18 

IBB 


6.53196E-15 

4.32394E-15 

N/A 

N/A 

19 

HS03 

— H 1 II 1 1 III 1 

4.87219E-13 

6.07823E-13 

N/A 

N/A 

20 

HS02 

5.81476E-1 7 

7.31318E-17 

8.16752E-17 

N/A 

N/A 

21 

SO 

1.24025E-14 

7.82809E-15 

4.90596E-15 

N/A 

N/A 

22 

HCO 

2.09232E-16 

1.78302E-16 

1.45277E-16 

N/A 

N/A 

23 

H2 

5.76045E-10 

1.35492E-09 

1 .49002 E-09 

N/A 

N/A 

24 

N2 

7.70872E-01 

7.70873E-O1 

7.70872E-01 

N/A 

N/A 

25 

AR 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

N/A 

N/A 

26 

ESI 


O.OOOOOE+OO 

O.OOOOOE+OO 

N/A 

N/A 

27 

N 

O.OOOOOE+OO 

0.00000E+00 

O.OOOOOE+OO 

N/A 

N/A 

28 

H2S04 

2.7 1298 E-10 

4.32577E-10 

5.97855E-10 

N/A 

N/A 

29 

S 

8.93357E-26 

8.93358E-26 

8.93357E-26 

N/A 

N/A 

Note: All < 

quantities are species mole fractions 

(S03+H2S04)/S0x 


6.017002% 

6.017683% 

N/A 


A from comb exit 


0.002375% 

0.003057% 

N/A 


(HN02+HN03)/N0y 

■HBQQQTgj 

0.047862% 

0.067032% 

N/A 

N/A 

A from comb exit 


0.030415% 

0.049584% 

N/A 

N/A 

CO/COx 


0.255167% 

0.254920% 

0.254850% 

N/A 

N/A 

A from comb exit 

0.000000% 

-0.000247% 

-0.000317% 

N/A 

N/A 


Table 5.5 Summary of results from HPT1 high fidelity modeling at cruise condition 
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Summary of 2-0 Simulation Results: 

Non-Uniform Max Power Condition 


Number 

Species 

Combustor Exit 
Initial Condition 

HPT 1 NGV Exit 

HPT1 R 

mixed-out steady 

otor Exit 

time ave. unsteady 

Nozzle Exit 

1 

CO 

1.15329E-05 

8.351 63E-06 

7. 38 236 E-06 

7.79052E-06 

5.03788E-06 

2 

C02 


4.781 35E-02 

4.78145E-02 

4.781 09E-02 

4.78168E-02 

3 

02 

1 .341 60E-01 

1 34167E-01 

1.34172E-01 

1.34181 E-01 

1.34182E-01 

4 

H20 

4.77659E-02 

4.77844E-02 



4.78231 E-02 

5 

H 


3.13144E-08 

1.9071 IE-08 

2.65527E-08 


6 

0 

5.61601 E-06 

3.61282E-06 

2.78809E-06 

3.09000E-06 

2.26319E-09 

7 

OH 

1.20450E-04 

8.40250E-05 

6. 1 4268E-05 

6.45515E-05 

1.76425 E-07 

8 

H02 

6.76771 E-07 

4.431 11 E-07 

4.21 141 E-07 

4.71 472E-07 

1.12166E-09 

9 

H202 

3.81 1 33E-08 

4.82831 E-08 

1 .07656 E-07 

8. 13361 E-08 


10 

NO 

3.03484E-04 

3.02503E-04 

3.01528E-04 

2.98333E-04 

2.87829E-04 

11 

N02 

2.8821 9E-06 

3.83386E-06 

4. 75487 E-06 


1 .09828 E-05 

12 

N03 

7.33999E-12 

1.97938E-11 

6.09685E-1 1 

9.39380E-1 1 

5.28522E-13 

13 

HNO 

3.88249E-10 

3.50353E-10 

6.70221 E-10 

4.79792E-10 

3.15691 E-10 

14 

HONO 

2.62278E-08 

5.56272E-08 

1.09335E-07 

1.20774E-07 

7.47313E-06 

15 

HN03 

7.55054E-11 

2.47416E-10 

5.20103E-10 

1 .41 816E-09 

1 .07858 E-07 

16 

S02 

1.01206E-05 

1 .01041 E-05 

1.00939E-05 

1.01001 E-05 

1 .0031 1 E-05 

17 

S03 

3.41779E-07 

3.58292E-07 



4.29548E-07 

18 

HOSO 

3.87086E-12 

2.08593E-12 

1.28268E-12 

1.62738E-12 

1 .92871 E-16 

19 

HS03 

6.59378E-12 

6.38182E-12 

8.06975E-12 

6.88370E-12 

2.1 1402E-12 

20 

HS02 

3.77560E-15 

2.22706E-1 5 

2.10835E-15 

2.05852E-1 5 

7.21 925E-1 8 

21 

SO 

6.86930E-1 1 

3.42928E-1 1 

1.46333E-11 

2.15939E-1 1 

6.72918E-16 

22 

HCO 

6. 085 98 E- 15 

1.91539E-15 

1.19428E-15 

1 .2761 8 E- 15 

1 . 44639 E-1 9 

23 

H2 

8.1 4864 E-07 

5.79348E-07 

4.3491 6E-07 

5.26470E-07 

1 .09997E-07 

24 

N2 

7.69808E-01 

7.69809E-01 

7.6981 0E-01 

7.69761 E-01 

7.69810E-01 

25 

AR 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

26 

■sfteM 

O.OOOOOE+OO 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

27 

N 

4.01390E-12 

4.00640E-1 2 

4. 00640 E- 12 

3.8651 IE-1 2 

4.00640E-12 

28 

H2S04 

3.05830E-1 1 

4.19789E-11 



1.80619E-09 

29 

S 

4.1 1071 E- 16 

4.10394E-16 

4.10395E-16 

3.99506E-1 6 

4.10395E-16 

Note: All i 

quantities a 

re species mole fractions 


(S03+H2S04)/S0x 
A from comb exit 


■H3 

3.522683% 

0.255658% 

3.429555% 

0.162531% 

4.106721% 

0.839697% 

(HNQ2+HN03)/N0y 
A from comb exit 


0.018236% 

0.009651% 

0.035854% 

0.027269% 

0.039538% 

0.030953% 

2.474273% 

2.465688% 

CO/COx 

A from comb exit 


BBSS 

0.015437% 

-0.008679% 

0.016292% 

-0.007825% 

0.010535% 

-0.013582% 


Table 5.7 Summary of results from HPT1 high fidelity modeling at max power condition 


5.5 Low Fidelity Modeling 

5.5.1 Low Fidelity Modeling Set Up 

Low fidelity models were use to model both the HPT1 and the HPT1 exit to nozzle exit flow paths. Several 
steps were required to prepare for these runs. First, the ID temperature, pressure, and velocity profiles were 
extracted from the high fidelity models of the HPT1 NGV and rotor blade rows. A mass-averaged profile as 
well as profiles along several streamlines were investigated. Next, the available engine performance data was 
used to generate similar ID profiles for the HPT1 exit to nozzle exit. A set of input files to CALCHEM for the 
max power condition is given in Appendix D as an example. 
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Streamlines 


Representative low and high temperature streamlines were also selected at the exit of the NGV and rotor 
domains. Flow and species quantities were extracted along those streamlines for both the max power and non- 
uniform max power cases. The high temperature streamline was between the blades near the pressure surface 
and the low temperature streamline was in the boundary layer near the suction surface for both the NGV and 
rotor. Figure 5.24 shows the streamlines selected for the max power case. The corresponding high and low 
temperature ID profiles extracted from those streamlines are shown in Figure 5.25. The streamline profiles 
were used to perform a ID analysis which would be expected to show the range in magnitude of species change 
due to 2D non-uniformity. 




Figure 5.24 High and low temperature streamlines selected for ID profiles 
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HPT1 Exit to Nozzle Exit ID Profile (DERA Data) 



Xll*l Dt»t*nc* 


HPTt Exit to Nozzl# Exit ID Proflto 



Figure 5.26 HPT1 exit to nozzle exit ID profile, available data (top) and curve fit (bottom) 

Figure 5.27 shows the final normalized HPT1 exit to nozzle exit ID profiles used for the low fidelity modeling. 
The magenta lines approximate the locations of the inlet and exit planes for the HPT2, LPT1, and LPT2 turbine 
stages. 
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HPT1 and Engine Flow-Through Time 


The flow-through times ranged from 0.24 ms for the mass averaged case to 0.32 ms for the high temperature 
streamline of the HPT1 at max power conditions. For the ID analysis the flow-through time was computed by 
integrating the axial velocity along the axial distance rather than the total velocity along the actual streamline 
path. These estimated flow-through times for the HPT1 were consistent with those from the high fidelity 
modeling (see 5.4. 2.2). The flow-through time for the HPT1 exit to nozzle exit was 2.77 ms for the max power 
condition, giving an overall flow-through time for the post-combustor gas path of about 3 ms. 


5 .5.2.2 Flow-Chemistry Results 

A summary of the important chemical species results from HPT1 low fidelity modeling is given in Table 5.6 
and Table 5.7. In each table, the column labeled “Nozzle Exit” is the result of the ID low fidelity modeling of 
the HPT1 exit to nozzle exit using the mixed-out HPT1 results as input to the low fidelity models. 

The results from the low fidelity model were plotted as mass fraction versus time for three subsets of species. 
The plots are divided into two halves, the left graph shows the species evolution within the HPT1, and the right 
graph shows the species evolution from the HPT1 exit to the nozzle exit. Notice that the axis scales are 
different on both graphs. The species subsets were divided by family: SO* (SO 2 , SO 3 , H 2 SO 4 , HSO 3 . SO), NO y 
(NO, N0 2 , HONO, HN0 3 , HNO, N0 3 ), and CO, OH, and O. The plots contain the three ID calculations for the 
mass average, high temperature streamline, and low temperature streamline profiles, as well as, the mass- 
averaged species from the high fidelity modeling (labeled “2D Mass Average ). As an example, the NO y plots 
are given in Figure 5.28 and Figure 5.29 for the max power and non-uniform max power case, respectively. For 
a complete set of detailed results for the max power conditions see Appendix C. 
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DERA Engine Exit Plane Calculation For Non-Uniform Combuitoi Profile Max Power Case 
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Figure 5.29 Example low fidelity modeling results, NOy species evolution in the HPTl and from the 
HPTl exit to nozzle exit for the non-uniform max power condition 

NOTE: The legend on the chart lists the species in the proper sequence for the lines (top to bottom), however, 
the line style is not correct (dashed, dotted, solid, etc is wrong). 

5.6 Implications of Results 

5.6.1 Effects Of Engine Operating Conditions 

The three operating points studied for the NASA/DERA engine test simulation are typical of an older, but 
currently in-use engine. They can be used to assess the effects of engine operating conditions on trace species 
evolution and engine emissions. The cruise case has a combustor exit temperature of 1 120 K and pressure of 
6.8 atm, while for the max power case they are 1360 K and 15.7 atm, and for the non-uniform max power case 
the temperature is nominally 1600 K. It was found that S0 3 and HONO concentrations increased with 
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the mass averaged T, P, and V x profiles. There was a significant difference between the ID and 2D 
computations. The ID computation under-predicted the SO3 and HONO by 0.25% and 11.25% for the max 
power case and 2.7% and 21.3% for the non-uniform max power case, respectively. This trend is general and 
the difference is expected to depend on the magnitude of the temperature, pressure, or residence time 
discrepancies captured in the 2D model, as well as, the scale of those non-uniformities. 

Another way to look at these effects is to consider the different pressure, temperature, and velocity profiles 
encountered by fluid elements traveling along different streamlines. The species concentrations along a high 
and low temperature streamline were also plotted versus axial distance. The results from the two streamlines 
differ significantly from the averaged values, as well as each other. These results further support the 
conclusions from previous work by Lukachko [1] that multi-dimensional non-uniformities have a considerable 
influence on trace chemistry evolution which indicates that some important features can not be captured in ID 
models. Furthermore, other non-uniformities such as endwall boundary layers/cooling (similar to blade 
cooling), cooling air mass addition (addition of 0 2 ), regions of separation or re-circulation (large residence 
time), or afterburner duct geometry (large residence time) would also be expected to affect trace species 
evolution. 


5.6.3 Effects Of The Unsteady Interaction Of Non-Uniformities 
With Downstream Stations 

The high fidelity model for the rotor of the HPT1 was used to investigate the effects of the persistence of non- 
uniformities on subsequent blade rows. Two cases were investigated, focusing on temperature non- 
uniformities; the max power case which has smaller scale non-uniformities mainly caused by the wakes from 
the upstream cooled stator blades, and a non-uniform combustor exit temperature profile case which has a larger 
scale non-uniformity caused by combustion variations from the discrete fuel injector nozzles as well as the non- 
uniformity from the upstream cooled stator blades. Qualitatively, local regions of non-uniformity can be seen to 
persist through the downstream blade row in Figure 5.22. The mass-averaged exit plane species concentrations 
were time averaged over one cycle for the unsteady cases. The (S0 3 +H 2 S0 4 )/S0 x is 0.0005% higher and the 
(HN0 2 +HN0 3 )/N0 y is 0.007% higher at the rotor exit plane for the max power case with the effects of the 
unsteady non-uniformities relative to the mixed-out case. The low reactivity of the S0 3 chemistry for this 
condition makes it difficult to make a definitive statement about the unsteady effects of the non-uniformity on 
SO3, although, it appears to enhance the oxidation of both trace species. Including the unsteady non- 
uniformities increases the HONO ratio by 26% relative to the mixed-out case. 

The (S0 3 +H 2 S0 4 )/S0 x is 0.09% lower and the (HN0 2 +HN0 3 )/N0 y is 0.004% higher at the rotor exit plane for 
the unsteady non-uniform combustor exit temperature max power case. This case shows that the effect of the 
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• A first of its kind analysis was completed on the post-combustor gas path of the NASA/DERA engine to 
support an engine test aimed at characterizing aviation emissions. 

• The engine is typical of an older, but currently in-use civil aircraft engine. Three conditions were 
simulated: cruise, max power, and non-uniform max power. 

• Pre-test simulations were used to guide the test plan. They helped to ensure that the operating conditions 

selected provide useful information about intra-engine chemistry. Also, they helped direct the 

measurement efforts by indicating important species to monitor and giving estimates for the 
instrumentation requirements. 

• The pre-test simulations suggested testing at the highest combustor exit temperature attainable, as well as, a 
lower power setting. Also, SO3 and HOHO were selected as important species to monitor. 

• A time scale analysis indicated that most SO 3 production occurs in the HPT1 at max power conditions 
(approximately 4%). The analysis indicated that there is several orders of magnitude less sulfur conversion 
after the HPT1, and modeling resources were thus concentrated on the HPT1. 

• PRE was used to generate grids of the NGV and rotor. A velocity triangle analysis was used to derive 
necessary intra-stage fluid parameters for the high fidelity modeling. The species initial condition for each 
condition were derived using the technique outlined in Section 2.3.2. 

• Two types of high fidelity modeling were employed, a mixed-out case and a wake model case. All 
solutions were deemed converged using the criteria from Section 4.4. Also, the flow features such as 
deviation angle, streamlines, and inlet/exit conditions were checked to ensure the quality of each solution. 

• The boundary layers and trailing edge geometry were not fully resolved for the grids used, however, they 
are not expected to have a significant impact on the trace chemistry results. 

• HPT1 NGV blade-row flow-through time ranged from 0.18 to 0.26 ms and HPT1 rotor blade-row flow- 
through time ranged from 0.10 to 0.15 ms. The flow-through times computed from the low fidelity 
modeling were similar. 

• A summary of all high fidelity modeling results is contained in Table 5.5, Table 5.6, and Table 5.7. 
Examples of some detailed results of the high fidelity modeling can be found in Section 5. 4. 2. 3 or 
Appendix C. 
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• Local regions of flow or species non-uniformity persist through the subsequent blade rows. 

• With the small scale non-uniformities (wakes from cooled blades) in the max power case, the 
(S0 3 +H 2 S0 4 )/S0 x is 0.0005% higher and the (HN0 2 +HN0 3 )/N0 y is 0.007% higher at the rotor exit 
plane with the effects of the unsteady non-uniformities relative to the mixed-out case. In other words, 
including the unsteady non-uniformities increases the HONO ratio by 26% relative to the mixed-out 
case. 

• With the large scale non-uniformities (combustor temperature profile and wakes from cooled blades) 
in the non-uniform max power case, the (S0 3 +H 2 S0 4 )/S0 x is 0.09% lower and the 
(HN0 2 +HN0 3 )/NO y is 0.004% higher at the rotor exit plane with the effects of the unsteady non- 
uniformities relative to the mixed-out case. In other words, mixing out the non-uniformities over 
predicts the conversion of S0 3 by about 60%. 


135 



6 Thesis Summary And Contributions 


This chapter will present a summary of the research discussed in this thesis and describe specific contributions 
of this research. 

6.1 Thesis Summary 

Trace species of aircraft emissions were shown to have potential impacts on the local and global atmosphere, 
including changes in radiative forcing and ozone depletion. A summary of prior work was provided, which has 
concluded that most trace species chemistry occurs early, thus the need exists for better characterization of trace 
species in the intra-engine environment. 

This thesis presented an overview of a modeling methodology developed to simulate the evolution of trace 
species in the post-combustor flow path of gas turbine engines. The modeling methodology incorporates 
various levels of detail to accurately and efficiently predict levels of intra-engine trace species by considering 
key parameters affecting their evolution, specifically temperature, pressure, residence time, and species 
concentration. The model is intended to improve the overall understanding of the fundamental physical 
processes that effect trace species evolution and to serve as a predictive design tool which can direct the 
development of new engine technologies which reduce undesirable aviation emissions. 

Modeling tools and improvements made to these tools were discussed. The existing tools consist of the PRE 
grid generation code, PROCESS grid pre-conditioner code, POST grid post processor code, CNEWT and 
CALCHEM flow chemistry solvers, and CNEWTVS6 flow visualization code. Improvements in modeling 
tools were made, which allow for more complicated geometries, multiple inlet and exits, improved execution 
time for solutions with chemistry, and a specification of a pressure, temperature, and species concentration 
deficits associated with a wake of an upstream blade or blade row. 

Several validation exercises were performed to benchmark these modeling tools. An attempt was made to 
model a VPFR experiment to demonstrate the accuracy of the modeling tools, however, a limitation to the flow 
solver prevented this ultimate objective. The validation exercise did nonetheless provide the opportunity to test 
several code improvements, investigate the limitations of the modeling tools, and analyze flow features of the 
reactor. A chemical mechanism study which investigated the effects of the chemical mechanism on intra- 
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• Time scale analyses are useful in formulating a overall modeling strategy by indicating the critical areas ot 
the engine flow path to investigate with higher fidelity modeling. 

• A time scale analysis is a simple parametric study which incorporates thermodynamic potential, chemical 
kinetics, and residence time. The output is a plot of severity parameter for a given species of interest, 
although, a preliminary estimate of species evolution can be obtained. 

• Time scale analyses can also be used to provide insight into the physical phenpmena influencing trace 
species evolution. A few example time scale analyses were used to show: 

• Blade cooling can enhance sulfur conversion. 

• Engine operating condition can influence sulfur conversion, take-off conditions can have higher sulfur 
conversion than cruise conditions. 

• The OH pathway to S0 3 can be more active than the O pathway. 

• A simple method for specifying chemistry initial conditions involves perturbing an equilibrium 
composition , re-balancing the species based on mass, and doing a stabilizing kinetics calculation. 

• High and low fidelity models involve integrating the equations of mass, momentum, and energy over a 
discrete domain coupled with chemical kinetics to simulate intra-engine chemistry. A high fidelity model 
can incorporate many multi-dimensional effects (i.e. non-uniformities) while a low fidelity model uses 
average quantities. 

• CNEWT and CALCHEM have been useful modeling tools for studying intra-engine trace species flow- 
chemistry problems. CNEWT is a code created by Lukachko [1] which was built upon the NEWT 
turbomachinery CFD code from Dawes [51]. CALCHEM is a simple, ID version of CNEWT. 


• Other supporting codes necessary for turbine chemistry modeling are PRE, PROCESS, POST , and 

CNEWTVS6. 

• The improvements made to existing modeling tools include: 

• A new CAD solid model based grid generator which allows for the simulation of more complicated, 
arbitrary geometries. The output is a truly unstructured grid. The tools were in a developmental stage 
and a few difficulties were encountered with the translation step and numerical issues in the volume 
mesher. The VPFR was successfully meshed with the new grid generator. 

• Addition of multiple inlet/exit capability which allows for simulation of mixing flows. A test case and 
the Princeton VPFR were used to validate the multiple inlet/exit code modifications. 

• Addition of parallel chemistry subroutines. A test case showed a factor of 3.4 improvement in time on 
6 compute nodes. The accuracy of the parallel chemistry modifications was verified. 


139 



• The pre-test simulations suggested testing at the highest combustor exit temperature attainable, as well as, a 
lower power setting. Also, S0 3 and HOHO were selected as important species to monitor. 

• A time scale analysis indicated that most S0 3 production occurs in the HPT l at max power conditions 
(approximately 4%), The analysis indicated that there is several orders of magnitude less sulfur conversion 
after the HPT1, and modeling resources were thus concentrated on the HPTl. 

• PRE was used to generate grids of the NGV and rotor. A velocity triangle analysis was used to derive 
necessary intra-stage fluid parameters for the high fidelity modeling. The species initial condition for each 
condition were derived using the technique outlined in Section 2.3.2. 

• Two types of high fidelity modeling were employed, a mixed-out case and a wake model case. All 
solutions were deemed converged using the criteria from Section 4.4. Also, the flow features such as 
deviation angle, streamlines, and inlet/exit conditions were checked to ensure the quality of each solution. 

• The boundary layers and trailing edge geometry were not fully resolved for the grids used, however, they 
are not expected to have a significant impact on the trace chemistry results. 

• HPTl NGV blade-row flow-through time ranged from 0.18 to 0.26 ms and HPTl rotor blade-row flow- 
through time ranged from 0.10 to 0.15 ms. The flow-through times computed from the low fidelity 
modeling were similar. 

• A summary of all high fidelity modeling results is contained in Table 5.5, Table 5.6, and Table 5.7. 
Examples of some detailed results of the high fidelity modeling can be found in Section 5.4. 2. 3 or 
Appendix C. 

• Low fidelity models were used to model the HPTl and HPTl exit to nozzle exit flow paths. A high and 
low temperature streamline, as well as the averaged profile, through the HPTl were investigated. The 
HPTl low fidelity model can be compared to the HPTl high fidelity model to investigate the impacts of 2D 
versus ID modeling. 

• The HPTl exit to nozzle exit low fidelity model profile was improved using a curve fit to results from the 
HPTl high fidelity modeling and the available data. 

• A summary of all low fidelity modeling results is contained in Table 5.6 and Table 5.7. Examples of some 
detailed results of the low fidelity modeling can be found in Section 5.5. 2. 2 or Appendix C. 
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Appendix A 


Princeton VPFR Validation Modeling 


The Princeton University’s Variable Pressure Flow Reactor (VPFR) is an experimental apparatus specifically 
designed to measure chemical kinetic data at the high temperatures and pressures typical of energy conversion 
systems. The VPFR is ideally suited to evaluate the capability of the modeling tools to represent many flow 
features of interest, mixing, and passively reacting flow. The approach was to model a published flow reactor 
experiment and compare the simulation data to the experimental data in an attempt to benchmark the accuracy 
of the modeling tools. The plan was to first simulate the flow-only, then flow with inert species, and finally 
passively reacting flow. In the end, the objectives of the flow reactor were reduced due to the restriction of the 
modeling tools to simulation of compressible flows since the flow-only solution did not accurately represent the 
actual VPFR experimental data. 

The objectives of the Princeton VPFR validation modeling effort can be found in Section 4.2.1 and a summary 
of results can be found in Section 4.2.2. 


VPFR Modeling Set Up 

Grid Generation 

The flow reactor is essentially a cylindrical tube with an inlet section, a baffle, an injector assembly, a diffuser, 
and a test section. Figure A. 1 shows the basic flow reactor geometry and identifies some important features. 
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generation procedure are detailed in Section 3.3.1 and [66]. The solid model and final flow reactor grid are 
shown in Figure A. 2. 



Figure A.2 VPFR solid model and grid detail near baffle and injector assembly 


As mentioned in Section 3.3.1, the new grid generator had difficulty meshing geometries with large aspect 
ratios. This was overcome in the flow reactor grid by scaling the radial dimension in the solid model by a factor 
of 2 to reduce the aspect ratio then re-scaling the resultant grid back to the original aspect ratio. 

To improve convergence time the grid is typically initialized with an approximate guess for the fluid variables 
(temperature, pressure, density, and velocity). In particular, an axial velocity in the direction of flow was 
estimated and used for the flow reactor. However, due to the complex nature of the flow path in the reactor it 
was found to be better to start with zero initial velocity. In the inlet and test section the flow is nearly axial, 
however, near the injection site it moves radially outward around the baffle and then radially inward into the 
diffuser. In other words, the baffle is perpendicular to the flow which means that the fluid near the baffle had a 
poor initial guess. Therefore, the grid with a prescribed axial velocity would cause the code to enter the “panic 
smoothing” routine which is triggered by an abnormally low pressure since the fluid on the downstream side of 
the baffle was moving away from the baffle causing a low pressure region. In most cases, the code w'ould 
eventually crash once the pressure became very small. Initializing the grid with a zero velocity resolved this 
problem at the expense of having a large transient in the beginning of the solution and long convergence time. 

Experimental Conditions And Modeling Assumptions 

Several aspects of the experimental operating conditions need to be considered relative to the modeling 
assumptions inherent in CNEWT when evaluating the applicability of the code to the flow reactor problem. A 
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results for Mach numbers as low as 0.2 to 0.3 [69] it is in general not well suited to simulate the experimental 
conditions typical of the VPFR operation. Modifying the code by changing the flow variables to double 
precision could extend the compressibility limits of the code by helping to resolve the pressure differences [69], 
however, it is uncertain how much this modification would help. 

In general, the error introduced by exceeding the compressibility limit was large and the fluid solution did not 
represent experimental or empirical features well. Initially, a similarity approach was attempted to improve the 
solution quality by using non-dimensional analysis. The Mach number could be increased (to 0.1 < M < 0.3) 
while maintaining the incompressibility assumption by altering the numeric artificial viscosity. The Reynolds 
number is defined as: 


Re D 


pVD 


(A. 1) 


Matching the Reynolds number of the simulation to the experiment is an important step to reproducing the flow 
features. However, it was found to be practically impossible to scale the reactor with complete similarity [59], 
[60]. Using the available independent variables, it is very difficult to match the non-dimensional numbers 
necessary to equate the various other parameters required to accomplish modeling with even partial similarity 
for this situation [59]. The fluid flow features, mixing, and chemical kinetics are all important to accurately 
represent the flow reactor, some of the difficulties encountered in the similarity approach for this situation will 
be detailed further in “Steady Non-Reacting Flow Solution. 


Other possibilities to overcome the compressibility problem with the flow reactor modeling effort was to use the 
double precision fluid variables, implement a pre-conditioner to the current solution algorithm [70] or adopt a 
entirely new solution algorithm. Using the double precision option was not expected to resolve the problem and 
the other options were not easily implemented and were outside the context of the research agenda. 


Expected Flow Reactor Results 

Diffuser Pressure Recovery Coefficient 

The pressure recovery coefficient is one method to measure the performance of a diffuser. It is defined as: 


c 


p 


= Pe-Pl 
Po'-Pt 


(A. 2) 


There are ample data in the literature relating diffuser design parameters and operating conditions to this 
pressure recovery coefficient. These pressure recovery maps are useful to diffuser designers. In the context of 
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VPFR operating point 

intersection of red lines (L D = 15, AR = 13.7) 
Fully developed two dimensional stall 

Initial separation form computation 

intersection of blue lines (LVD - 8.5, AR - 2.5) 
Near line of appreciable stall 


Figure A.4 Diffuser flow regime map [39] 


Velocity Profiles 

The most informative method to evaluate the quality of the computational solution was to compare the velocity 
profiles to experimental data and empirical correlations. The VPFR was designed to operate such that the 
velocity profiles transition quickly to fully-developed turbulent velocity profiles [26]. These profiles have a 
characteristic “top hat” shape (see Figure A. 5 and Figure A. 6). 

Empirical Correlation 

The Power Law is one empirical correlation commonly used to predict velocity profiles in pipe flow r : 


u 

V~ 


1 — 


R 


(A. 3) 


Where V c is the centerline velocity, u is the time average x component of velocity, and R is the pipe radius. For 
the flow r reactor case, n, which is a function of Reynolds number, is about 6. Figure A. 5 shows the expected 
velocity profile for the VPFR derived from the Power Law. 
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Figure A. 5 Typical velocity profiles from Power Law [38] 


Experimental Data 


For comparison, some velocity profile data specific to the Princeton VPFR is given in Figure A. 6. 



Another point of comparison are other experimentally measured velocity profiles for conical diffusers, such as 
that given in Figure A. 7. The red labels in the figure indicate the profiles for diffusers with design parameters 
similar to those of the Princeton VPFR. This figure shows that the profiles become more asymmetric and 
peaked with higher AR and L/D. 
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Note: Both plots for 20 = 8°, reactor 20 = 1 0 a 
Flow reactor entrance l 0 /D 0 = 3 



Note: Both plots for AR = 4, reactor AR = 13.7 
/ 



Figure A.7 General velocity profiles in diffusers [40] 


Steady Non-Reacting Flow Solutions 

Several modeling runs were made in an attempt to simulate the flow reactor and characterize the capability and 
limitations of the CNEWT code. Two cases were selected for presentation in this thesis to illustrate the major 
findings of the flow reactor validation effort. The first case, labeled “Run 1 ” has the lowest Reynolds number 
in the test section for which a converged solution could be obtained. The second case, labeled “Run 12/ has a 
higher Reynolds number. Further differences between the cases and details of the solutions follow in “Low 
Reynolds Number Solution (Run 1)” and “High Reynolds Number Solution (Run 12)”. 

Genera! Results 


A summary of results from the two selected cases is given in Table A. 1. The Reynolds number for both 
solutions is several times larger than the experimental conditions of the VPRF. 
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Run 1 

Run 12 

AP (kPa) 

1 

30 

V* max (m/s) 

-85 

-452 

V x min (m/s) 

11 

163 

Mmax 

0.25 

1.7 

p, 

93,004 

55,181 

Pot 

96,677 

122,028 

Pe 

96,994 

106,128 

Cp 

1.09 

0.76 

V x average exit (m/s) 

ESHI 

-16 

Reo 

31,400 

116,700 

Mass flow in (kg/s) 

2.7e-3 

8.4e-3 

Mass flow out (kg/s) 

2.5e-3 

1.1 e-2 

Mass flow error 

7.2% 

Injectors ?? 

Shock position (cm) 

none 

5 


Table A.l Summary of selected results for VPFR validation runs 


Figure A. 8 shows a cross-section of the flow reactor with locations of nine axial stations from 15 to 95 cm. In 
the simulations presented here, the flow travels in the negative X direction (right to left in the figure). 



Figure A.8 Axial position X (cm) for flow reactor modeling 
Low Reynolds Number Solution (Run 1) 

Run 1, a low Re solution, was solved on a preliminary coarse grid which had about 60,000 cells and no injector 
flow. It was started with a AP across the reactor main inlet to exit of 30 kPa which yielded a Mach number of 
1.23 at the throat of the diffuser and a Reynolds number of 82,400 in the test section. The convergence history, 
which was typical of most flow reactor runs, is plotted in Figure A. 9. 
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Figure A. 9 Flow reactor Run 1 convergence history 

The peak suction (lowest pressure in the domain) converged to about 87 kPa, log(RMS) (the root mean square 
residual) was below 1.0E-8, and the continuity (mass flow error) oscillated between 10% and 60%. The mass 
flow error and oscillations were larger for lower Mach number flows, and Run 1 represents the edge of solution 
stability. The large mass flow error was eventually attributed to the codes inability to compute incompressible 
flows. 

As mentioned previously, a similarity approach was planned as a potential method to overcome the 
compressibility issue and continue the flow reactor validation efforts. In an attempt to match the experimental 
conditions for the Princeton VPFR, the inlet pressure was sequentially reduced until the Mach number at the 
throat of the diffuser was less than 0.3, which occurred at a pressure difference of 1 kPa. Contours of axial 
velocity are given in Figure A. 10. The average Mach number was 0.013 in the test section. 



Figure A. 10 Axial velocity from -85 to 11 m/s (note: contours clipped for presentation) 
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The diffuser and test section were divided into nine axial stations (see Figure A.8) and velocity profiles were 
obtained at each. The profiles are given in Figure A.l 1 for Run 1. 



Figure A. 11 Velocity profiles at nine axial stations 

The velocity profiles for Run I are not similar to the “top hat” profiles typical of fully developed turbulent flow. 
They are peaked and the curvature is opposite to that which is expected. The results from Run 1 exhibit no 
recirculation, however the infective velocity profiles indicate the possibility of separated jet flow. Furthermore, 
the pressure coefficient was unrealistically greater than unity. Considering the very low Mach number in the 
bulk of the test section (M ~ 0.01 « 0.3), it was suspected that the majority of the error encountered was 
associated with using the code beyond the compressibility limit. However, potential discrepancies caused by 
the mismatched Reynolds number or boundary layer grid resolution were still investigated. 

Carrying on with the plan to match the Reynolds number by increasing the artificial viscosity, three additional 
solutions were started from Run 1. The laminar and turbulent viscosity were increased by a factor of 10, 20, 
and 50 for these solutions. The results of this study are summarized in the following table: 
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Run 1 

Run 2 

Run 3 

Run 4 

Description 

lx p 

lOx p 

20x p 

50x p 

V* max (m/s) 

-85 

-73 

-72 

-61 

V x min (m/s) 

11 

10 

9 

7 

Mmax 

0.247 

0.210 

0.206 

0.176 

p, 

93,004 

94,444 

95,314 

96,401 

Pot 

96,677 

96,898 

97,018 

97,211 

Pe 

96,994 

96,999 

96,999 

96,999 

c D 

1.09 

1.04 

0.99 

0.74 

V x average exit (m/s) 

-4.7 

-3.8 

-3.2 

-2.3 

Reo 

31,400 

2,590 

1,080 

300 

Mass flow in (kg/s) 

2.7E-3 

2.3E-3 

1.9E-3 

1.3E-3 

Mass flow out (kg/s) 

2.5E-3 

1.8E-3 

1.5E-3 

1.5E-3 

Mass flow error 

7.2 % 

25.1 % 

29.8 % 

12.7% 


Table A.2 Summary of Reynolds number matching study 


The Reynolds number for the lOx case was already beyond the goal of reaching 3,000 to 15,000. The profiles 
for that case were very similar to Run 1 and are not repeated here. Instead the velocity profiles for Run 4, the 
extreme case with 50x viscosity, are shown in Figure A. 12. As the numerical viscosity increased, the velocity 
profiles became more inflected and still remained peaked, thus matching the Reynolds number did not move 
toward a more realistic solution. In fact, increasing the viscosity lowered the flow Mach number pushing the 
flow more into the incompressible regime which would make the error associated with the incompressibility 
higher. 


£ 

£ 


Velocity Profiles; 5Qx vise 
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Figure A.12 Velocity profies for 50x viscosity case (Run 4) 
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Also, it was initially thought that the wall functions or grid boundary layer resolution could be causing the poor 
solution quality. The grid post processor was used to refine the grid near the outer walls of the diffuser and test 
section. The first node from the outer wall for the coarse grid was 5.9 mm away. Several levels of refinement 
were used until a grid with nearly 670,000 cells was created with the first node located 0.6 mm from the outer 
wall, see Figure A. 13. The results of this study are summarized in the following table: 




Run 5 

Run 6 

Grid 

coarse 

2 level refinement 

3 level refinement 

Cells 


205,682 

668,681 

V x max (m/s) 

-85 

-71 

-85 

V x min (m/s) 

11 

10 

11 

Mmax 


0.209 

0.245 

p. 


94,438 

93,338 

Po. 



96,543 

Pe 

96,994 

97,000 

97,000 

c p 

1.086 

1.006 

1.143 

Vx average exit (m/s) 

-4.7 

-3.9 

-4.3 

ReD 

31,400 

26,300 

28,600 

Mass flow in (kg/s) 

2.7E-3 

2.3E-3 

2.5E-3 

Mass flow out (kg/s) 

2.5E-3 

2.2E-3 

2.6E-3 

Mass flow error 

7.2 % 

5.9 % 

3.4 % 


Table A.3 Summary of boundary layer grid refinement study 


The velocity profiles from Run 6 with the highest degree of boundary layer refinement are shown in Figure 
A. 14. Again, the velocity profiles did not move in the direction of resembling the experimentally or empirically 
derived profiles. The velocity profiles were still peaked/in flee ted and the pressure recovery coefficient was 
unreasonable, lending further evidence that the main source of error in these flow reactor simulations is due to 
the compressibility limits of the code. 
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| Coarse mesh: 60,166 cells, about 5.9 mm to first node 



1 Mesh with two levels of refinement: 205,682 cells, about 1 .5 mm to first node 



Mesh with three levels of refinement: 668,681 cells, about 0.6 mm to first node 


Figure A.13 Mesh refinement in flow reactor test section 
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Velocity Profiles 1 x vise. 2 level refinment in test section & 1 level in diffuser 
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Radial Distance (m) 


Figure A. 14 Velocity profiles for case with refined boundary layer grid (Run 6) 

High Reynolds Number Solution (Run 12) 

The second case to be discussed in some detail, Run 12, was solved on a finer grid which had about 270,000 
cells (see Figure A. 2). It was maintained at a higher AP of 30 kPa from the main inlet to exit and had flow at 
the two injector inlets. The injector inlets were set at a pressure ratio of about 9 to 1 relative to the pressure at 
the main exit. In this case the mass flow error was not computed correctly because the injector mass flow could 
not be easily accounted. Figure A. 15 shows the axial velocity contours plotted on streamlines for a solution 
dumped at 10,000 and 35,000 iterations. 
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1 0,000 iterations: 




Figure A. 15 Axial velocity from -452 to 163 m/s at 10,000 (top) and 35,000 (bottom) iterations (note: 

contours clipped for presentation) 


The average Mach number in the test section was 0.05. This solution shows that the flow separates in the 
diffuser and recirculation zones begin to form. The point of separation is roughly 23 cm downstream of the 
injection site. This point corresponds to a L/D of about 8.5 and a AR of about 2.5 and was anticipated based on 
the flow regime map presented in Figure A.4. 
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Figure A. 16 Axial velocity profiles at nine axial stations for 10,000 and 35,000 iteration solutions 
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The velocity profiles were peaked prior to the recirculation, followed by flat profiles characteristic of the 
experimental data. The flat profiles after the stall cells are likely due to the large scale mixing, exchange of low 
momentum fluid at the wall with high momentum fluid in the core. The profiles show a deficit near the 
centerline which is typical of flow around a bluff body (the injector). These observations would suggest that the 
diffuser is operating with large scale unsteady stall similar to the expected operation based on the flow regime 
map. The flat experimental velocity profiles of the Princeton VPFR could be due to this averaged unsteadiness 
rather than high turbulence at the diffuser inlet. However, the code was being run in steady mode with a non- 
uniform time step, therefore, further investigation into the unsteadiness was required (see Unsteady Non- 
Reacting Flow Solution). 

In the higher AP solution. Run 12, the diffuser was choked and there was a shock positioned 5 cm into the 
diffuser. Figure A. 17 shows the supersonic region where the flow accelerates from the diffuser throat to the 
shock where the maximum Mach number was 1.7. 



Figure A. 17 Mach number contours near diffuser throat showing subsonic to supersonic transition 

The pressure recovery coefficient for this case was 0.76. However, this lower value was due to the losses 
associated with the shock rather than a better representation of the diffuser operation at the desired flow 
conditions. 

The presence of the shock was an additional complication when considering the plan for a similarity solution. It 
was apparent from Run 1 that having a Mach number of about 0.2 to 0.3 solely in the throat of the diffuser was 
not sufficient to overcome the compressibility limits of the code. The bulk of the flow field would have to be 



above a Mach number of 0.2 to obtain an accurate solution. The flow reactor geometry is such that the throat of 
the diffuser chokes well before the Mach number in the test section increases to a reasonable level for the flow 
solver. From a compressible flow table for an ideal gas (with y = 1.4), assuming isentropic flow with an area 
ratio of 13.7, the Mach number at the diffuser exit would only reach about 0.04 when the throat chokes [38]. Or 
in other words, to use geometric scaling and assuming a desired Mach number in the test section of less than 
0.25, the maximum diffuser area ratio would be 2.4. With this Mach number and area ratio, the velocity in the 
test section is much higher and the length of the test section must be increased to capture the same flow reactor 
residence time. The test section becomes prohibitively long (about 25 m), making lab measurements difficult 
and the geometry harder to mesh (worse aspect ratio). 

Figure A. 18 shows a detail of the Z-direction velocity contours overlaid on streamlines near the injectors. This 
figure gives qualitative evidence that the multiple inlet/exit code modifications are functioning properly. At this 
point, no attempt has been made to quantitatively evaluate the ability of the code to simulate the mixing process 
associated with the injected fluid. 



Figure A. 18 Z direction velocity component near injectors 


Matching the momentum ratio of the injection jet and the opposing primary flow is another complex problem 
with regard to performing a similarity solution. 
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Unsteady Non-Reacting Flow Solution 

After discovering the flow separation in Run 12 and knowing that the diffuser may operate in a unsteady 
fashion, it was decided to continue Run 12 using a uniform time step to get a proper unsteady solution. The 
choice of time step was constrained by the size of the smallest cell. Choosing too large a time step would cause 
the code to crash because the flow could not be resolved on the scale of that smallest cell. The fine mesh near 
the injector limited the time step to 0.2455E-06 seconds per iteration. 

This unsteady case gives more insight into the operating characteristics of the flow reactor and the capability of 
the code to simulate unsteady flow. Figure A. 19 shows contours of axial velocity overlaid on streamlines at 
60,000 iterations, just prior to the time when the large initial transient leaves the domain. 



Figure A. 19 Axial velocity just before 60,000 iterations 

Figure A.20 shows the instantaneous velocity profiles for the unsteady Run 12 at 100,000 iterations. The shape 
of the velocity profiles are similar to the steady case and are no more representative of the “top hat” profiles 
expected of fully developed turbulent flow. The deficit near the centerline from the bluff body and inflected 
profile shape are still evident. In fact, the negative velocity near the outer radius is indicative of reversed flow 
and signifies separation. 
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Velocity Profiles run 12 with 100 000 iterations 



Radial Distance (m) 


Figure A.20 Velocity profiles for unsteady Run 12 at 100,000 iterations (not time averaged) 

The unsteady Run 12 solution was continued for one cycle after the initial transient left the domain. The 
limitations on time step meant that nearly 30 days of computational time were required to compute that one 
cycle. A recirculation vortex was observed to form near the outer wall at the point of first separation then travel 
downstream while growing in size until it is shed into the test section, see Figure A.21. The periodic vortex 
shedding was estimated to occur at a frequency of 151 Hz by computing the time elapsed between one cycle in 
the unsteady computation. 



Figure A.21 Streamlines with axial velocity contours for unsteady Run 12 showing one cycle of the 

periodic vortex shedding 

An animation of the initial transient leaving the domain and one cycle of the periodic vortex shedding can be 
found in Figure F.9 and Figure F. 10, respectively (see Appendix F). 
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Acoustic Resonance Calculation 


The VPFR was suspected to have an acoustic resonance at certain operating conditions [26]. It was 
hypothesized that these acoustic pressure waves could be linked to the vortex shedding frequency. Using the 
wave equation, the resonant frequency for the modeled portion of the flow reactor was estimated as: 

/ = - = — = ( 34 ° Hfo - L = 1 76 Hz (A. 4) 

/l 2 L 2- (0.964 m) 

This theoretical frequency compared well to the value of 151 Hz predicted in the unsteady simulation. The 
simulation did not include the entire length of test section of the real flow reactor. Using L = 1.7 m for the 
actual length from the baffle to the exit of the flow reactor, the vortex shedding frequency for the experimental 
facility is predicted to be about 100 Hz. 

VPFR Modeling Summary 

Limitations of the flow solver hampered the ultimate objective of validating the modeling tools through a 
successful reacting flow simulation of the Princeton VPFR. Although the validation exercise was not 
completed, many code improvements were tested, a considerable amount was learned about the limitations of 
the tools, and some potential issues with the flow reactor design were identified. 

Refer to Section 4.2.2 for a more detailed summary of the results from the Princeton VPFR validation modeling 
exercise. 
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Appendix B 


Turbine Chemistry Mechanism 


Several turbine chemistry chemical mechanisms were evaluated for use in the calculations. The mechanism 
referred to as the “Mueller et al. (2000)-truncated mechanism” was selected for the NASA/DERA engine 
simulations. The mechanism has 29 species and 73 reactions. For more information regarding the selection of 
chemical mechanism see Section 4.3. This appendix lists the output from CHEMKIN which describes the 
mechanism. 


CHEMKIN INTERPRETER OUTPUT: CHEMKIN-II Version 3.9 Aug. 1994 

DOUBLE PRECISION 


ELEMENTS 

CONSIDERED 

ATOMIC 

WEIGHT 

1. C 

12.0112 

2. H 

1.00797 

3. 0 

15.9994 

4. S 

32.0640 

5. N 

14.0067 

6. AR 

39.9480 


C 

P H 
H A 
A R 


SPECIES 

CONSIDERED 

S 

E 

G 

E 

MOLECULAR 

WEIGHT 

TEMPERATURE 
LOW HIGH 

ELEMENT COUNT 
C H O S N 

AR 

1 . 

CO 

G 

0 

28.01055 

300 

5000 

1 

0 

1 

0 

0 

0 

2. 

C02 

G 

0 

44.00995 

300 

5000 

1 

0 

2 

0 

0 

0 

3. 

02 

G 

0 

31.99880 

300 

5000 

0 

0 

2 

0 

0 

0 

4. 

H20 

G 

0 

18.01534 

300 

5000 

0 

2 

1 

0 

0 

0 

5. 

H 

G 

0 

1.00797 

300 

5000 

0 

1 

0 

0 

0 

0 

6. 

O 

G 

0 

15.99940 

300 

5000 

0 

0 

1 

0 

0 

0 

7 . 

OH 

G 

0 

17.00737 

300 

5000 

0 

1 

1 

0 

0 

0 

8. 

H02 

G 

0 

33.00677 

200 

3500 

0 

1 

2 

0 

0 

0 

9. 

H202 

G 

0 

34.01474 

300 

5000 

0 

2 

2 

0 

0 

0 

10. 

NO 

G 

0 

30.00610 

300 

5000 

0 

0 

1 

0 

1 

0 

11. 

NO 2 

G 

0 

46.00550 

300 

5000 

0 

0 

2 

0 

1 

0 

12. 

NO 3 

G 

0 

62.00490 

300 

5000 

0 

0 

3 

0 

1 

0 

13. 

HNO 

G 

0 

31.01407 

300 

5000 

0 

1 

1 

0 

1 

0 

14. 

HONO 

G 

0 

47.01347 

300 

5000 

0 

1 

2 

0 

1 

0 

15. 

HNO 3 

G 

0 

63.01287 

300 

5000 

0 

1 

3 

0 

1 

0 

16. 

S02 

G 

0 

64.06280 

300 

5000 

0 

0 

2 

1 

0 

0 

17. 

SO 3 

G 

0 

80.06220 

300 

5000 

0 

0 

3 

1 

0 

0 

18. 

HOSO 

G 

0 

65.07077 

300 

5000 

0 

1 

2 

1 

0 

0 
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19. 

H0S02 

G 

0 

81.07017 

300 

5000 

20. 

HS02 

G 

0 

65.07077 

300 

5000 

21. 

SO 

G 

0 

48.06340 

300 

5000 

22. 

HCO 

G 

0 

29.01852 

300 

5000 

23. 

H2 

G 

0 

2.01594 

300 

5000 

24. 

N2 

G 

0 

28.01340 

300 

5000 

25. 

AR 

G 

0 

39.94800 

300 

5000 

26. 

C(S) 

S 

0 

12.01115 

300 

5000 

27. 

N 

G 

0 

14.00670 

200 

6000 

28. 

H2S04 

G 

0 

98.07754 

300 

5000 

29. 

S 

G 

0 

32.06400 

300 

5000 


0 13 10 0 

0 12 10 0 

0 0 110 0 

1110 0 0 

0 2 0 0 0 0 

0 0 0 0 2 0 

0 0 0 0 0 1 

1 0 0 0 0 0 

0 0 0 0 1 0 

0 2 4 1 0 0 

0 0 0 1 0 0 


REACTIONS CONSIDERED 


(k = A T**b exp( -E/RT) ) 
A b E 


1. 

H2+M=H+H+M 



4.57E-K19 

1.4 

104400.0 


H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1 . 200E+01 





CO 

Enhanced by 

1 . 900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7.500E-01 




2. 

0+H2-H+OH 



5 . 08E+04 

2.7 

6290.0 

3. 

0+0+M=02+M 



6.16E+15 

0.5 

0.0 


H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1 .200E+01 





CO 

Enhanced by 

1.900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7 . 500E-01 




4 . 

H+02=0+0H 



1.91E+14 

0.0 

16440.0 

5. 

H+02 { +M } =H02 ( +M) 



1 .48E+12 

0.6 

0.0 


Low pressure limit: 

: 0 . 34820E+17 

-0.4 11 OOE+OO 

-0 . 11 1 50E+04 



TROE centering: 

0 . 50000E + 00 

0 . 10000E-29 

0 . 10000E + 31 



H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1.200E+01 





CO 

Enhanced by 

1 . 900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7.500E-01 




6. 

H+0+M=0H+M 



4.71E+18 

1.0 

0.0 


H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1 .200E+01 





CO 

Enhanced by 

1 . 900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7.500E-01 




7. 

0H+H2=H20+H 



2 . 16E + 08 

1.5 

3430.0 

8. 

H20+0=0H+0H 



2.97E+06 

2.0 

13400.0 

9. 

H202 ( +M) =OH+OH { +M) 



2.95E+14 

0.0 

48430.0 


Low pressure limit: 

0 . 12000E+18 

0. 00000E+00 

0 . 4 5500E+05 



TROE centering: 

0 .50000E+00 

0 . 10000E- 89 

0 . 10000E+91 



H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1 . 200E+01 





CO 

Enhanced by 

1 .900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7 . 500E-01 




10. 

OH+H+M=H20+M 



2.21E+22 

2.0 

0.0 


H2 

Enhanced by 

2 . 500E+00 





H20 

Enhanced by 

1 . 200E+01 





CO 

Enhanced by 

1 . 900E+00 





C02 

Enhanced by 

3 . 800E+00 





AR 

Enhanced by 

7.500E-01 




11. 

H02+0=02+0H 



3 . 2 5E + 1 3 

0.0 

0.0 

12. 

H02+H=H2+02 



1 .66E+13 

0.0 

823.0 

13. 

H02+H=0H+0H 



7 . 08E+1 3 

0.0 

295.0 

14. 

H02+OH=H20+02 



2.89E+13 

0.0 

-497.0 

15. 

H02+H02=H202+02 



4 . 2 0E+ 14 

0.0 

11982.0 


Declared duplicate 

reaction. . . 


- 
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16. H02+H02=H202+02 



1.30E+11 

o 

o 

-1629.0 

Declared duplicate 

reaction . . . 





17. H202+0=0H+H02 



9 . 55E+06 

2.0 

3970.0 

18. H202+H=H20+-0H 



2.41E+13 

0.0 

3970.0 

19. H202+H=H02+H2 



4.82E+13 

0.0 

7950.0 

20. H202+0H=H20+H02 



1 . 00E+12 

0.0 

0.0 

Declared duplicate 

reaction . . . 





21. H202+0H=H20+H02 



5 . 80E + 14 

0.0 

9557.0 

Declared duplicate 

reaction . . . 





22. HNO+H=NO+H2 



4 . 4 0E + 1 1 

0.7 

650.0 

23. NO+O ( +M ) =N02 ( +M) 



1.30E+15 

-0.8 

0.0 

Low pressure limit: 

0 . 47200E+25 

-0 . 28700E+01 

0 . 15510E+ 

04 


TROE centering: 

0 . 95700E+00 

0 . 10000E- 89 

0. 83320E+ 

04 


A R 

Enhanced by 

7.500E-01 




24. NO+H ( +M) =HNO { -t-M ) 



1 . 52E+15 

-0.4 

0.0 

Low pressure limit: 

0 .31000E + 20 

-0 . 13200E+01 

0 . 73 520E + 

03 


TROE centering: 

0 . 82000E+00 

0 . 10000E-89 

0 . 10000E+ 

91 


AR 

Enhanced by 

7 . 500E-01 




25. NO+OH { +M) =HONO ( +M ) 



1 . 99E+12 

-0.1 

-721 .0 

Low pressure limit: 

0 . 50800E+24 

-0 . 25100E+01 

- 0 . 67600E+ 

02 


TROE centering: 

0 . 62000E+00 

0 . 10000E- 89 

0 . 10000E + 91 


AR 

Enhanced by 

7 . 500E-01 




26. N02+H2-H0N0+H 



7.33E+11 

0.0 

28810.0 

27. N02+0=02+N0 



1 .05E+14 

-0.5 

0.0 

28. N02+0 ( +M) =N03 { -»-M) 



1 .33E+13 

0.0 

0.0 

Low pressure limit: 

0 . 14900E+29 

-0.40800E+01 

0 . 24670E+04 


TROE centering: 

0 . 82600E+00 

0 . 10000E- 89 

0 . 31910E+04 


AR 

Enhanced by 

7.500E-01 




29. N02+H=N0+0H 



1 .32E+14 

0.0 

362.0 

30. N02+0H(+M) =HN03 (+M) 



2.41E+13 

0.0 

0.0 

Low pressure limit: 

0 . 64200E+33 

-0. 54900E+01 

0 . 23500E + 04 


TROE centering: 

0 . 83700E+00 

0 . 10000E- 89 

0 . 16570E+04 


AR 

Enhanced by 

7 . 500E-01 




31. H02+N0=N02 +OH 



2.11E+12 

0.0 

-479.0 

32. N02+N02=N03 +NO 



9.64E+09 

0.7 

20920.0 

33. N02+N02=2N0+02 



1.63E+12 

0.0 

26120.0 

34. HN0+O=OH+NO 



1 . 81E + 13 

0.0 

0.0 

35. HN0+0H=H20+N0 



1 . 30E+07 

1.9 

-956.0 

36. HN0+N02 -HONO+NO 



6 . 02E + 11 

0.0 

1987.0 

37. H0NO+0=0H+N02 



1.20E+13 

0.0 

5961.0 

38. H0N0+0H=H20+N02 



1 .70E + 12 

0.0 

-520.0 

39. HCO+M=H+CO+M 



1.86E+17 

-1.0 

17000.0 

H2 

Enhanced by 

2 . 500E+00 




H20 

Enhanced by 

1 . 200E+01 




CO 

Enhanced by 

1.900E+00 




C02 

Enhanced by 

3 . 800E+00 




40. HC0+02=C0+H02 



7 . 58E+12 

0.0 

410.0 

41. HC0+0=C0+0H 



3.02E+13 

0.0 

0.0 

42. HCO+H=CO+H2 



7.23E+13 

0.0 

0.0 

43. HCO+OH=CO+H20 



3.02E+13 

0.0 

0.0 

44. CO+O ( + M) =C02 { +M} 



1.80E+10 

0.0 

2384.0 

Low pressure limit: 

0 . 13500E+25 

- 0 . 27 880E+01 

0.41910E+04 


H2 

Enhanced by 

2 . 500E+00 




H20 

Enhanced by 

1 .200E+01 




CO 

Enhanced by 

1 .900E+00 




C02 

Enhanced by 

3 . 800E+00 




45. C0+02=C02+0 



2 . 53E+12 

0.0 

47700.0 

46. C0+0H=C02 +H 



1.40E+05 

1.9 

-1347.0 

47. C0+H02=C02+0H 



3 . 01E+13 

0.0 

23000.0 

48. NO+HCO=HNO+CO 



7.23E+12 

0.0 

0.0 

49. N02+HC0=C0+H0N0 



1.26E+23 

-3.3 

2354.0 

50. N02+HC0==H+C02+N0 



8 . 43E+1 5 

-0.8 

1927.0 

51. N02+C0=C02+N0 



9.03E+13 

0.0 

33780.0 

52. S03+0=S02+02 



4 .40E+11 

0.0 

6100.0 

53. S03 + SO= S02 + S02 



1.00E+12 

0.0 

4000.0 

54. S02+0 { +M ) =S03 ( +M) 



9.20E+10 

0.0 

2384 . 0 

Low pressure limit: 

0 . 400Q0E+29 

- 0 . 40000E+01 

0 . 52500E+04 


N2 

Enhanced by 

1 . 300E+00 
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H20 

Enhanced by 

1.000E+01 




55. 

S02+0H ( +M } =H0S02 (+M) 



1 .21E+12 

0.0 

0.0 


Low pressure limit: 

0 . 1 8700E+32 

-0 . 46100E+01 

0 . 20500E+04 



TROE centering: 

0 . 3 5000E+00 

0 . 1Q000E-29 

0 . 10000E+31 



H20 

Enhanced by 

1.000E+01 




56. 

SO2+0H=SO3+H 



4 . 90E+01 

2.7 

23800.0 

57. 

SO+O ( +M) =S02 ( +M) 



3 . 20E+13 

0.0 

0.0 


Low pressure limit: 

0 . 29000E+25 

-0 . 29000E+01 

Q.0000QE+00 



TROE centering: 

0 . 55000E+00 

0 . 10000E-29 

0 . 10000E+ 

31 



N2 

Enhanced by 

1 .500E+00 





H20 

Enhanced by 

1 .000E+01 




58. 

S0+0H=S02+H 



5 . 20E+13 

0.0 

0.0 

59. 

SO+OH+M=HOSO+M 



8 . 00E+21 

-2.2 

830.0 

60 . 

S0+02=S02+0 



6 . 20E+03 

2.4 

3050.0 

61. 

H0S0+M=S02+H+M 



5 . 90E + 3 4 

-5.7 

50900.0 

62. 

H0S0+0H=S02+H20 



1. 002+12 

0.0 

0.0 

63. 

HC5S0+02 = S02 + H02 



1 .OOE + 12 

0.0 

1000.0 

64. 

HS02+M=S02+H+M 



1 . 20E+28 

-4.1 

18900.0 

65. 

HS02+M=H0S0+M 



1.10E+21 

-2.0 

29900.0 

66. 

H0S02=H0S0+0 



5.40E+18 

-2.3 

106300.0 

67. 

H0S02+M=S03+H+M 



3 . 20E+16 

-0.8 

53700.0 

68. 

H0S02+H=S02+H20 



1. OOE + 12 

0.0 

0.0 

69. 

H0S02+0=S03+0H 



5. OOE+12 

0.0 

0.0 

70. 

HOS02 + 0H=SO3 + H20 



1. OOE+12 

0.0 

0.0 

71. 

H0S02+02 = S03 +H02 



7.80E+11 

0.0 

656.0 

72. 

S02+N02=S03+N0 



6.30E+12 

0.0 

27000.0 

73 . 

S03 +H20=H2S04 



7.23E+08 

0.0 

0.0 


NOTE: A units mole- cm- sec -K, E units cal/mole 


NO ERRORS FOUND ON INPUT ... CHEMKIN LINKING FILE WRITTEN. 

WORKING SPACE REQUIREMENTS ARE 
INTEGER: 1569 

REAL: 1675 

CHARACTER: 35 
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Appendix C 


HPT1 Velocity Triangle Analysis 


Due to the limited data available for the NASA/DERA engine, a velocity triangle analysis was used to obtain 
the intra-stage data needed for the HPT1 high fidelity modeling. This type of analysis is common in 
turbomachinery textbooks, and further details can be found in [7], [36], or [37]. This appendix presents the 
specific procedure along with some example results for the velocity triangle analysis referred to in Section 0: 


The axial stations of interest were labeled as follows: 

Station A: Combustor exit (8) (numbers refer to DERA data designation) 
Station B: HPT1 rotor inlet (9) 

Station C: HPT1 rotor exit 

Station “Overall HPT exit”: HPT2 rotor exit (11) 


The subscripts a, b, and c are used to designate these stations in this analysis, as well as, Figure C. and Table C.. 
The annulus area was calculated at each station. The properties of air R-288.7 J/kgK and y= 1 .34 were assumed 
to be constant over the temperature range of interest, 800-l,300K. 

The fluid density was calculated using the ideal gas law and the relationship between total and static 
pressure/temperature by using an initial guess for the Mach number as follows: 


P a = 



T.= 



Pa = 


RT a 


(C.1) 


(C.2) 


(C.3) 
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Then, the axial velocity and Mach number were calculated as follows: 


U 


Xtt 


m a 

PaK 


(C.4) 


M 


a 


4r RT . 


(C.5) 


An iteration was performed using equations (C.l) to (C.5) until the Mach number was consistent. The 
combustor exit flow was assumed to enter the NGV with zero incidence angle, thus: 


V a =u xa (C.6) 

A similar iterative procedure was used to compute u* b . However, since the total temperature and pressure were 
not explicitly given within the stages, negligible total pressure loss was assumed to occur through the HPT1 
NGV and thus, the “Overall HPT exit” values (which were similar to the combustor exit) were used for this 
intra-stage station. Also, the Mach number was computed using the absolute velocity, V b . The flow was 
initially assumed to have zero deviation, thus the blade metal angle at the trailing edge, (5 b , and u xb were used to 
compute V b as follows: 




\xb 


cos 


(ft) 


(C.7) 


M» = 


Jr-R-n 


(C.8) 


Vector addition was used to subtract the tangential velocity of the rotor, car, from V b to get the velocity, V b ’, and 
inlet angle, p b ’ , relative to the rotor reference frame as follows: 


v * = 


(C.9) 


= v b -cor 

(C.10) 

v= 

4 u lb +v 2 

(C.l 1) 


/ A 


5 * 

if 

> 1 : 

T 

P 

a 

(C.12) 
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Again, since the HPT1 rotor exit condition was not given explicitly the turbine was assumed to have a 50% 
reaction, thus the work per stage was assumed to be equal. A power balance using Station A and the "Overall 
HPT exit” values was done assuming a polytropic efficiency of 0.90 to obtain the total temperature and pressure 
at station C, T Tc and P Tc respectively. An iteration was performed to obtain T Tc and the turbine temperature and 
pressure ratios were computed to obtain P Tc as follows: 


C P ( r T, -Tj = m z c p {T T2 -T Tc ) => T Tc 


T 


1 1 



(C. 13) 
(C. 14) 
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/l 



(C. 15) 


P Tc — 7l JX * P Ta (C.16) 

Finally, an iterative proceedure over equations (C.5) to (C.12) was used to compute the velocity vectors at the 
HPT1 rotor exit. Again a zero deviation was assumed (p c ‘ = rotor blade metal angle at the tri ling edge) and the 
Mach number was computed using the absolute velocity, V c \ as follows: 


y 9 = 

cos {fi f ’) 


v r = v ( . —CD • r 


V , ”V“i + V r ! 


(C. 17) 


(C.18) 

(C.19) 

(C.20) 


f .. \ 


P c = tan 


-1 




M . = 


u r ,. 


r Jr-R-T, 


(C.21) 


(C.22) 


The results of this analysis were reasonable since the fluid flow angle at the exit of the stage is nearly axial as is 
typical in the design of a turbine stage, P c = 3.7°. However, the incidence angle to the rotor, although 
reasonable, was fairly large p b ’ = 37.1°. Therefore, a deviation of the flow from the blade metal angles w'as 
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assumed. The deviation was initially set at 2° for both blades. This lowered the rotor incidence angle to -2.2° 
and changed the rotor exit angle to 28.7° As an example, the resultant velocity triangles and underlying data 
from this analysis for the mamimum power condition with 2° assumed deviation are given in Figure C. and 
Table C.. 
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NASA/DERA Engine Test 
HPT1 Velocity Triangles 
Max Power Condition 
2° of Deviation 


Station A 


— U x a 


HPT1 NGV 




Station B 



Station C 


Dashed lines and X' quantities represent the relative frame 


HPT1 Rotor 


Figure C.l Velocity triangle diagram for maximum power condition with 2 deviation 
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Station overall HPT exit 
HPT2 rotor 0x11(11) 


NASA/P6R A Engine To it. HPT1 Velocity Trianolei. Max Power Condition. Aeaumed 2 degreea of Devietjon 


Station A 

Combuotor exit condition (8) 


R * 

270.38 mm 

N blades = 

80 

t = 

0 mm 

DO a 

22 88 in 

Di = 

19.1 m 

delta D - 

3 76 in 

A a 

123 9119541 m*2 
0 079943036 rrv^ 


Ma = 

<M 920406W ; 

Y = 

1.34 

Pta = 

1 834 1 50 0 Pa 

Tta = 

1,365 2 K 

Pa = 

1,594,021 3 Pa 

Ta * 

1,358 8 K 

R = 

288 7 JAq-K 



4 070 kfllm*3 


m dot = 

45 479 kQ/s 

uxa = 

136.77 nV» 


Station 0 

HT1 rotor inlet (10) 


Nblades = 

100 


t * 

0 mm 

0 

Do = 

22 7 .n 

il! 

Di = 

19 74 in 

^BEES 

delta 0 = 

2 96 in 

■H41HW 

A = 

98 66360224 in*2 
0 06365381 ^2 

110 3023 
0 071163 

A ave = 

0 06740822 m*2 



0.S9298S768 


Plb a 

1,634.150.0 Pa 

Ttb = 

1.325 2 K 

Pb = 

1,299 940 8 Pa 

Tb = 

1.250 4 K 


3 601 


m dot = 

4 8 8348 Kq/s 

uxb - 

20035 nVe 


Station C 
HPT1 rotor exit 

t = 

0 mm 

Do = 

2352136781 in 

Di = 

19 25626812 in 

delta D = 

4 263099693 in 

A = 

143 238089 in" 2 
0 092410195 m*2 


Me = 

BL'Jf-HW -7*. » BHBB 


Ptc = 

1,034.667 2 Pa 

Tic = 

K 

Pc = 

985.627 1 Pa 

Tc = 

1.1855 K 


EZ 


2 880 Kty'm*3] 


m dot = 

46 6348 kg/a 

uxc = 

iae.74 m* 




0da fl 1 




QQSBD 


niwrffB 

■Mill ■ 1 — 

fiassma 



Ma = 

0 192941463 

delta Ma = 

0 00000067 





■,vx.*i(f!nv7rKmiH 


mmEEwssm 


mi 





'imi 

«■ 


Vb = 

412.40 nv» 

vb = 

380.48 nVe 

Mb = 

0 592963765 

delta Mb = 

1 0.00000000 



Mlltl 

qsxhhbbh 

■nHHPv'H jJTTTWi 




MBHiiF' E E3BB 




— JSMgm 










HBHKEBXSuHi 





TtUBSEHHHI 

■ ihiiini ■ 




Vc = 

304 90 nV« 1 

VC = 

24408 r»Wa 1 


VC = 


Vc « 


Me = 

0 270039781 

delta Me * 

4>.ooooqqot 


5Z 3 


Stage 2 



EE3BEOZBB1 


BDtUSQI^^I 

CXSEHBi 

mrrxrrrvrm 


MBSSnEiTZMl 






Stage 1 



^K&CE377Bi 



EHB^^BB 

■KUXSH 


B^H3£EBB 


Idelta(W) a [ 0W I 


tl poly * 

09 

T 11 * 

0 879105 

X tl = 

0 633153 


Table C.l Data for the maximum power 2° deviation velocity triangle diagram of Figure C. 

In an attempt to obtain a more accuate estimate of the deviation an empirical correlation, Carter's Rule, was 
employed as follows [36]: 


<!>- P\+ Pi 

c 

a = — 

s,=-A= 

8 -Jo 

This resulted in a deviation of 7.5° and 7.3° for the NGV and rotor, respectively. Using these deviation values, 
the rotor incidence angle was 7.0° and the rotor exit angle was -14.8°. 

Another possibility to increase the accuracy of the velocity trianlgle analysis was investigated. A total pressure 
loss coefficient (“profile loss”) at Station B was attempted rather than assuming a negligible total pressure loss 
across the NGV. However, using empirical cascade correlations for a stator which are a function of the solidity 


(C.23) 

(C.24) 

(C.25) 








(a = c dx j a i/s) the pressure loss coefficient was found to be w = 0.03 and the effect on the result for P-^ and (3b was 
found to be less than 27c [36]. 

Overall, the velocity triangle analysis provided reasonable fluid property data for the high fidelity modeling. 
Including the profile loss had a negligible effect on the results and thus it was ignored to keep the analysis 
simple. Also, it is difficult to justify a particular deviation angle. The NASA/DERA engine test simultaions 
used a deviation angle of 2° since it was thought to be a conservative approach and gave a nearly axial exit flow 
angle. As another point of reference, the HPT2 NGV has a chord angle of 41° versus 45° for the HPT1 NGV. 
This would imply that the blades are both roughly designed to have a zero inlet swirl, supporting the notion that 
the exit flow angle be nearly axial. 
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Appendix D 


Detailed Results from NASA/DERA Engine 
Simulation of Max Power Case 

This appendix contains a complete set of calculation results for the max power condition. From the high 
fidelity modeling effort the following items are given in order for the mixed-out case, first for the NGV and 
then for the rotor: 

• Fluid and chemistry convergence history 

• Range and mass averaged quantities of fluid and chemistry variables at exit plane 

• Total temperature and total pressure contours 

• Mass fraction contours of trace species 

Then, for the high fidelity rotor wake model case: 

• Time series of S0 3 mass fraction for one blade passing cycle (see Appendix D for an animation) 

And finally, from the low fidelity modeling effort the following items are given: 

• SO* species evolution 

• NO y species evolution 

• CO, OH, and O species evolution 

The CNEWT and CALCHEM input file used to obtain these results are given in Appendix E. 
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range and mass averaged quantities at cutting plane 
cutting plane at 0.04515 


Num 

Spec 

min 

max 

mass ave 

1 

CO 

1 . 0960670E-05 

1 .1788260E-05 

1 . 1291178E-05 

2 

C02 

7 . 2901480E-02 

7 . 2902328E-02 

7 . 2902240E-02 

3 

02 

0.1480127 

0.1480136 

0.1480131 

4 

H20 

2 . 9883171E-02 

2 . 9883521E-02 

2 . 9883305E-02 

5 

H 

2 . 8244 520E-12 

9 . 4324088E-12 

8 . 3788480E- 12 

6 

0 

2.6128790E-08 

4 . 4815678E-08 

3.4125186E-08 

7 

OH 

1.9860620E-06 

3 . 9784268E-06 

3.6788304E-06 

8 

H02 

2 . 4688800E-08 

5 . 8078822E-08 

4 . 2947221E-08 

9 

H202 

5 . 44 99472E-09 

9 . 8093 857E-08 

1 . 6948434E- 08 

10 

NO 

3 . 1216681E-04 

3 . 1549390E-04 

3 . 1493246E- 04 

11 

NO 2 

4 . 9967889E-06 

6 . 4257251E-06 

5 . 7360417E-06 

12 

NO 3 

5 . 9383180E- 12 

3 . 1245451E-11 

1 . 0588782E - 11 

13 

HNO 

7 . 8105612E-11 

1.7963479E-10 

1 . 2943577E- 10 

14 

HONO 

4 . 0774090E-08 

4 . 1504150E-06 

2 . 0222627E- 07 

15 

HNO 3 

- 1 . 4242200E-10 

5 . 3347769E-08 

9 . 1435254E- 10 

16 

S02 

2 . 1805350E-05 

2 . 1853461E- 05 

2 . 1 843829E- 05 

17 

S03 

1 . 744 4 370E-06 

1 .7997299E-06 

1.7563391E-06 

18 

HOSO 

2.7650279E-14 

2 . 5342201E-13 

1 . 44 87009E - 1 3 

19 

HS03 

5 . 0049240E- 12 

1 .1248480E 10 

8. 1400182E-12 

20 

HS02 

1 .2148050E-16 

5 . 7409781E-16 

1 . 5103668E- 16 

21 

SO 

9 . 1956160E- 14 

3 . 9928561E- 13 

3 . 4998060E- 13 

22 

HCO 

2. 5679660E-17 

5 . 7935929E- 17 

4 . 5071751E- 17 

23 

H2 

1 . 1504340E-10 

1 .0300690E-09 

7 . 3 81 991 IE- 10 

24 

N2 

0.7488416 

0.7488421 

0.7488415 

25 

AR 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

26 

C(S) 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

27 

N 

0.0000000E+00 

O.OOOOOOOE+OO 

0 . OOOOOOOE+OO 

28 

H2S04 

3 . 7797829E-10 

6.7893300E 09 

5. 722225 5E-10 

29 

S 

2 . 4445289E-21 

2 . 4445451E-21 

2 . 4445346E-21 


X/X 

7 . 1622198E-03 

5 . 9027892E-02 


Y/R 

0.2616538 

0.2712400 


Z/T 

- 6 . 3 572727E- 02 

2 . 0361030E- 02 


Tt 

1156.810 

1386.379 

1363.339 

dens 

3.119411 

4.331536 


Ps 

10.55461 

16.76656 

12 . 99278 

Vx 

-123.4683 

450.0187 

201.4363 

Vt 

-430.1193 

261.4255 

361.1647 

Vr 

-84 . 10556 

25.99604 

- 5 . 2535836E-02 

TK 

0.5282233 

285292.6 


TE 

12 . 38410 

6 . 1178208E+10 


Ts 

1083.624 

1374.872 

1287.489 

Mach 

3 . 7381984E-03 

0.9269373 


Pt 

13.64641 

21.36677 

16.27730 

Pred 

10.55461 

16.76656 


S 

0.6861496 

0.8642535 
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Y(OH): 1.9861e-6 to 3.9784-6 Y(H02): 2.4689e-8 to 5.8079e-8 Y(H202): 5.4499e-9 to 9.8094e-8 

MassAve.: 3.6788e-6 MassAve.: 4.2947e-8 MassAve.: 1.6948e-8 
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VPFR operating point 

intersection of red lines (LD = 15, AR = 13.7) 
Fully developed two dimensional stall 

Initial separation form computation 

intersection of blue lines (UD - 8.5, AR - 2.5) 
Near line of appreciable stall 


Figure A.4 Diffuser flow regime map [39] 


Velocity Profiles 

The most informative method to evaluate the quality of the computational solution was to compare the velocity 
profiles to experimental data and empirical correlations. The VPFR was designed to operate such that the 
velocity profiles transition quickly to fully-developed turbulent velocity profiles [26]. These profiles have a 
characteristic “top hat” shape (see Figure A. 5 and Figure A.6). 

Empirical Correlation 

The Power Law is one empirical correlation commonly used to predict velocity profiles in pipe flow: 


u 

K 


1 -- 


R 


(A3) 


Where V c is the centerline velocity, u is the time average x component of velocity, and R is the pipe radius. For 
the flow reactor case, n, which is a function of Reynolds number, is about 6. Figure A.5 shows the expected 
velocity profile for the VPFR derived from the Power Law. 
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Y(S02): 2.1805e-5 to 2.1853e-5 Y(S03): 1 .7444e-6 to 1 ,7997e-6 Y(HOSO): 2.7650e-14 to 2.5342e-13 

MassAve.: 2.1844e-5 MassAve.: 1.7563e-6 Mass Ave.: 1.4487e-13 
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Y(HNO): 0.5744e-10 to 0.2084e-9 Y(HONO): 0.7931 e-7 to 0.4635e-6 Y(HN03): 0.1836e-9 0.5146e-8 

Mass Ave.: 1 .487081 6e-10 Mass Ave.: 2.3786764e-7 MassAve.: 1.1930059e-99 



Y(S02): 0.2185e-4 to 0.2185e-4 Y(S03): 0.1732e-5 to 0.1750e-5 Y(HOSO): 0.2800e-14 to 0.2575e-13 

Mass Ave.: 2.1851994e-5 MassAve.: 1.7476095e-6 MassAve.: 1 ,4743281e-14 
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range and mass averaged quantities at cutting plane 
cutting plane at 0.01060 


num 

species 

min 

1 

CO 

1.1028110E-05 

2 

C02 

7 . 2902203E-02 

3 

02 

0.1480128 

4 

H20 

2 . 9883269E-02 

5 

H 

3 .0306G2CE-12 

6 

O 

3 . 0920891E-08 

7 

OH 

1 . 8885230E-06 

8 

H02 

2. 373809 IE- 08 

9 

H202 

1 . 0050790E-08 

10 

NO 

3 . 1211879E-04 

11 

N02 

5 . 7360421E-06 

12 

NO 3 

8 . 3264983E-12 

13 

HNO 

1 . 2877401E -10 

14 

HONO 

8 . 1516284E--08 

15 

HNO 3 

-3 . 1286479E-10 

16 

S02 

2 . 1808110E -05 

17 

S03 

1. 756339 QE- 06 

18 

HOSO 

2 . 1189459E-14 

19 

HS03 

5.8176632E-12 

20 

HS02 

1 .3123580E-16 

21 

SO 

8. 1447761E-14 

22 

HCO 

2 . 3288410E-17 

23 

H2 

7.3819911E-10 

24 

N2 

0.7488415 

25 

AR 

0. OOOOOOOE+OO 

26 

C{S) 

0. OOOOOOOE+OO 

27 

N 

0. OOOOOOOE+OO 

28 

H2S04 

4 . 2082671E-10 

29 

S 

2. 444528 9E -21 


max mass ave 

1.1291180E-C5 1.1155897E-05 

7 . 2902679E-02 7 . 2902434E • 02 

0.1480133 0.1480130 

2 . 9883681E-02 2 . 98834 01E-02 

1 . 1951770E-11 9 . 4046255E-12 

5. 6032619E-08 4 . 3273968E-08 

3 . 7138709E-06 3 . 3753893E-06 

4 . 5069211E-08 3 . 68824 4 0E- 0 8 

8 . 7000387E-08 3 . 9920529E-08 

3 . 1495970E-04 3 . 1454579E - 04 

6 . 32113 10E-06 5 . 9863410E-06 

3 . 6301271E-11 1 . 4460729E-11 

1 . 8293370E-10 1 . 550404 0E - 10 

4 . 033 3 512E-06 5 . 5106256E-07 

6. 190 371 2E- 08 2 . 4 6052 82E-09 

2 . 184 3 83 IE- 05 2 . 1838596E-05 

1 .7951180E-06 1 . 7627525E -06 

1 . 6339089E-13 9 . 3304000E - 14 

1 . 5686860E-10 1 . 0739637E-11 

7 . 0853009E-1 6 1 . 9591 100E - 16 

3 . 6073450E-13 3. 118533 8E -13 

4 . 5266589E-17 4 . 0807966E- 17 

8 . 7 823 260E- 1 0 8 . 1180129E-10 

0.7488420 0.7488419 

0 . OOOOOOOE+OO 0 . OOOOOOOE+ 00 
0 . OOOOOOOE+OO 0 . OOOOOOOE+OO 
0 . OOOOOOOE + OO 0 . OOOOOOOE+OO 
9 . 0190957E -09 7 . 570201 IE - 1 0 

2 . 4445469E -21 2 . 4445352E-21 


x/x 

-1 . 1852840E -02 

2 . 2079229E- 02 


Y/R 

0.2671618 

0.2712400 


T/Z 

-1 . 2139130E-G2 

3.3451378E-02 


Tt 

1124.372 

1333.011 

1300.225 

dens 

2.046658 

3.929755 


Ps 

5.856664 

14.15914 

9.808423 

Vx 

-72.06285 

712.5885 

271.0783 

Vt 

-484.9583 

490.9988 

-379.8698 

Vr 

-52.12230 

43.46152 

-0.2326850 

TK 

0.5275032 

527503.2 


TE 

20.16434 

6 . 5262518E+11 


Ts 

851.8360 

1315.987 

1202.610 

Mach 

3 . 7721284E-03 

1.512067 


Pt 

9.385468 

23.24667 

13.33437 

Pred 

5.856664 

14.15914 


S 

0.6160989 

0.8636632 
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Y(HNO): 1.2877e-10to 1.8293e-10 Y(HONO): 8.1516e-8 to 4.0334e-6 Y(HNO3):-3.1286e-10 to 6.1904e-8 

Mass Ave.: 1.5504e-10 MassAve.: 5.5106e-7 MassAve.: 2.4605e-9 



Y(HS03): 5.8177e-12 to 1.5687e-10 Y(HS02): 1 3124e-16 to 7.0853e-16 Y(SO): 8.1448e-14 to 3.6073e-13 

MassAve.: 1.0740e-11 MassAve.: 1.9591e-16 MassAve.: 3.1185e-13 
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DERA Engine Exft Plane Calculation For Max Power Case 


4 SOx for HPT1 

| 

q i 

* r 


10" 

-7 

10 



10 '" 



0 15 02 0.25 

Time [ms) 


0.35 


SOx for HPT1 Exit to Nozzle Exit 




V 


\ 




— - S02 

— S03 
H2S04 

- NS03 

— so 




High Temperature Streamline 
Low Temperature Streamline 


0.5 1 1.5 2 2.5 3 

Time (mi) 


NOTE: The legend on the chart lists the species in the proper sequence for the lines (top to bottom), however, 
the line style is not correct (dashed, dotted, solid, etc is wrong). 
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DERA Engine Exit Plane Calculation For Max Power Case 


NOv tor HPT1 



0 0.05 0.1 0 15 0.2 0.25 0 3 0.35 

Time (ms) 


NOy for HPT1 Exit to Nozzle Exrt 



0 0.5 1 1-5 2 25 3 


Time (ms) 


NOTE: The legend on the chart lists the species in the proper sequence for the lines (top to bottom), however, 
the line style is not correct (dashed, dotted, solid, etc is wrong). 
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DERA Engine Exit Plane Calculation For Max Power Case 


CO. OH. and O for HPTl 



0 0.05 0 1 0 15 0 2 0 25 0.3 0 35 

Time (ms) 


CO, OH, and O for HPTl Exit to Nozzle Exit 



0 0.5 1 1.5 2 2 5 3 


Time (ms) 


NOTE: The legend on the chart lists the species in the proper sequence for the lines (top to bottom), however, 
the line style is not correct (dashed, dotted, solid, etc is wrong). 
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Appendix E 


CNEWT and CALCHEM Input Files For Max Power 
Condition 

This appendix contains the CNEWT run files used for the high fidelity modeling and CALCHEM input files 
used for the low fidelity modeling of the NAS A/DERA engine test max power condition. 


There are three CNEWT run files listed in the following order: 

• HPT1 NGV - Max pow-er 

• HPT1 rotor - Max power - Mixed-out 

• HPT1 rotor - Max power - Unsteady wake model 

Followed by four CALCHEM input files listed in the following order: 

• HPT1 - Max power - Mass averaged 

• HPT1 - Max power - High temperature streamline 

• HPTi - Max power - Low temperature streamline 

• HPTI exit to nozzle exit - Max power - Mass Averaged 
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CNEWT: mid -span of DERA HPTl NGV:max power condition 

3550 0 2 1 1 0 < change KNORM=l for stronger blade be 

1.00000 0.00000 0.50000 0.00200 0.25000 

1634150. 1365.20000 00.00000 0.00000 000.00000 1299940.8 1299940.8 998.99700 

1138. 1.34 0.0000546 0000.0 0.031 

1 2 3 3 -4 

3 o 1 29 1 1 < MASKIN=0; code will not read lines after MIN 

25 80 3 


0.26924 

0.27124 


1634150.0 

1634150.0 

1634150 

1365.2 

1365.2 

1365. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1 . 

1 . 

1 . 


1. 178826 Oe- 05 1 . 1788260e- 05 1 . 1788260e-05 
7 . 2901490e -02 7 . 290 1490e- 02 7 . 2901490e-02 
1 . 4801360e-01 1 . 48013 60e-01 1 . 4801360e-01 
2. 9883310e-02 2 . 98833 lOe- 02 2 . 9883310e -02 
7. 8778370e-12 7 . 8778370e-12 7 . 8778370e-12 
2 . 6128790e - 08 2 . 6128790e- 08 2 . 6128790e-08 
3 . 743 8540e- 06 3.7438540e 06 3 . 7438540e-06 
5 . 4143820e- 08 5 . 4143820e- 08 5 . 4143820e-08 
6 . 6700590e-09 6.6700590e 09 6 . 6700590e-09 
3 . 1549390e-04 3 . 1549390e- 04 3 . 1549390e-04 
4 . 9967900e- 06 4 . 9967900e- 06 4 . 9967900e-06 
6 . 0388730e-12 6.0388730e 12 6 . 03 8873 0e-12 
7 . 8105610e- 11 7 . 8105610e- 11 7 . 8105610e-ll 
7 . 8798570e-08 7 . 8798570e- 08 7 . 8798570e-08 
3 . 0977400e-10 3 . 0977400e- 10 3 . 0977400e-10 
2 . 1853460e-05 2 . 1853460e-05 2 . 1853 460e-05 
1 . 7444370e-06 1 . 7444370e - 06 1 . 7444 370e-06 
2. 291683 0e- 13 2 . 2916830e- 13 2 . 2916830e-13 
5 . 94 82370e - 12 5 . 9482370e- 12 5 . 94 82370e-12 
1.2286700e-16 1 . 2286700e - 16 1 . 2286700e-16 
3 . 603 001 Oe-- 13 3.6030010e 13 3 . 6030010e-13 
5 . 5508550e - 17 5 . 5508550e- 17 5 . 5508550e-17 
1 . 1504340e-10 1.1504340e 10 1 . 1504 340e-10 
7 . 4884180e-01 7 . 4884180e-01 7 . 4884180e -01 
0 .0000000e+00 0 . 0000000e+00 0 . 0000000e+00 
0 . 0000000e+00 0 . 0000000e+00 0 . 0000000e+00 
0 . 0000000e+00 0 . 0000000e+00 0 . 0000000e+00 
4 . 3699380e-10 4 . 3699380e-10 4 . 3699380e- 10 
2.4445370e-21 2.4445370e 21 2 . 4445370e-21 


-0.028 

0.343 

0. 

< boundary 

staionary wrto. 

blade 


1 . 

1 . 

1 . 

1 . 1 . 

0.999 

0.995 

0.992 

0.99100 

0.992 

0.995 

0.999 1. 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 1 . 

0.999 

0.995 

0.992 

0.99100 

0.992 

0.995 

0.999 1. 

1 . 

1 . 

1 . 

1 . 
5. 

5. 

4 . 

3. 2. 

2. 

2 . 

2. 

2. 

2. 

2 . 

2. 3. 

4. 

5. 

5. 

5. 

5. 

4 . 

3. 2. 

2. 

2. 

2. 

2. 

2. 

2 . 

2. 3. 

4. 

5 . 

5. 
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CNEWT: mid- span of DERA HPT1 rotor: max power condition, mixed- out 

4800 0 2 1 1 0 < change KNORM=l for stronger blade be 

1.00000 0.00000 0.50000 0.00200 0.25000 

1627730.0 1363.33900 -60.84980 0.00000 361.16470 985627.1 985627.1 956.20700 

1138. 1.34 0.0000546 -8863.8 0.019 

1 2 3 3 4' 

3 o 129 11 < MASKIN=0 ,* code will not read lines after MIN 

25 80 3 

0.26924 0.27124 

1627730.0 1627730.0 1627730.0 
1363.339 1363.339 1363.339 

-60.850 -60.850 -60.850 

0. 0. 0. 

361.165 361.165 361.165 

1. 1. 1. 

1.1291178E-05 1.1291178E-05 1 . 1291178E-05 
7 . 2902240E-02 7 . 2902240E-02 7 . 2902240E-02 
0.1480131 0.1480131 0.1480131 

2 . 9883305E-02 2 . 98833 05E-02 2 . 9883305E-02 
8 . 3788480E-12 8 . 3788480E- 12 8 . 3788480E-12 
3 . 4125186E-08 3 . 4125186E-08 3 . 4 125186E-08 
3 . 6788304E-06 3 . 6783304E-06 3 . 6788304E-06 
4 . 2 947221E-08 4 . 2 947221E- 08 4 . 2947221E- 08 
1 . 694 8434 E-08 1 . 6948434E-08 1 . 694 8434E- 08 
3.1493246E-04 3 . 14 93 246E- 04 3 . 1493246E- 04 
5 . 7360417E-06 5 . 736041 7E- 06 5 . 73 604 1 7E- 06 
1 . 0588782E-11 1 . 0588782E- 1 1 1 . 0588782E- 1 1 
1 . 2943577E-10 1 . 2 943 577E - 1 0 1 . 2943577E-10 
2 . 0222627E-07 2 . 0222627E- 07 2 . 0222627E-07 
9. 1435254E-10 9 . 143 5254E - 10 9 . 14 3 5254E-1 0 
2 . 1 84 3829E- 05 2 . 184 3 329E- 05 2 . 184 3 82 9E- 0 5 
1 .7563391E-06 1 . 7563 3 91E - 06 1 . 7563 3 9 IE- 06 
1 . 4487009E-13 1 . 44 87 009E- 1 3 1 . 44 87009E-13 
8.1400182E-12 8 . 1400182E-12 8 . 1400182E-12 
1 . 5103 668E- 1 6 1 . 5103668E-16 1 . 5103668E -16 
3 . 4998060E-13 3 . 4998060E-13 3.4998060E 13 
4.5071751E-17 4 . 5071751E-17 4 . 507175 IE - 17 
7 . 3819911E-10 7 . 3819911E-10 7 . 3819911E-10 
0.7488415 0.7488415 0.7488415 

0. OOOOOOOE+OO 0 . OOOOOOOE+OO 0 . 0000000E+00 
0 . OOOOOOOE+OO 0 . OOOOOOOE+OO 0 .OOOOOOOE+OO 
0 . OOOOOOOE+OO 0 . OOOOOOOE+OO 0 .OOOOOOOE+OO 
5 . 7222255E-10 5 . 72222 55E- 10 5 . 72222 5 5E- 1 0 
2 . 4445346E-21 2 . 4445346E-21 2 . 4445346E-21 


-0.028 

0.343 

0. 

< boundary 

staionary wrto. 

blade 


1 . 

1 . 

1 . 

1. 1 . 

0.999 

0.995 

0.992 

0.99100 

0.992 

0.995 

0.999 1. 

1 . 

1 . 

1 . 

1. 

1 . 

1 . 

1. 1 . 

0.999 

0.995 

0.992 

0.99100 

0.992 

0 . 995 

0.999 1. 

1 . 

1 . 

1. 

1 . 
5. 

5. 

4 . 

3. 2. 

2. 

2 . 

2. 

2. 

2. 

2. 

2. 3. 

4 . 

5. 

5. 

5. 

5. 

4 . 

3. 2. 

2. 

2 . 

2. 

2. 

2. 

2. 

2. 3. 

4 . 

5. 

5. 


201 



HPT1 rotor - Max power - Unsteady wake model: 


CNEWT: mid -span af DERA HPTl rotor:max power conditions pitch inlet wake 

10000 0 2 1 1 0 < change KNORM=l for stronger blade be . 

-2.01000 0.00000 0.50000 0.00200 0-25000 

1627730.0 1363.33900 -60.34980 0.00000 361.16470 985627.1 985627.1 956.20700 

1138. 1.34 0.0000546 -8863.3 0.019 

1 2 3 3 -4 

1 1 < -- MASKIN=70, MIN=0 


121 75 3 

-0.042507 0.269035 8863.8 

1 . 6319600E+06 1 . 633 40 97E^06 1 . 6344 153E+06 1 . 634 82 85E+06 
1 . 6458172E+06 1 . 64 306 68E+ 06 1 . 644 S275E+06 1 . 6254299E+06 
1 . 5864542E+06 1 . 613 6267E+06 1 . 6284102E+06 1 . 633 9775E+06 
1 . 6334097E+06 1 . 63441 53E+06 1 . 6348285E+06 1 . 63 3 9272E+06 
1 . 6480668E+06 1 . 644 527SE+06 1 . 6254299E+06 1 . 60193 41E+06 
1 . 6136267E+06 1 . 62 84 102E+06 1 . 63 3 9775E+06 1 . 634 004 1E+06 
1 . 6344153E+06 1 . 63 4 8285E> 06 1 . 63 3 9272E+06 1 . 6326136E+06 
1 . 6445275E+06 1 . 6254299E+06 1 . 6019341E+06 1 . 5646283 E+06 
1 . 6284102E+06 1 . 63 3977 5E+06 1 . 6340041E+06 1 . 6323421E+06 
1 . 3733000E+03 1 . 3736851E+03 1 . 3734623E+03 1 . 373 3692E+03 
1 . 3602751E+03 1 . 3468093 E+03 1 . 3327657E+03 1 . 3 123256E+03 
1 . 3204067E+03 1 . 3591487E+03 1 . 3768992E+03 1 . 3778210E+03 
1 . 3736851E+03 1 . 373 4 623 E+03 1 . 373 3692E+03 1 . 37447 91E+03 
1 . 3468093E+03 1 . 33 27657E+03 1 . 3123256E+03 1 . 2930020E+03 
1 . 3591487E+03 1 . 376 8 992E+03 1 . 377 82 10E+03 1 . 37513 60E+03 
1 . 3734623E+03 1 . 373 3692E+ 03 1 . 374 4 7 91E+03 1 . 376353 OE+03 
1 . 3327657E+03 1 . 3 1232 56E+03 1 . 293 0020E+03 1 . 2666420E+03 
1 . 3768992E+03 1 . 377 82 10E+03 1 . 3751360E+03 1 . 373 5196E+03 
1 . 1344700E-05 1.1317641E 05 1 . 128557SE-05 1 . 125S040E- 05 
1 . 1116077E-05 1 . 1103603 E 05 1 . 1099339E- 05 1 . 1104153E-05 
1 . 1350943E- 05 1.1357249E 05 1 . 136 3 4 10E - 05 1 . 1372702E- 05 
1 . 1317641E-05 1 . 1285575E -05 1 . 1255040E-05 1 . 122 37 97E- 0 5 
1 . 11036C3E-OS 1 . 1099339E-05 1 . 1104 1 53E- 05 1 . 1146548E- 05 
1 . 1357249E-05 1.1363410E 05 1 . 1372702E-05 1 . 1 37 64 3 8E-05 
1 . 1285575E-05 1.1255040E 05 1 . 1 22 37 97E- 05 1 . 11947 30E- 05 
1 . 1099339E-05 1.1104153E 05 1 . 114 6548E- 05 1 . 11973 13E-05 
1 . 1363410E-05 1 . 1372702E- 05 1 . 137 64 38E - 05 1 . 1374067E-05 
7 . 2902292E-02 7.2902292E 02 7 . 2902292E 02 7 . 2 9022 92E- 02 
7.2902292E-02 7.2902292E 02 7 . 2902292E-02 7 . 2902292E- 02 
7 . 2902292E-02 7.2902292E 02 7 . 2902292E-02 7 . 2 902292 E-02 
7 . 2902292E-02 7.2902292E 02 7 . 29022 92E- 02 7 . 2 902292 E-02 
7 . 2902292E-02 7.2902292E 02 7 . 2902292E- 02 7 . 2902292 E-02 
7 . 2902292E-02 7 . 2902292E • 02 7 . 2 902292 E - 02 7 . 2902292 E-02 
7 . 2902292E-02 7.2902292E 02 7 . 2902292E-02 7 . 29022 92 E-02 
7 . 2902292E-02 7.2902292E 02 7 . 2902292E-02 7 . 29022 92 E- 02 
7 . 2 902292E-02 7.2902292E 02 7 . 2902292E- 02 7 . 2 9022 92E-02 
1 .4801300E-01 1.4801300E 01 1 . 4801300E-01 1 . 4 801300 E- 01 
1 .4801300E-01 1 . 4801300E 01 1 . 4 801300E- 01 1 . 4 801300E- 01 
1.4801300E 01 1.4801300E 01 1 . 4801300E 01 1 . 4 801300E-01 
1 .4801300E-01 1.4801300E 01 1 .4801300E-01 1 . 4 801300E-01 
1 .4801300E-01 1.4801300E 01 1 . 4 8013 OOE- 01 1 . 4 801300E- 01 
1 . 4801300E-01 1.4B01300E 01 1 .4801300E-01 1 . 4801300E-01 
1 .4801300E-01 1 . 4801300E- 01 1 . 4801300E-01 1 . 48013 OOE- 01 
1 . 4801300E-01 1 . 4 8 013 00E 01 1 . 48013 00E-01 1 . 48013 OOE- 01 
1 .4801300E-01 1 .4801300E- 01 1 . 4801300E-01 1 . 4 8013 OOE- 01 
2 . 9883300E-02 2.9883300E 02 2 . 9883300E-02 2 . 9883300E-02 
2 . 9883300E-02 2 . 98 83300E- 02 2 . 9883300E-02 2 . 9883300E-02 
2 . 9883300E-02 2 . 98 B33O0E- 02 2 . 9883300E-02 2 . 98833 0OE-O2 
2 . 9883300E-02 2.988330GE 02 2 . 9883300E-02 2 . 98833 OOE- 02 
2 . 9883300E-02 2.9883300E 02 2 . 9883300E-02 2 . 98833 OOE-02 
2 . 9883300E-02 2 . 9883300E- 02 2 . 9883300E-02 2 . 9883300E-02 
2 . 9B83300E-02 2 . 98833O0E 02 2 . 98833 00E-02 2 . 9883300E- 02 
2 . 988330GE-02 2 . 98833O0E- 02 2 . 98833 Q0E -02 2 . 98833 OOE-02 
2 . 9883300E-02 2 . 9883300E- 02 2 . 98833 00E-02 2 . 98833 OOE- 02 
8 . 4637200E-12 8 . 5439389E- 12 8 . 5969422E-12 8 . 63545 14E- 12 
8 . 6063023E-12 8 . 4093940E- 12 8 . 2426977E - 12 7 . 9595924E-12 
7 . 8549552E- 12 8 . 1559254E- 12 8 . 2472808E-12 8 . 23646 94E-12 
8 . 5439389E- 12 3.5969422E 12 8 . 63545 14E- 12 8 . 6 S83739E- 12 
8.4093940E-12 8 . 2426977E- 12 7 . 9595924E - 12 7 . 57057 33 E- 12 
8 . 1559254E-12 8.2472808E 12 8 . 2364694E-12 8 . 2 33 24 6 9E- 12 
8 . 5969422E- 12 8 . 63 54S14E- 12 8 . 6883739E-12 8 . 7451264E-12 
8.2426977E-12 7 . 9595924E 12 7 . 5705733E-12 7 . 3 1634 17E- 12 
8 . 2472808E-12 8 . 2364694E- 12 8 . 23324 69E- 12 8 . 2851215E-12 

3 . 3 813600E- 08 3 . 4291738E- 08 3 . 4 6423 3 IE 08 3 . 4 864049E-08 
3 . 6222523E-08 3 - 6740560E- 08 3 .7039767E-08 3 . 72813 42E- 08 
3 . 3929308E-08 3 . 2220627E- 08 3 . 1432577E-08 3 . 1 552902E- 08 
3 . 4291738E-08 3 . 4 6423 31E- 08 3 . 4864049E-08 3 . 50178 3 9E- 08 


1 . 6339272E+06 1 . 6326136E+06 1 . 6342565E+06 1 . 6400197E+06 
1.60193 41E+06 1 . 5646283E + 06 1 . 55303 12E+06 1 . 5615181E+06 
1 . 6340041E+06 1 . 6323421E+06 1 . 63 2 1 173E+06 1 . 63 19600E+06 
1 .6326136E+06 1 . 6 3 42 56 5E+0 6 1 . 64001 97E+06 1 . 6458172E+06 
1 .5646283E+06 1 . 5 5 3 0 3 12E+0 6 1 . 561 51 81E+06 1 . 5864542E+06 
1 . 6323421E+06 1 . 63 21 17 3 E+0 6 1 . 631 9600E+06 1 . 6334097E+06 
1 . 6342565E+06 1 . 6400197E+06 1 . 6458172E+06 1 . 6480668E+06 
1.55303 12E+06 1 . 5615181E+06 1 . 5864542E+06 1 . 6136267E+06 
1 . 6321173E+06 1 . 6 3 1 9 600E+0 6 

1.37447 91E+03 1 . 3763530E+03 1 . 3763 4 3 5E+03 1 . 3711324E+03 
1 . 2930020E+03 1 . 2 6 66 42 0E+0 3 1 . 263 67 17E+03 1 . 2829062E+03 
1 . 37513 60E+03 1 . 3 7 3 5 1 96E + 0 3 1 . 37 3 1996E + 03 1 . 373 3 OOOE+03 
1.376353 0E+03 1 . 3 7 63 4 3 5E+0 3 1 . 3711 324E+03 1 . 3 602751E+03 
1 . 2666420E + 03 1 . 263 67 17 E + 0 3 1 . 2829062E + 03 1 . 3204067E+03 
1 .3735196E + 03 1 . 3 7 3 1 996E + 0 3 1 . 37 3 3 QO0E+O3 1 . 373 6 851E+03 
1 . 3763435E + 03 1 . 3 7 11 3 2 4 E + 0 3 1 . 3 60275IE+03 1 . 3468093E + 03 
1 .2636717E+03 1 . 2829062E+03 1 . 3204 067E+03 1 . 3 5914 87E+03 
1.37319 96E + 03 1 . 3733 OOOE + 03 

1 . 1223797E-05 1 . 1 1 947 3 0E - 05 1 . 1 16067QE -05 1 . 114 1 8 55E - 05 
1 . 1146548E- 05 1 . 11 97 3 1 3E - 05 1 . 1260912E -05 1 . 13 24297E- 05 
1 . 1376438E-05 1 . 1 3 7 4 067E - 0 5 1 . 13 63 370E-05 1 . 1344700E-OS 
1 . 1194730E-05 1 . 1160670E-05 1 . 1141855E-05 1 . 1116077E-05 
1 . 1197313E-05 1 . 1260912E-05 1 . 1324297E -05 1 . 1350943E-05 
1 . 1374067E-05 1 . 1 3 63 3 7 0E - 0 5 1 . 1344700E-05 1 . 13 1764 IE -05 
1 . 1160670E-05 1 . 114185SE -05 1 . 1116077E-0S 1 . 1 103 603E -05 
1 .1260912E -05 1 . 1324297E-05 1 . 13 5094 3E-05 1 . 13 57249E -05 
1 . 1363370E-05 1 . 1 3 4 4 7 OOE- 05 

7 . 2902292E- 02 7 . 2 9 0 2 2 92E - 02 7 . 2902292E-02 7 . 2902292E-02 
7 . 2902292E-02 7 . 2 90 2 2 92E- 02 7 . 2902292E- 02 7 . 2 902292E-02 
7 . 2902292E-02 7.2902292E -02 7 . 2902292E 02 7 . 2902292E-02 
7 . 2902292E -02 7 . 2 902 2 92E- 02 7 . 2902292E -02 7 . 2902292E-02 
7 . 2902292E - 02 7 . 2902292E-02 7 . 2902292E- 02 7 . 2902292E -02 
7 . 2902292E -02 7 . 2902292E-02 7 . 2902292E- 02 7 . 2902292E- 02 
7 . 2902292E- 02 7 . 2902 292E- 02 7 . 2902 292E- 02 7 . 2902292E-02 
7 . 2902292E- 02 7 . 29 02 292E - 02 7.2902292E 02 7 . 2 902292E-02 
7 . 2902292E- 02 7 . 2902 292E- 02 

1 .4801300E-01 1 . 4801300E-01 1 . 4801300E-01 1 . 4 801 300E- 01 
1 . 4801300E-01 1 . 4801300E-01 1 . 4801300E -01 1 . 4801300E-01 
1 . 4801300E-01 1 . 4801300E-01 1 . 4801300E-01 1 . 4801300E-01 
1 . 4 8013 OOE - 01 1 . 4801300E-01 1 . 4801300E-01 1 . 4801300E-01 
1 . 4801300E-01 1 . 4801300E-01 1 . 48013 OOE- 01 1 . 4 8013 OOE- 01 
1 . 4801300E-01 1 . 4801300E-01 1 . 4B013 OOE- 01 1 . 4 801 3 OOE- 01 

1 . 4 801 3 OOE- 01 1 . 480130QE- 01 1 . 4801300E- 01 1 - 4 801 300E-01 
1 . 4801300E-G1 1 . 4801300E-01 1 . 4801300E-01 1 . 4 801300E- 01 

1 . 4 8013 OOE- 01 1 . 4801300E-01 

2 . 9883300E-02 2 . 98 8 3 3 OOE- 02 2 . 9883300E- 02 2 . 98833 OOE- 02 
2 . 98833 OOE- 02 2 . 9883 300E-02 2 . 9883300E-02 2 . 98833 OOE- 02 
2 . 9883300E02 2 . 98 8 3 3 OOE- 02 2 . 9883300E-02 2 . 98833 OOE- 02 
2 . 9883300E -02 2 . 98 8 3 3 OOE- 02 2 . 9883300E-02 2 . 98 833 OOE-02 
2 . 9883300E-02 2 . 98 8 3 3 OOE- 02 2 . 9883300E-02 2 . 98833 OOE- 02 
2 . 9883300E-02 2 . 9883 3 00E-02 2 . 98833 OOE- 02 2 . 9833300E-02 
2 . 9883300E-02 2 . 9883 300E-02 2 . 9883 300E - 02 2 . 9883300E-02 
2 . 9883300E-02 2 . 988 3 300E-02 2 . 9883300E-02 2 . 9883300E-02 
2 . 9883300E-02 2 . 9883300E-02 

8 . 6383739E-12 8 . 7 4 5 12 64E- 12 8 . 7 81 1106E -12 8 . 7 5615 93E - 12 
7 . 570 573 3Z- 12 7 . 3163417E 12 7 . 2633166E-12 7 . 52 01166E- 12 
8.2332469E 12 8 . 2851215E-12 8.3733926E 12 8 4637 200E - 12 
8 . 7451264E-12 8 . 7 8 1 1 1 06E - 12 8 . 7561593E-12 8 . 6063 023E-12 
7 .3163417E-12 7 . 2633 166E- 12 7 . 5201166E-12 7 . 8549552E -12 
8 . 2851215E-12 8 . 37 33 926E- 12 8 . 4637200E-12 8 . 54 3 93 89E - 12 
B .7811106E-12 8.7561593E-12 8 . 6063 02 3E -12 8 . 4093 940E- 12 
7 . 2633166E-12 7 . 5201166E-12 7 . 8549552E-12 8 . 1559254E -12 
8 . 3733926E-12 8 . 4637200E-12 

3 . 5017839E-08 3 . 5 1 3 2 84 6E- 0 8 3 . 53 5153 3E-08 3.5603954E 08 
3 . 7219986E- 08 3 . 6967498E- 08 3 . 6469722E-08 3 . 53 094 51E-08 
3 . 1978402E- 08 3 . 2608269E-08 3 . 3260385E- 08 3 . 3813 600E-08 
3 .5132846E-08 3 . 53 51 533E-08 3 . 5603954E - 08 3.6222523E 08 


202 



3 . 6740560E-08 3 . 7 03 9767E - 08 
3 . 2220627E-08 3 . 143 2577E - 08 
3 . 4 6423 3 IE-03 3 . 4864049E- 08 
3 . 7039767E-08 3 . 7231342E- 08 
3 . 1 432577E- 08 3.15529Q2E 38 
3 . 6 878900E- 06 3 . 7 0802 13E- 06 
3 . 7 690451E - 06 3.6930851E 06 
3 . 4920339E-06 3-6311205E 06 
3 . 7080213E-06 3.7277736E 06 
3 . 6930851E-06 3.6294911E 06 
3 . 6311205E-06 3.6818206E 06 
3 . 727773 6E- 06 3.7479838E 06 
3 . 6294911E-06 3.5229465E 06 
3 . 6818206E- 06 3.6812374E 06 
4 . 2494400E- 08 4.2789645E 08 
4 . 4586922E-08 4 . 3546491E- 08 
4 . 0514640E-08 4.2704330E 08 
4 . 2789645E- OB 4.3188823E 08 
4 . 3 546491E -08 4.2736072E 08 
4 . 27 0433 OE- 08 4.3773S99E 08 
4 . 3188823E-08 4 . 3644891E- 08 
4.273 6Q72E- 08 4.1526559E OS 
4.3773 599E-08 4 . 3 552703 E- 08 
1 . 4940000E-08 1 . 5278168E- 08 
1 . 8140504E-08 2 . 1484439E- 08 
2.3321 9 89E- 0 8 1.4368399E 08 
1 . 5278168E-08 1.5490873E 08 
2 . 1484439E- 08 2.4347309E 08 
1 . 4368399E-08 1.1843122E 08 
1 . 5490873E-08 1 . 5513094E ■ 08 
2.4 3473 09E- 08 2 . 9 84 6094E- 08 
1 . 1843122E-08 1.2080676E 08 
3 . 1502400E-04 3.1499689E 04 
3 . 1475S65E- 04 3.1467365E 04 
3 . 1485724E- 04 3.1500566E 04 
3 . 14 9 96 8 9E ■■ 04 3.1496590E 04 
3 . 1467365E-04 3.1461182E 04 
3 . 1 50 0 5 66E - 04 3.1505820E 04 
3 . 1496590E- 04 3.1493705E 04 

3 . 1 4 6 1 1 82E- 04 3 . 14516B2E • 04 
3 . 1505820E- 04 3.1506870E 04 
5.623 9200E- 06 5.6602822E 06 
5 . 9992466E -06 6 . 0721S84E -06 
5 . 7684544E- 06 5.6547121E 06 
5 . 6602822E-06 5 . 70307 55E- 06 
6. 0721584E-06 6.1363353E 06 
5. 6547121E-06 5 . 6017010E 06 
5 . 7030755E- 06 5.7473146E 06 
6. 1363353E-06 6-2009919E 06 
5 . 6017010E-06 5.5844275E 06 
9 . 9303600E-12 1.0133041E 11 
1 . 1930260E-11 1 . 283 3108E 11 
1. 1850625E-11 9.8119282E 12 
1 . 0133041E -11 1 . 03 15558E 11 
1 . 2833108E-11 1 . 3693774E 11 
9 . 8119282E -12 8.9716894E 12 
1 . 0315558E-11 1 . 0459090E -11 
1 . 3693774E-11 1.4720614E 11 
8 . 9716894E-12 8.9745790E 12 
1 . 2691400E- 10 1.2515727E 10 
1 . 2517224E-10 1.3132024E 10 
1.5083575E 10 1.3571490E 10 
1 . 2515727E-10 1.2381082E 10 
1 . 3132024E-10 1.3749699E 10 
1 . 3571490E-10 1 . 283 5028E 10 
1 . 2381082E-10 1.2254941E 10 
1. 3749699E-10 1.4588913E 10 
1 . 2835028E-10 1.2782938E 10 
1 . 7317200E-07 1.7721829E 07 

2. 14 9 08 6 5 E- 07 2.5504971E 07 
2 . 9083002E-07 1.7581928E 07 
1 . 7721829E-07 1.7970502E 07 
2 . 5504971E-07 2.9995857E 07 
1 .7S81928E-07 1.41B5696E 07 
1 . 7970502E- 07 1 . 80042 58E - 07 
2 . 9995857E-07 3.7428040E 07 
1 . 4185696E-07 1.4128161E 07 
7.579 1800E- 10 7 . 811 1140E 10 
1 . 0077276E-09 1.2096770E 09 
1 . 3072919E-09 7 . 77663 58E 10 
7.811 114 OE- 10 7 . 986 1544E 10 
1 . 2096770E-09 1.4373165E 09 
7 .77663 58E- 10 6 . 2276265E - 10 
7.986 1544E- 10 8.0693481E 10 


3 . 7281342E-08 3.7219986E 08 
3 . 1552902E-08 3 . 197 8402E - 08 
3 . 5017839E -08 3.5132346E 03 
3 .7219986E-03 3 . 6967498E- 0 3 
3 . 197 84 Q2E-08 3.2608269E 03 
3 . 72777 36E- 06 3 . 747983 BE-06 
3 . 6294911E-06 '3 . 52294 65E-06 
3 . 6313206E-06 3 . 6812374E-06 
3 . 747983 8E-06 3 . 7766963 E- 06 
3 . 52294 65E- 06 3 . 3766549E-06 
3 . 6312374E-06 3 . 6729106E- 06 
3 .7766963E-06 3 . 807 3 647 E- 06 
3 . 3766549E -06 3 . 2 8 04 4 93 E-- 06 
3 . 6729106E -06 3 . 66893 27E- 06 
4 . 3 188823E-08 4 . 3 644 891E- 08 
4.2736072E-08 4 . 1526 559E- 0 8 
4 . 3 77 3 599E-08 4 . 3 5 52703E - 08 

4 . 3 644S91E-08 4 . 4264 53 OE - 0 8 
4 . 1526559E-08 3 . 9883 5 13E -08 
4 . 3552703E-08 4 . 30S0069E-08 

4 . 4 264530E -0 8 4 . 4 912680E-08 
3 . 9883513E-08 3 . 37 35799E-08 
4 . 3050069E 08 4 . 2603 554E- 08 
1. 5490873E-08 1 . 5513094E--08 
2 . 4347309E -08 2 . 9846094E-08 
1 . 1843122E-08 1 . 2080676E-08 
1 . 5513094E-08 1 . 5253305E-08 
2 . 9846094E 08 3 . 7765007E-08 
1 . 2080676E -08 1 . 2780012E- 08 
1 . 5253305E-08 1 . 4886596E- 08 
3 . 7765007E-08 4 . 263 2635E - 08 
1 . 2780012E-- 08 1 . 36534 57E-08 
3 . 1496590E-04 3.1493705E 04 
3 . 1461182E-04 3 . 14 51682E - 04 
3 . 1505820E-04 3 . 1506870E- 04 
3 . 14937OSE-04 3 . 149144 IE- 04 
3 . 1451682E-04 3 . 1442564E-04 
3 . 15Q6870E-04 3.1506562E 04 
3 . 1491441E-04 3 . 1489664E-04 
3.1442564E-04 3 . 14 39SS6E-04 
3 . 1 506562E-04 3 . 1505710E-04 
5 . 7030755E-06 5 . 7473 146E-06 
6 . 1 3633 53E-06 6 . 2009919E-06 
5 . 6017010E -06 5 . 5844275E-06 
5 . 7473146E- 06 5 . 7 3617 05E-06 
6 . 2 009919E-06 6 . 2247980E-06 
5 . S844275E-06 5 . 5824495E-06 
5 . 7 861705E-06 5 . 82163 38E-06 
6 . 2247980E-06 6 . 2046267E-06 
5 . 5824495E-06 5.5869516E 06 
1 . 0315S58E- 11 1 . 0459090E- 11 
1. 3693774E-.il 1 . 4720614E-11 
8. 9716894E 12 8 . 974 57 9 OE* 12 
i . 0459090E-11 1 . 0512444E-11 
1 . 4720614E-11 1 . 5492479E-11 
0 . 9745790E-12 9 . 203 0030E- 12 
1 . 0512444E-11 1.0515709E-11 
1 . 5492479E-11 1 . 627 5552E- 1 1 
9 . 203003 OE- 12 9 . 4 6034 98E- 12 
1 . 2381082E-10 1.2254941E 10 
1 . 3749699E-10 1 . 4588913E- 10 
1 . 2835028E-10 1 . 27 82 93 8E- 10 
1 . 2254 94 IE 10 1 . 2080321E-10 
1 .4588913E-10 1 . 5407723E-10 
1 . 2782938E-10 1 . 2876088E- 10 
1 . 2080821E-10 1 . 1883032E-10 
1 . 5407723E-10 1 . 6582781E-10 
1 . 2876088E-10 1 . 2905563E-10 
1 .7970502E-07 1 . 80042 5BE-07 
2 .9995857E-07 3 . 7428040E-07 
1 .4185696E-07 1 . 4 128161E-07 
1 . 8004258E-07 1 . 7648628E- 07 
3.742 8040E-07 4 . 5420162E- 07 
1.4128161E07 1 . 4811355E-07 
1 .7648628E-07 1 .709S664E- 07 
4 .54201 6 2E - 07 5 . 7 1 8 9776E- 07 
1.4811355E-07 1 . 572213 8E - 07 
7.9 861544E-10 8 . 0693 4 81E - 10 
1.4373165E 09 1 . 8106757E-09 
6 . 2276265E-10 6 . 1713 629E- 10 
8.069348 IE- 10 7 . 9701261 E- 10 
1.8106757E 09 2 . 20327 95E- 09 
6 . 1713629E- 10 6 . 4497725E-10 
7.9701261E 10 7 . 7744 679E- 10 


3 .69674 9 8E 03 3 . 6469722E-08 
3 .2603269E -08 3 . 3260385E-08 
3 . 535153 3E 08 3 . 5603954E - 08 
3 . 6469722E- 08 3 . 5309451E-08 
3 . 3260335E 03 3 . 3313600E-08 
3 . 7766963E-06 3 . 8073647E -06 
3 . 3766549E-06 3 . 2 8 044 9 3 E - 06 
3 . 6729106E-06 3 . 6689327E-06 
3 . 9Q73647E-06 3 . 93 03 082E - 0 6 
3 . 2804493E-06 3 . 2578806E-06 
3 . 6689327E-06 3 . 67 4 3 286E-06 
3 . 8308032E-06 3 . 8246800E-06 

3 . 2578806E-06 3 . 3545979E-06 
3 . 6743286S-06 3 . 637 8900E - 06 

4 . 4264530E-08 4 . 4 9 12 680E- 0 8 
3 .93835132 -08 3 . 37 3 57 9 9E - 0 8 
4 . 3050069E -08 4 . 2603 554E - 08 
4 . 491268 CE- 0 3 4 . 5437244E-08 
3 . 873 5799E-08 3 . 82 32 551E- 08 
4 . 2603554E-08 4 . 2413847E-08 
4 . 5437244E-08 4 . 53 877 8 8E - 08 
3 . 823255 IE -08 3 . 894 9 1 1 9 E - 08 
4 . 2413847E -08 4 . 24 94 400E- 08 
1 . 5253305E-08 1 . 48 B6596E -08 
3 . 7765007E- 08 4 . 2 63 2 63 5E - 0 8 
1 . 2780012E-08 1 . 3 653457E-08 
1 . 4886596E-08 1.4851887E 08 
4 .263263 5E-08 4 . 2347613E-08 
1 . 3653457E-08 1 . 44 1 8 4 8 5E - 08 
1.4851887E-08 1 . 5485478E-08 
4 .2347613E-08 3 . 3 5 4 917 2E- 0 8 
1.4418485E 08 1 . 4 9 4 000 OE - 0 8 
3 . 1491441E-04 3 . 1489664E-04 
3 . 144 2 5 64 E - 04 3 . 14 3 9 5 5 6E - 0 4 
3 . 1506562E-04 3 . 1 50 57 1 OE - 04 
3 . 1489664E-04 3 . 14 8 672 9E - 04 
3 . 1439556E-04 3 . 144 6 5 9 IE - 04 
3 . 15O5710E- 04 3 . 1504266E- 04 
3 . 148672 9E - 04 3 . 1 4 8 3 7 0 OE - 04 
3 . 1446591E-04 3 . 1 4 67 01 6E - 04 
3 . 15042 66E- 04 3 . 1502400E-04 
5 . 78617 05E- 06 5 . 82 16 3 3 8E- 06 
6 . 2247980E-06 6 . 2046267E-06 
5 . 58244 95E - 06 5 . 5869516E-06 
5 . 82 163 3 8E - 0 6 5 . 864 87 0 IE - 06 
6 . 2046267E-06 6 . 08 65 67 3 E - 06 
5 . 5869516E- 06 5 . 59 9 5 50 8E - 06 
5 . 8648701E-06 5 . 9257070E-06 
6 . O865673E-06 5 . 93 06 47 6E- 06 
5 . 599 550 8E 0 6 5 . 62 3 9 200E- 06 
1 .0512444E 11 1.0515 7 09E 1 1 
1 . 549247 9E - 11 1 . 627 5 552E- 1 1 
9 . 203003 OE 12 9 . 46 03 4 9 8E- 12 
1 . 0515709E -11 1 . 0654593E-11 
1 .6275552E-11 1 . 57 4 07 53E- 1 1 
9 . 4603498E-12 9 . 7 05 5 53 IE- 12 
1 . 0654593E-11 1 . 1131837E- 11 
1.5740753E 11 1 . 41 5 3 02 2E - 1 1 
9 . 7055531E-12 9 . 93 0 3 6 OOE - 1 2 
1 . 2080821E-10 1 . 1883032E-10 
1 .5407723E-10 1 . 6 5 827 8 IE - 10 
1 . 2876088E-10 1 . 290 5 5 6 3E - 10 
1.188303 2 E-10 1 . 1 8 07 2 8 9E - 10 
1. 65827 81E-10 1 . 6942488E- 10 
1 . 290S563E -10 1 . 2839354E-10 
1. 1807289E -10 1 . 2 0 1 92 09E- 10 
1. 6942488E-10 1 . 6411289E -10 
1 . 2839354E-10 1.2691400E 10 
1. 764862 3S- 07 1 . 7095664E-07 
4 . 542 0162E - 07 5 . 7 1 897 7 6E - 07 
1 . 4811355E-07 1 . 5722138E -07 
1 .7095664E-07 1 . 707 5 8 9 6E - 07 
5 . 718977 6E -07 5 . 67 24 4 0 3E - 07 
1 . 5722133E- 07 1 . 6620917E- 07 
1 .7075896E-07 1 . 8447604E-07 
5 . 67244 0 3E- 07 4 . 5 5 90 03 OE- 07 
1 .6620917E-07 1 . 7 3 17 20 OE - 07 
7 . 970126 IE- 10 7 . 7744679E- 10 
2 . 2032795E- 09 2 . 7677292E - 09 
6 . 4497725E- 10 6 . 83 7 6 7 2 6E- 1 0 
7 .7744679E-10 7 . 828042 OE - 10 
2.7677292E -09 2 . 694 83 57E- 09 
6 . S37672SE-10 7 . 2404 147E- 10 
7 . 3280420E 10 8 . 5455 809E- 10 


3 . 53 0945 IE- 0 S 3 . 39293 08E - 08 
3 . 3813600E- 08 3 . 429173 8E- 08 
3 . 6222 52 3E-QS 3 . 6740 560E - 08 
3 .3 92930 8E 08 3.2220627E 08 

3 . 8308082E- 06 3 . 8246800E 06 
3 . 2578806E -06 3 . 3545979E- 06 
3 . 6743286E -06 3.6878900E 06 
3 . 8246800E-06 3.7690451E 06 
3 . 354S979E- 06 3 . 4 9203 39E - 06 
3 . 637 8900E - 06 3 . 7080213E-06 
3 . 7690451E- 06 3 . 6930851E-06 
3 . 4920339E- 06 3 . 63 11205E- 06 

4 . 54 3 724 4E- 08 4 . 53 8778SE- 08 
3 . 823 2551E-08 3 . 8949119E 08 
4 . 2413847E-03 4 . 24 944O0E- 08 

4 . 53 87788E - 08 4 . 4 586922E - 08 
3 . 8949119E-08 4 . 05 14640E-08 
4 . 2494400E-08 4 . 2789645E-08 
4 . 4586922E-08 4 . 3 546491E - 03 
4 .051464QE-08 4 . 2704330E- 08 

1 . 4851887E-08 1 . 5485478E-08 
4.2347613E 08 3 . 3 549 172E- OB 
1 . 4418485E-08 1 . 4 9 40000E- 08 
1 . 5485478E- 08 1 . 81 40 504E - 08 
3 .3549172E-08 2 . 3 3 21 989E- 08 
1.4940000E 08 1 . 527B16 8E- 08 
1 . 8140504E 08 2 . 14 8443 9E- 08 
2 . 3321989E- 08 1 . 4368399E-08 

3 . 14 86729E 04 3 . 14 S3700E - 04 
3 . 144659 IE- 04 3 . 1 4 67 016E- 04 
3 . 1504266E-04 3 . 1502400E- 04 
3 . 1483700E-04 3 . 147 5565E- 04 
3 . 1467016E- 04 3 . 148S724E - 04 
3. 1502400E-04 3 . 1499689E- 04 
3 . 147 5565E - 04 3.1467365E 04 
3 . 1485724E -04 3 . 1500566E-04 

5 . 8648701E-06 5 . 92 57070E- 06 
6 . 0865673E-06 5 . 93 06476E- 06 
5 . 599550 8E- 06 5 . 62 39200E-06 
5 . 9257070E- 06 5 . 9992466E- 06 
5 . 930647 6E - 06 5 . 7684 54 4E- 06 
5 . 6239200E- 06 5 . 6602 322E- 06 
5 . 9992466E- 06 6 . 07 21S84E- 06 
5.7684544E-06 5 . 6547121E- 06 

1.0654593E-11 1 . 11 31B37E- 11 
1 . 5740753E- 11 1 . 41 53022E- 11 
9. 705553 IE 12 9 . 93 03 600E- 12 
1 . 1131837E-11 1 . 1930260E -11 
1 . 4153022E-11 1 . 1850625E-11 
9 . 9303600E -12 1 . 0133 04 IE- 11 
1 . 1930260E-11 1 . 2833108E-11 
1 . 1850625E-11 9 . 8119282E-12 

1 . 1807289E -10 1 . 20 19209E- 10 
1 . 6942488E-10 1 . 64 11289E- 10 
1 . 2839354E 10 1 . 2691400E - 10 
1.2019209E 10 1.2517224E-10 
1 . 6411289E-10 1. 5083575E-10 
1 .2691400E-10 1. 2515727E- 10 
1 .2517224E-10 1 . 31 32024E- 10 
1 . 5033575E-10 1 . 3571490E-10 

1 . 7075896E-07 1 . 84 47 604 E - 07 
5 . 67 24403E-07 4 . 5 59003 OE -07 
1 . 6620917E -07 1 . 73 17200E-07 
1 . 8447604E-07 2 . 1490 865E - 07 
4 . 559003 OE- 07 2 . 9083002E-07 
1 . 7 317200E- 07 1 . 7721829E- 07 
2 . 149086 5E- 07 2 . S504971E-07 
2 . 9083002E-07 1.75B1928E 07 

7 . 8280420E-10 8 . 54 55 809E- 10 
2 . 6948357E-09 2 . 1 1 52 B27E- 09 
7 . 2404 147E- 10 7 . 57 91 800E- 10 
8 . 54 55809E- 10 1 . 0077 27 6E- 09 
2 . 1152827E 09 1 . 3 072919E - 09 
7 . 57 91 800E- 10 7 . 81 11 14 OE- 10 
1 . 0077276E-09 1 . 2096770E-09 


203 



1 . 4373165E-09 
6.227626SE- 10 

2 . 1 84 62 Q0E- 05 
2 . 1B40660E-05 
2 . 1839086E-05 
2 . 1846130E-05 
2 . 1836863E-05 
2 . 1844319E- 05 

2.184 5767E- C5 
2 . 1833167E-05 
2 . 1846388E-05 

1 . 7 533500E-06 

1 . 7 602275E-06 
1 . 7621491E-06 

1 . 7 53 53 21E - 06 
1 . 7649223E-06 
1 . 7558S44E-06 

1 . 7 539587E- 06 
1 . 7695277E-06 
1 . 7532033E-06 
1.4539800E-13 
1 . 4625191E-13 
1 . 2473472E-13 
1 . 4567049E-13 
1 . 3752474E- 13 
1 . 4545968E-13 
1 . 465S058E- 13 
1 . 2939141E- 13 
1 . 5651731E-13 
7 . 6656300E- 12 
8 . 6580485E-12 
9 . 57 5843 8E- 12 
7 . 7597340E-12 
9 . 3292611E-12 
7 . 8104252E 12 
7 . 83 50000E- 12 
1 . 0014639E-11 
7 . 1 3 5 17 4 8E- 12 
1 . 4670400E- 16 
1 . 593 743 5E- 16 
1 . 6404270E- 16 
1 .4845194E 16 
1 . 6607707E- 16 
1 . 4563494E- 16 
1 . 4969004E-16 
1 . 7246348E-16 
1 . 3777645E- 16 
3 . 5487000E- 13 
3 . 5972269E-13 
3 . 1119822E- 13 
3 . 5725470E- 13 
3 . 4649583E -13 
3 . 3777830E 13 
3 . 5927767E-13 
3 . 3296046E-13 
3 . 497B063E-13 
4 . 506 5100E- 17 
4 . 5572423E 17 
4 . 377 8982E- 17 
4 . 5188662E- 17 
4 . 503 5212E- 17 
4 .4630397E 17 
4 . 5312424E-17 
4 . 44974 91E- 17 
4 . 5262162E- 17 
6 . 8239300E-10 
9 . 4358217E-10 
6 . 2161061E-10 
7 . 2047271E-10 
9. 3642017E-10 
6 . 4334177E-10 
7 . 5922901E- 10 
9 . 2326712E- 10 
6. 5119063E- 10 
7 . 4884200E- 01 
7 . 4884200E-01 
7 . 488420QE-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4384200E-01 
7 . 4884200E-01 
0. OOOOOOOE+OO 


1 . 81067 57 E- 09 
6 . 1713629E 10 
2 .1846130E- 05 

2 . 1 83 6 863 E- 05 
2 . 1844319E 05 

2 . 184 57 67 E - 05 
2.1833167E 05 
2.184S388E 05 
2 . 1845456E- 05 
2. 1823942E 05 
2 . 1846740E 05 
1 .7535321E 06 

1 . 7 64 92 23 E- 06 
1 . 7 558544 E 06 
1 .7539S87E 06 
1 . 7 69 5277E 06 
1.7532033E 06 
1.7543418E 06 
1 . 7747588E 06 
1 .7528148E 06 
1 . 4567049E 13 
1.3752474E 13 
1 . 454 5968E 13 
1 . 4655058E 13 
1 . 293 9141E 13 
1.5651731E 13 
1 . 4 819564E 13 
1 . 1973194E 13 
1 . 5637891E 13 
7 . 7597340E 12 
9 . 3292611E 12 
7 . 8104252E 12 
7.8350000E 12 
1 . 0014639E 11 
7.1351748E 12 
7 . 8765823E 12 
1 .0954885E 11 
7.1126737E 12 
1.4845194E 16 
1 .66Q7707E -16 
1 .4563494E -16 
1 . 4 969004E 16 
1 . 7246348E 16 
1.3777645E 16 
1 . 5039305E 16 
1 . 8030028E -16 
1 .3753474E 16 
3 . 5725470E 13 

3 . 4 649588E 13 
3 . 3777830E 13 
3 . 5927767E 13 
3 . 3296046E 13 
3 .4978063E 13 
3 . 6174416E 13 
3 . 1469997E 13 

3 . 5237724E 13 

4 . 51 88662E 17 
4 . 5035212E 17 
4 .4688397E 17 
4 . 53 12424 E 17 
4 . 4497491E 17 
4 . 52621 62E 17 
4 . 5466088E - 17 
4 . 3760097E 17 
4.S372587E 17 
7 . 2047271E 10 
9 . 3642017E -10 

6 -4334177E 10 
7 . 5922901E 10 
9 .2326712E 10 
6.5119063E 10 
7 . 9880871E 10 
8 . 9538815E 10 
6 .4223399E 10 
7 . 4884200E 01 
7 . 4884200E 01 
7 . 4884200E 01 

7 . 4884200E 01 
7 . 488 4200E 01 
7 . 4884200E 01 
7 . 488 4200E 01 
7 . 488 4200E - 01 
7 . 4884200E 01 
0 . OOOOOOOE+OO 


2.203 2795E- 09 
6 . 4497725E - 10 
2 . 1045767E-O5 
2 . 1833167E-05 
2.18463 88E-05 

2 . 184 54 56E-05 

2.182 3942E-05 

2 . 184 67 4 OE- 0 5 

2 . 184 5528E-05 
2 . 1826232E-05 

2 . 184 6521E-05 
1 . 7539587E-06 
1 .7695277E-06 
1.7532033E-06 
1 .7543413E-06 
1.7747 588E-06 
1 . 7528148E-06 
1 .7542017E-06 
1 .77B0490E-06 
1 .7S30174E-06 
1 .46SS0S8E-13 
1 .2939141E-13 
1 . 5651731E- 13 

1 . 4 819564E- 13 
1 .1973194E-13 
1 . 5637 891E - 13 
1 .5120238E-13 
1 .1136375E-13 
1 . 5301333E-13 
7 . 8350000E-12 
1 .0014639E-11 
7 .1351748E -12 
7 . 8765823E -12 
1.0954885E-11 
7 . 1126737E-12 
7 . 847 8174E-12 
1. 1845243E-11 
7 . 2343980E-12 

1 . 4 969004E-16 
1 .7246348E-16 
1 . 3777645E-16 
1 . 50393 05E-16 
1 . 8030028E-16 
1 .3753474E-16 
1.5028102E-16 
1 . 8701680E-16 
1 . 3925336E-16 
3. 5927767E-13 
3 . 3296046E- 13 
3 .4978063E--13 
3.61744 16E- 13 
3 . 1469997E- 13 
3 . 5237724E-13 
3 . 6590718E-13 
2 . 9769566E- 13 
3 . 5304730E- 13 
4 . 5312424E-17 
4.44 974 91E- 17 
4 . 5262162E-17 
4 . 5466088E-17 
4 . 3760097E-17 
4 . 5372 S87E- 17 
4 . 57 0 5 8 99E - 17 
4 . 3075285E-17 
4 . 5281166E-17 
7 . 5922901E-10 
9 . 2326712E-10 
6 . 5119063E-10 
7 . 9880871E - 10 
8 . 9538815E-10 
6 . 4223399E-10 
8 . 4029218E-10 
8 . 5087584E-10 
6 . 3592993E-10 
7 .4884200E-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 488420QE-01 
7 . 4884200E-01 

7 . 4 884200E-01 
7 . 4884200E-01 

7 . 4 884200E-01 
0. OOOOOOOE+OO 


2 . 7 677292E 09 
6 .8376726E-10 
2 . 1845456E-05 

2.182 8942E-05 
2 . 1846740E-05 

2 . 184 5528E-05 

2 . 182 62 32E-05 

2. 184 6521E- 05 
2 . I845912E-05 
2 . 18243 03E-05 
2 . 1846337E-05 
1 .7543418E-06 
1 .774758SE-06 
1 .7528148E-06 
1 . 7542017E-06 
1 . 7780490E- 0 6 
1 .7530174E-06 
1 .7537326E-06 
1 .7803816E-06 
1 .7532016E-06 
1 .4819564E-13 
1 . 1973194E-13 
1.5637891E-13 
1 . 512023 8E- 13 
1 . 1136875E-13 
1. 5301333E-13 
1 . 549 9828E-13 
1 .0056010E-13 
1. 493B816E-13 
7 . 876 5823E- 12 
1 . 0954 8 85E- 11 
7 . 112 67 37 E - 12 
7 . 847 817 4E- 12 
1. 1845243E-11 
7 . 2343980E-12 
7 . 7807566E- 12 
1 . 3052217E-11 
7 . 3829861E - 12 
1 . 503 93 05E - 16 
1 . 803 002 8E - 16 
1 .3753474E-16 
1 . 5028102E- 16 
1 . 8701680E- 16 
1 . 39253 36E-16 
1 . 4973209E-16 
1 . 9535543E -16 

1 . 4166861E-16 

3 . 6174416E-13 
3 . 14 6 9997E-13 
3 . 5237724E-13 
3 . 6590718E-13 
2 . 9769566E- 13 

3 . 5304730E-13 
3 . 7111213E-13 

2 - 7454122E- 13 
3 . 53 15192E-13 

4 . 546 608 8E-17 
4 . 3760097E-17 
4 . 5372587E-17 
4 . 5705899E- 17 
4 . 3 07 52 85E- 17 
4 . 5281166E - 17 
4 . 599842 4E- 17 
4 .2159348E-17 
4 . S131644E-17 
7 . 9880871E-10 
8.9538815E-10 
6 . 42233 99E- 10 
8 . 4029218E-10 
8 . 5087584E- 10 
6 . 3592993E-10 
8 . 8244201E- 10 
7 . 623 63 4 3 E - 1 0 
6 . 3882588E - 1 0 
7 . 4884200E-01 
7 . 4834200E-01 
7 . 4884200E-01 
7 . 4884200E-01 

7 . 4 88 4200E- 01 
7 . 4884200E 01 

7 . 4 884200E-01 
7 .488 4200E-01 
7 .4884200E-01 
0. OOOOOOOE+OO 


2 . 6948357E - 09 
7 . 2404147E- 10 
2 . 184 5528E- 05 
2 .1826232E-05 
2 . 1846521E-05 
2 . 1845912S-05 
2 . 1S24303E-05 
2 . 1846337E 05 
2 .18457C8E-C5 
2 . 1827146E-05 
2.18463 Q0E-05 
1 .7542017E-06 
1 .77B0490E-Q6 
1 .7530174E-06 
1 .7537326E-06 
1 . 7303816E-06 
1.75320165-06 
1 . 7540143E-06 
1 . 7763499E- 06 
1 .7532878E-0S 
1 .5120239E-13 
1 .11363752-13 
1 . 53013332-13 
1 . 54998232-13 
1 .0056010E-13 
1 . 4938316E-13 
1 .5679555E-13 
9 .9187928E -14 
1 .4664500E-13 
7 . 8 47 8174 E - 12 
1 . 1045243E-11 
7 .2343980E-12 

7 . 7 807566E- 12 
1 . 3052217E-11 
7 .3829861E-12 
7 . 8141748E-12 
1 . 2967442E-11 
7 . 5347583 E - 12 
1.5028102E-16 
1 . 87Q1680E-16 
1 .3925336E-16 
1 . 4973209E-16 
1 . 9535543E -16 
1.4166861E- 16 
1 . 50265 19E -16 
1 . 9407493E-16 
1 . 4431868E- 16 
3 . 6590718E-13 
2 . 9769566E -13 
3 .5304730E-13 
3 . 7 11 12 13E - 13 
2 . 7454122E-13 
3 . 5315192E-13 
3 .7383858E-13 

2 . 7 0697 58E-13 
3 . 5351372E-13 
4 . 5705899E-17 
4 . 3075285E- 17 
4 .5281166E-17 
4 .5998424 E -17 
4 . 21593 4 8E - 17 
4 . 5131644E- 17 
4 . 6160301E-17 
4 . 2146912E- 17 
4 . 5041686E- 17 
8.40292 18E - 10 
8 . 5087584E- 10 
6 . 3592993E- 10 
8 . 8244201E - 10 
7 . 62363 4 3 E-10 
6 - 3882588E- 10 
9.1852370E-10 
6 . 7960463E 10 
6 . 538 3493E-10 
7 . 4894200E-01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4834200E-01 
7 . 4884200E-01 
7 . 48842O0E -01 
7 . 4884200E-01 
7 .4884200E- 01 
7 . 4884200E-01 
0. OOOOOOOE+OO 


2. 1152827E-09 
7. 5791800E-10 
2.1S45912E-05 
2. 1824303E - 05 
2. 1846337E-05 
2. 1845708E-05 
2 . 1827146E-05 
2 . 1846 300E - 05 
2.1343922E-05 
2 . 1 832525E -05 
2 . 1346200E-05 
1 . 7537326E-06 
1 . 7803816E-06 
1.7532016E-06 
1 . 754014 3E-06 
1 . 7768499E-06 
1 . 7532878E-06 
1 . 756207 9E- 06 
1 . 77Q2612E-06 
1.7533500E-06 
1 . 549982 0E-13 
1 . 0056Q10E-13 
1 . 493 881 6E -13 
1 . 567 955 5E -13 
9. 9187923E-14 
1.46645 0 OE - 1 3 

1 . 53 66 882E - 13 
1 . 0731573E -13 

1 . 4 5 3 9 80 OE - 1 3 
7 . 7807566E-12 

1 . 3052217E-11 
7 . 3829861E- 12 
7 . 8141743E-12 
1 . 29 67 44 2E - 1 1 
7 . 5347 583E- 12 
8 . 1030285E-12 
1 . 16 6 1 57 7E - 1 1 
7 . 66563G0E-12 
1 . 4973209E- 16 
1 . 9535543E-16 
1 . 4166861E 16 
1 . 5026519E-16 
1 . 94 07 493E-16 
1 . 4431868E-16 
1. 5352456E-16 
1 . 8293547E-16 
1 . 4 67 04 0 OE 16 
3 . 7111213E-13 

2 . 7454122E-13 
3. 5315192E-13 
3 . 73 83858E-13 
2 . 7069758E-13 

3 . 5 3 513 7 2E - 1 3 
3 . 7008274E-13 
2 . 842565 OE- 1 3 

3 . 54 87 00 OE- 1 3 
4 . 59984 2 4E- 17 
4 . 215 9 3 4 8E- 17 
4 . 5131644E-17 
4 . 61603 01E-17 
4 . 2146912E-17 
4 . 5041686E-17 
4 . 6004914E-17 
4.2769709E-17 
4 . 50 65100E-17 
8 . 8244201E-10 
7 . 6236343E-10 
6. 38825B8E- 10 
9 . 1852370E-10 
6.7960463E-10 
6 . 5383493E-10 
9 . 3977602E-10 
6 . 3015099E-10 

6 . 8239300E-10 
7 . 4884200E-01 
7 . 4884200E -01 
7 . 4884200E-01 
7 . 4834200E-01 
7 . 4884200E- 01 
7 . 4884200E-01 

7 . 4 8 34 200E - 0 1 
7 . 488420 OE - 01 

7 . 4 8 84 20 OE -01 
0. OOOOOOOE + OO 


1 . 3 072919E- 09 

2. 1845708E-05 
2 . 1827 14 6E- 05 
2. 1846300E-05 
2. 1843922E-05 
2. 183252 5E-05 
2 . 1846200E -05 
2 . 1840660E-05 
2 . 1839086E- 05 

1 . 7 540143E-06 
1 . 7768499E-06 
1 . 7 532 878E- 06 
1 . 7562079E - 06 
1 . 7702612E- 06 
1.7 533 500E- 06 
1 . 760227 5E- 06 
1 .7621491E -06 

1. 5679555E-13 
9. 9187 92 8E - 14 
1 . 4664500E-13 
1. 5366882E- 13 
1. 0731573E-13 
1 . 4539800E-13 
1 . 4625191E-13 
1 . 247 347 2E- 13 

7 . 8141748E- 12 
1 .2967442E-11 
7 . 5347 58 3E- 12 
8 - 103028 5E- 12 
1. 1661577E-11 
7. 665630OE 12 
8 . 6580485E-12 
9. 5758438E-12 

1 . 502651 9E - 16 
1 . 9407493E-16 
1. 4431B68E-16 
1. 5352456E-16 
1. 8293547E-16 
1 . 4670400E- 16 
1 . 593743 5E - 1 6 
1 . 6404270E- 1 6 

3 .7383858E-13 
2.7069758E-13 
3 . 5351372E-13 
3 .70Q8274E- 13 
2 . 8425650E 13 
3 . 5487000E-13 
3 . 5972269E- 13 
3 . 1119822E- 13 

4 . 6160301E-17 
4 . 2146912E-17 
4 . 504 1 686E- 17 
4 . 6004 914E- 17 
4 . 2769709E-17 
4. 5065100E-17 
4 . 5572423E- 17 
4. 3778982E 17 

9. 1852370E-10 
6 . 7960463E- 10 
6. 5383493E- 10 
9. 3977602E 10 
6 . 3 015099E- 10 
6 . 8239300E 10 
9.4358217E-10 
6 . 216106 IE - 10 

7 . 4884200E 01 
7 . 4 884 200E- 01 
7 . 4884200E 01 
7 . 4 884 200E- 01 
7 . 4884200E - 01 
7 . 4884200E-01 
7 . 4884200E-01 
7 . 4884200E-01 

0. OOOOOOOE+OO 


7 . 77 663 53E 10 

2 . 1343922E-05 
2 . 183 2525E - 05 
2 . 1846200E-05 

2 . 1840660E-05 
2 . 1B3 9086E • 05 
2 . 1946130E - 05 
2 . 1836863E-05 
2. 1844319E- 05 

1 .7562079E- 06 
1 . 7702612E-06 
1 . 7533 500E-06 
1 . 7602275E-06 
1 .7621491E-06 
1 .7535321E-06 
1 .7649223E-06 
1 .7S58544E-06 

1 . 53 66882E-13 
1 . 0731573E- 13 
1 . 45 3 9800E- 13 
1.4625191E 13 
1.2473472E 13 
1 .4567049E-13 
1 . 3752474E-13 

1 . 4 54596SE-13 

8 . 103028SE-12 
1 . 1661577E-11 
7 . 6656300E - 12 
8 . 6580485E 12 
9 . 575843 BE 12 

7 . 7 597 340E- 12 
9 .3292611E-12 
7 . 81042 52E - 12 

1 . 5352456E 16 
1 . 8293S47E 16 
1 . 4670400E-16 
1 . 59374 3 5E- 16 
1 . 6404270E- 16 
1. 4845194E-16 
1 . 6607707E-16 
1 .4563494E 16 

3 . 7 008274E - 13 
2 . 8425650E- 13 

3 . 5487000E-13 
3 . 5972269E-13 
3 . 1119822E- 13 
3 . 5725470E- 13 
3 . 4649588E 13 
3 . 3777830E 13 

4 . 6004 914E- 17 
4 .2769709E-17 
4 . 5065100E 17 
4 . 5572423E-17 
4 . 3778982E 17 
4 . 5188662E-17 
4 . 5035212E-17 
4 . 4688397E-17 

9 . 3977602E-10 
6 . 3 015099E - 10 
6 . 82393O0E -10 
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HPT 1 - Max power - Mass averaged: 
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HPT1 - Max power - High temperature streamline: 
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1275488 .125 
1283159.944 
1239454.385 

1294101.381 
1296849.155 
1294999.942 
1291706.919 
1287430.792 
1232214.178 
1276026.816 

1269183.12 
1261622 . 65 
1253352.085 
1244372.688 
1234636.376 
1224062.804 
1212436.959 

1199319.479 

1184454 . 479 
1167675.024 

1148599.814 

1127118.716 

1103821 . 814 
1078615.991 
1G49886 .253 
1014461.748 
971575.4713 
928452.9769 
912710.0383 
999538.5689 

1024556.716 
1020755.464 
1007409.083 
994491 . 3791 
984280.3104 


213 .7535496 
228.3462013 
216.0797396 
214 . 4186074 
219.0556362 
220.4491074 
220.1306449 
218.4710914 
218.6341176 
219.2939185 
219.5469849 
218.0146573 
214 . 8767426 
210.4504239 
205.4279798 
201.5639105 
199.6967323 
199.901098 
203.5542228 
208 . 5842319 
213.9890612 
219.7755159 
225.650816 
231 . 5872877 
237 . 5970562 
243.6246136 
249.6141627 
255.5826876 
261.5443532 
267.5323968 
273.7337962 
280.3494671 
287.5356723 
295.4885152 
304 . 1021445 
313.0604302 
322 .4144125 
332.8202322 
346.5393325 
366.4237725 
393.2397397 
416.661042 
366.9117658 
335.9278348 
321.0972363 
315.0504912 
311.2279984 
307.8753852 


0.000974026 
0.000974026 
0.000974026 
0.000974026 
0.000974026 
0.000974026 
0.000974026 
0.000778679 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0 . 000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0.000583333 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 
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0 
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HPT1 - Max power - Low temperature streamline: 


1 80 
1 .17883E-05 
7 . 29015E- 02 
1.48014E-01 
2 . 98833E-02 
7 . 87784E-12 
2 . 61288E- 08 
3 . 74385E-06 
5 . 41438E-08 
6 . 67Q06E-Q9 
3 . 15494E-04 
4 . 99679E' 06 
6 . 03887E- 12 
7.81056E-11 
7 . 87986E-08 
3 . 09774E-10 
2 . 18535E-05 
1 . 74444E-06 
2 . 291 68E- 1 3 
5 . 94824E-12 
1 . 22867E- 16 
3 . 60300E-13 
5 . 5 5086E- 17 
1 .15043E-10 
7 . 4 8 842E-01 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . 00000E+00 
4.36994E-10 
2.44454E -21 
1361.493333 
1361.466192 
1361.604002 
1361.804629 
1362.145143 
1362.594095 
1363.232964 
1364.042053 
1365.162053 
1366.795742 
1369.272409 
1364.526187 
1291.527273 
1277.678082 
1275.978669 
1269.402691 
1257.998873 
1245.793181 
1232.620969 
1218.498533 
1203.778588 
1189.737367 
1177.237046 
1167.504991 
1163.923057 
1168.342227 
1172.833301 
1173-29684 
1167.346347 
1159.752545 
1153.419329 
1146.829021 


29 0 

0. OOOOOE+OO 
0 . OOOOOE+OO 
2 . 3 3000E-01 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
7. 67000E-01 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 

1578546.667 164.839656 

1578425.56 164.5952474 
1579053.754 163.5284478 
1579966.982 162.0800594 
1581521.66 159.5112592 
1583581.799 156.0765876 
1586517.578 150.9349891 
1590234.226 144.1284106 
1595407.78 134.1615632 
1602985.596 118.3095485 
1614682.489 89.43887755 
1584822.24 51.07783472 
1329161.844 315.6673148 
1299682.65 377.7369451 
1326068.804 370.1731122 
1309584.362 387.5991185 
1285392.64 403.5493714 
1258396.257 417.5255673 
1227732.865 429.6271387 
1193735.771 438.8576732 
1157823.321 441.8781725 
1125553.566 434.522286 

1099213.408 417.5474688 
1083412.816 388.1184497 
1089289.441 341.2596273 
1125630.029 290.5231776 
1169261.064 253.0616594 
1204656.195 231.0850281 
1218667.702 219.5558085 
1222097.051 209.8811117 
1226853.325 199.6633722 
1232944.565 188.3923005 


287.03 0 


0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 
0.000974026 0 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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1 
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1143.312013 
1140.402407 
1140.036382 
1175.145197 
1211.306202 
1207.12102 
1206.773496 
1208.699917 
1287.330753 
1287.762684 
1288.820598 
1290.589866 
1293.468976 
1297.547167 
1303.602496 
1164.931084 
1144.175412 
1131.419472 
1123.69347 
1116.675521 
1109.91928 
1102.738785 
1098.641442 
1094.594253 
1083.994215 
1073 . 370783 
1062.381619 
1051.045407 
1039.160949 
1031.110655 
1026.345079 
1028.250405 
1040.345257 
1058.673726 
1078.47177 
1092.282331 
1100.246523 
1104.8363 
1106.83783 
1105.5863 
1105.345968 
1105.586062 
1106.31671 
1106.954414 
1175.93931 
1185.101995 
1180.331047 
1174 . 886129 


1243698.304 
1258885.435 
1231944.708 
1335021.695 
1369922.242 
1337729.59 
1319203.566 
1307136 
1267688.4 
1269355.101 
1273488.41 
1230431.457 
1291761.396 
1307914.898 
1332182.473 
870116.5997 
853904 . 2104 
839969.4267 
834659.5046 
828255.7264 
820438.2077 
810998.8817 
759476.4425 
784471.2937 
761175.8426 
736268.411 
710342.4874 
682847.2009 
654660 . 4653 
631251.7175 
616296.2965 
616917.7426 
643772.7946 
689551 . 6292 
746780.4128 
793869.6333 
828547 . 8406 
855583.6909 
876502 . 0781 
886148.2441 
893562.7729 
913599.3989 
929963.5037 
941075.4732 
1076786.092 
1075679.495 
1046363.434 
1C23109 . 349 


175 . 4889752 
161.0597535 
145.2191088 
90.38037157 
147.4152696 
165.8133053 
175.0043011 
180.816468 
212.7221961 
211.6905576 
206.779503 
197.7692586 
182.5206961 
158.9508226 
118.4879717 
253 . 5501099 
385.1477051 
445.4349092 
481.7317531 
508.2110508 
530.0886496 
545.8099993 
558.8569498 
571.6215983 
585.1416047 
594 . 5327905 
599.7291891 
599.3197894 
595.8591248 
582.3609096 
560.2919372 
525.6422008 
476.6042655 
422.8428595 
375 . 6378799 
337.067325 
306.7501579 
284.8229632 
269.4715436 
257 . 312107 
244 . 0794738 
229.8527779 
214.9652901 
202.9952434 
54 . 5923099 
120 . 8449247 
160 . 4974738 
179 . 9368539 


0.000974026 

0.000974026 

0.000974026 

0.000974026 

0.000974026 

0.000974026 

0.000974026 

0.000778679 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 

0.000583333 
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0 
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0 

0 

0 

0 
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HPT 1 exit to nozzle exit - Max power - Mass Averaged: 


1 39 

1 . 1183 140e-05 
7 . 2 902 4 lOe -02 
1. 4801330e-01 
2 . 9383340e -02 
9 . 6843060e-12 
4 . 4202050e-08 
3 . 4493320e-06 
3 . 7568700e-08 
3 . 8295280e-08 
3 . 1464130e-04 
5 . 9002 840e-06 
1 . 4215190e-ll 
1 . 534832 Oe -10 
4 . 8974110e-07 
2 . 167 872 Oe- 09 
2 . 1842 190e-05 
1 . 7583280e- 06 
9 . 7461800e-14 
1 . 0319860e-ll 
1 , 9406000e-16 
3 . 2799880e-13 
4 . 1740790e-17 
8 . 12 90550e - 10 
7 . 4884 160e- 01 
0 . 0000000e+00 
0. 0000000e+00 
0 . 0000000e+00 
7 . 1511560e-10 
2 . 4445390e-21 
1202 . 61 
1184 . 993301 
1146.54476 
1111 . 539513 
1079 . 720221 
1050 . 839058 
1024 .657707 
1000.94736 
979 .4887209 
960.0720028 
942 . 4969292 
926.5727337 
912.11816 
898.961462 
886 . 9404038 
875 . 9022595 
865.7038134 
856.21136 
847.3007039 
838 . 8571599 
830.7755529 
822 . 9602179 
815.325 
807 . 7932547 
800.2978473 
792.7811535 
785 . 1950591 
777 . 50096 
769 . 6697622 
761 . 6318819 
753 . 5272454 
745.2052892 
736.72496 
728.1047145 
723 . 7504913 
719.3725196 
714.9757813 
710 . 5658524 
706.1489037 


29 

0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOCOE+OC 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . 000QQE+00 
0. OOOOOE+OO 
0.000005+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE + OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOCOE+OO 
980842 . 3 
955524.7629 

831071 . 8 
723225 . 8061 
630506.8605 
551500.6221 
484858.3293 

429296 . 8 
383598 . 4317 
346611.2013 
317248.6653 

294489.9597 

277379. 8 
265028 . 4813 
256611 . 8781 
251371 . 4445 
248614 .2141 

247712 . 8 
248105.3949 
249295.7709 
250853.2797 
252412 . 8525 

253675 
254405 . 8125 

254436 . 9597 
253665 . 6909 
252054 8349 

249632 . 8 
246493 . 5741 
242796 . 7245 
238767.3981 

234696 . 3213 

230939.8 
227919.7197 

226835.616 
226123 . 5453 
225854 6875 

226104 . 3213 
226951 . 824 


260 

264 . 378224 
283 . 547 
303.839296 
324 . 569528 
345.052112 
364 . 601464 
382.532 
398.158136 
410.794288 
419.754872 
424.354304 
423 . 907 
417.727376 
405.129848 
385.428832 
357.938744 
321 . 974 
276 . 849016 
247.8614208 
239.9685248 
232 . 1420512 
224.382 
216.6883712 
209.0611648 
201.5003808 
194.0060192 
186.57808 
179.2165632 
171.9214688 
164.6927968 
157.5305472 
145.4 
148.9824 
176.7896 
230.0912 
315 
437.6288 
604 . 0904 


237.03 


0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0. 02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.01 
0.01 
0.01 
0.01 
0.01 
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Appendix F 


Media Containing Suplemental Data 


Several of the computational results obtained in this research are time dependent and thus are best presented 
with an animation. The attached Compact Disc (CD) contains several figures which were referenced in this 
thesis. The figures were formatted using HTLM which can be viewed form a typical web browser such as 
Netscape Communicator 4.7 or Internet Explorer 5.0. To view these figures simply insert the CD into a 
personal computer, launch the web browser, and open the file “appendix_f.html. ’ 


An index of the figures is given below for convenience: 

4 Validation 

4.4 Convergence Criteria Study 
4.4.3 Results 

Figure F. 1 Convergence time evolution for mass fraction of O 
Figure F.2 Convergence time evolution for mass fraction of OH 
Figure F.3 Convergence time evolution for mass fraction of SO3 

5 Modeling the NASA/DERA Engine Test 

5.4 High Fidelity Modeling 

5.4.2 High Fidelity Modeling Results 

Figure F.4 Static temperature for max power wake model case 

Figure F.5 Mass fraction of SO3 for max power wake model case 

Figure F.6 Mass fraction of S0 3 for max power wake model case (full stage view) 

Figure F.7 Static temperature for non-uniform max power wake model case 
Figure F.8 Mass fraction of S0 3 for non-uniform max power wake model case 

Appendix A Princeton VPFR Validation Modeling 
“Unsteady Non-Reacting Flow Solution” 

Figure F.9 Initial transient 

Figure F. 10 One vortex shedding cycle 

Finally, the files “SMThesis.doc” and “SMThesis.pdf 1 is contains the text of this thesis in Microsoft Word 97 
and Adobe PDF format, respectively. 
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